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I. INTRODUCTION 


This report describes the analytical techniques and computer program 
developed in the fully-coupled rotor vibration study. 

The rotor blade natural frequency and mode shape analysis is implemented in 
a digital computer program designated DF1758. 

The program computes collective, cyclic, and scissor modes for a single 
blade within a specified range of frequency for specified values of rotor 
RPM and collective angle. 

The analysis includes effects of blade twist, eg offset from reference axis, 
and shear center offset from reference axis. Coupled inplane, out-of -plane, 
and torsional vibrations are considered. Normalized displacements, shear 
forces and moments may be printed out and Calcomp plots of natural fre- 
quencies as a function of rotor RPM may be produced. 

The analysis of this report was taken from "Natural Modes of a Helicopter 
Blade", an unpublished report by the late N. 0. Myklestad, Professor of 
Mechanical Engineering at the University of Texas at Arlington. 

This documentation is a revision of the original document by R. J. Brumbaugh 
(BHC Report 299-099-^91). 
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II. ANALYSIS 


A. The Physical Blade System: 

A sketch of the typical blade system and the axis system is shown 
in Figure 1. The blade system is divided into two major parts, 
the hub portion and the blade portion. The hub portion is un- 
affected by collective angle changes. 

Both the hub and blade are divided into a suitable number of sections 
(20 total) which need not be of equal length. A sketch of one section 
and its orientation to the blade axis system is shown in Figure 2. The 
point where two sections meet is referred to as i station. 

Three sets of boundary conditions of the blade system at the axis of 
rotation are considered. The collective mode boundary conditions are 
those of a hinged attachment with axis of rotation about the Y axis; 
the cyclic node, of a hinge with its axis on the X-axis; and the scis- 
sors mode, of rigid attachment. In all three mode types, rotation of 
the blade about the Z axis is constrained by the control system. 
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B. Approximations to the Blade System. 

The following physical properties of each section are approximated as 
being constant over the section: 

p* - Weight per inch of section 

r^,r c - Cross section eg offset from Z axis (Fig. 2) 

S b ,S c ~ Shear center offset from Z axis ^Fig. 2) 

b - Angle of twist between C, B, Z axis system and X, Y, Z axis 
system. (The 9 of a section is calculated as the average of 
the fe f s of the stations at each end.) 


EI^ - Beamwise bending stiffness 

El - Chordwise bending stiffness 
c 

GJ - Torsional stiffness 


2 

- fpb dA (over the cross sectional area) 

I = fpc^dA 
cc 


The control system restraint on blade rotation about the Z axis is 
approximated by a torsional spring constant acting at the origin. 
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C. The Mathematical Model 

1. Calculation of natural frequencies 

One-half of each section's mass and inertial properties is concen- 
trated at each er.d of the section. Each section is then a massless 
elastic structure and each station has one-half the mass and inertial 
properties from the section on each side of it. 

Deflections, moments, and shear forces at the inboard end of a 
section are calculated from the deflections, moments and shear 
forces at the outboard end (see Appendix A). All deflections, 
moments, and shear forces are calculated as coefficients times 
the deflections at the tip of the blade, finally resulting in 
the deflections, moments, and shear forces at the origin. The 
boundary conditions are then calculated as a function of the tip 
deflections leading to the equation: (See Page 43). 


lC (u>)j 


where uu is the frequency of vibration, deflections are as shown in 
Figure 3, and C (jj) is a 5x5 matrix (4x4 if torsion is ignored). 

The natural frequencies are those values of w that satisfy the 
boundary ecuations i.e. for which the determinant of |C| is equal 
to zero. 

The natural frequencies are found by calculating |C| at even increments 
of O) over the frequency range of interest. If the determinant of 

^ as a different sign than the determinant of LC{u)k + then 

a natural frequency is between uu^ and + Aw, and a parabolic itera- 
tion scheme is used to converge to the natural frequency. 

If two natural r 'equencies lie between the sample frequencies then no 
sign change will occur, so if any three consecutive determinants 
have the same sign and the absolute value of the middle determinant 
is the smallest of the three, then smaller frequency increments are 
taken in this range to bracket two roots. 


Ay. . 

J tip 

. 

tip 

Ax . 
tip 

^ tip 

0 . 

t ip 


= {o| 
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2* Boundary Conditions 

The boundary conditions used in calculating the blade natural 
frequencies are as follows: 


Collective 


Out-of -Plane 
Boundary Condition 


{3(0) = 0 


Inp lane 

Boundary Condition 
Q(0) = 0 


Y(0) = 


20*10 VMASS*NB 2 
R*\ OFT " 386.4 * 


M(0) - kgP(O) = 0 


x(0) = 0 


♦(0) * 0 


y(0) = 0 


x(0) = 


20X10 HMASS*NB 2 
R*HS0FT " 286.4 W 


Scissor 


8 ( 0 ) = 0 

y(0) = 0 


x(0) = 0 
*i0) = 0 


The functional notation designates the radial location at which the 
boundary condition is applied. 

The quantities VSOFT, VMASS, HSOFT, HMASS can be obtained from the 
rotor support system shown in Figures 4a and 4b by the following 
equations : 

'i r , i . r\ 6 


VSOFT = 


VMASS 


HSOFT = 


20*10 

R*k 

op 

^HUB , o . p . 
NB 

20 * 1 0 6 


HMASS = 


^HUB. i.p. 


where VSOFT is the out-of-plane restraint elasticity 

VMASS i the effective hub mass in the out-of-plane 
di ction per blade 

HSOFT is the inplane restraint elasticity 
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HMASS is the effective hub mass in the innlane direction 
per blade 

R is the blade radius (in) 

NB is the numbc of blades 

k is the apparent spring rate of the out-of-plane support 
P system (lb^/in) 

M T . n the apparent mass of the out-of-plane support 

HW ’° P system (lb ) 

m 

k. is the apparent spring rate of the inplane support 
system as seen in the rotating coordinate system 
(lb f /in) 

^HUB -d a PP arent m ass o'* the inplane support system 

’ ~ as seen in the rotating coordinate system (lb ) 

k is the olade flapping spring per blade (ft-lb /deg). 
o 1 

The torsional boundary condition for all modes is 

«m) - 

CK 

where CK is the control system stiffness. 


3. Calculation of Mode Shapes 
Assuming that: 

Ay = 1. 

and using 4 (3 if torsion i c ignored) of the boundary condition equa- 
tions, the tip deflections are calculated. Since all deflections, 
moments and forces are already known as functions of the deflections, 
the mode shapes can then be calculated. The mode shapes are then 
normalized to the largest linear deflection or 10 deg. twist. 

Since completely uncoupled modes often result in the above set of 
equations producing a near singular matrix, rode shape calculations 
for uncoupled motifs may be inaccurate. Uncoupled modes result from 
input of untwisted blade at 0 deg. root collective with no eg or 
shear center offset. If the program option to add uncoupled modes 
to plot is used the program actually calculates for a case with very 
small coupling compared to original case (twist angle divided by 
1000, etc.) 
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III. COMPUTER PROGRAM DF1758 


A. Program Description 

The natural frequencies and mode shapes of a blade are found within 
a specified range of frequency and for specified ranges of rotor RPM 
and root collective angle. Rotor RPM and root collective angle are 
taken at even increments between a minimum and a maximum value. Tor- 
sional terms may be calculated or ignored. 

A listing of program DF1758 is in Appendix D. A brief description of 
its function appears at the start of the listing of each subroutine. 
In addition the following system subroutines and Calcomp subroutines 
are called. 

System Subroutines: 

DATE (ndafe) 

ndate - 8 byte variable returned as date 
in form "01/02/71" (2 Jan 71). 


SETIME (tine) 


time - Initializes TIMEX to "time" in minutes. 
TIMEX (tu.dt, tl) 

tu - Time (in min.) since call to SETIME 

dt - Time (in min.) since last call to 
TIMEX OR SETIME 

tl * "time" - "tu" 
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Name- 

Card No. Columns Description Units List 

51-55 Initial frequency in sweep [Default value = . 1*RPMA] (/rev) PSQR 

56-60 Delta frequency in sweep [Default value = 10* max (/rev) DP 

(RPMA, RPMB, RPMC)3 
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C. Output Format 

1. A summary of all input is printed out, 

2. If input requested printout of mode shapes, then one page is 
printed for each natural frequency. Normalized values are 
printed for deflections in the x-y plane, shear forces and moments 
in the beamchord plane, and for torsional displacements and moments 
if torsion is used. 

3. A summary of all natural frequencies is printed. 

4. If input requested a plot, then natural frequency is plotted as 
a function of rotor RPM. The fan plots also show the forcing 
function frequencies as a function of hub type an* number of 
blades. 


The maximum Reflection plane angle shown on output is the arctan of Ay /Ax 
where Ax +Ay is at maximum value. 

See Appendix B for sample of output. 


5. If the input data requests the output to be punched out, the 
following cards are produced: 


CARD DESCRIPTION 

1 Identification Card (Same as input Card 3) 


FORMAT 


2-4 Blade Mass Distribution and Tip Weight 
5-7 Beamwise Mass Moments of Inertia 
8-10 Chordwise Mass Moments of Inertia 
Dp to 6 Collective Blade Modes (13 cards/mode) 

Up to 6 Cyclic Blade Modes (13 cards/mode) 

Up to 6 Scissor Blade Modes (13 cards/mode) 

6 Cyclic Detuning Cards for Collective Modes 
6 Cyclic Detuning Cards for Cyclic Modes 
6 Cyclic Detuning Cards for Scissor Modes 


7F10.0 

7F10.0 

7F10.0 

6F10.0 

6F10.0 

6F1G.0 

6F10.0, 2F5.0 
6F10.0, 2F5.0 
6F10.0, 2F5.0 
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The format for each of the first 10 cards for each blade mode is as follows: 


Column 1-10 

(Out-of-plane)^ 

11-20 

(Inplane)^ 

21-30 

(Torsion)^ 

31-40 

(Out-of -Plane 

41-50 

(Inplane) i+1 

51-60 

(Torsion) 

The 11th card is 


Column 1-10 

(Ovt.-of -Plane V . 

tip 

11-20 

(Inplane) tip 

21-30 

(Torsion). . 

tip 

31-40 

(Natural Frequency/RPM) 

41-50 

(Mode Type Indicator) 

51-60 

(’jamping Factor) 


The damping factor is set to zero for the rigid body flapping mode and 
0.02 for all other modes. 

The 12th card has only one value at present 

Column 1-10 The change in inplane slope across an element with 
zero chordwise stiffness. 

The 13th card will be blank. It may be replaced by the appropriate 
"cyclic detuning" card. 

The first 12 cards contain the following additional information for 
identification purposes: 

Column 61-66 NAME as input on Card 3 

67-68 Mode number of this mode type (1 to 6) 

69-71 Mode type indicator 

1 for collective modes 
-1 for cyclic modes 
0 for scissor modes 
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72-75 Value of root collective pitch 

76- 80 Value of rotor RPM 

The "cyclic detuning" cards at the end of the mode shapes have the fol- 
lowing format: 

Columns 1-10 Natural frequency at low RPM and low pitch angle 
(cpm) 

11-20 Natural frequency at low RPM and high pitch angle 
(cpm) 

21-30 Natural frequency at high RPM and 1 ow pitch angle 
(cpm) 

31-40 Natural frequency at high RPM and high pitch angle 

(cpm) 

41-50 1/2 (high pitch angle - low pitch angle) 

51-60 1/2 (high RPM - low RPM) 

61-65 Average collective pitch 

66-70 Average RPM 

71-76 NAME from input card 3 

77- 78 Mode number 

79-80 Mode type indicator 

2 

The transfer term, pRC is added to chordwise mass moment of inertia, 

I , before the inertia is punched out. 
cc’ 

The mode shape components have the following units 

Component Units 

inplane feet 

out-of-plane feet 

torsion degrees 

If unequal segment lengths are used, the following changes are made to 
the punched output. 

1. The mass and inertias are recomputed to be represented with 
20 equal segments. 

2. Linear interpolation is used to obtain mode shapes corres- 
ponding to 20 equal segments. 
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APPENDIX A 

CALCULATION OF DEFLECTIONS, SHEAR FORCES, AND MOMENTS 
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A. Elastic Coefficients 
1. Definitions: 

dpg = Beamwise deflection of elastic axis due to 

unit force 

d_„ * Chordwise deflection due to unit force 

d._, d w _ = Deflection due to unit moment 
Mb MC, 

v , v„ * Angular deflection due to unit force 
r b ft. 

v mt>> * Angular deflection due to unit moment 

Mb MC» 

W,p ^ Angle of twist due to unit torque 

Wp = Angle of twist due to unit force 

z « Segment length 

2* Beam-Chord Axis Elastic Coefficients 

- v F - 2 2 /2EI 

d p * 7 3 /3EI 

w T a Y/GJ 

3. X-Y-Z Axis Coefficients 

The elastic coefficients are taken from Fig. 10 through Fig. 17. 
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Figure 10. Linear Deflections Due to Unit Torque. 
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Figure 11. Linear Deflection Due to Twist Caused 
by Unit Force. 
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«L = S x«T * d TY 

d LX * ( d PB * d Pc) SM • C0S * - S y d TY 
Ut: d BC * ( d PB ‘ d Fc) SI" * COS# 

d ut * d BC " S y d TY 

d L Y * d pg COS 2 « * d pc SW 2 0 ♦ s x d Tt 

Figure 12. Linear Deflections Due to Unit 
Out-of-Plane Force. 
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W D ‘ - S y«T * d TX 


d nx * d PB SIN2 8 * d PC cos2 *- Vi* 


d DY * d BC “ S 3f S y W T * d LX 


Figure 13. Linear Deflections Due to Inplane Unit Force. 
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%IY S ^LY = U TB C0S ^ + V FC SIn2 ® 

Fig. 14. Angular Deflection Due to Out-of -Plane Unit Force and 
Linear Deflection Due to Out-of-Plane Unit Moment. 
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Figure 15. Angular Deflection Due to Inplane Unit Force and 
Linear Deflection Due to Inplane Unit Moment. 
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*>Mx * ( y MB “ y MC ) S1N 9 C0S ^ 
y MY * ‘'MB 0082 9 + r MC SIm2 ^ 

Figure 16. Angular Deflection Due to Out-of-Plane 
Unit Moment. 
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v * 

C> 

sine 

. T/mc cose 



V * "MB SIf|2 9 * "mc c ® 2 ® 

"qy * ("mb - "mc) sw »a» »* "MX 


Figure 17. Angular Deflection Due to Inplane Unit 
Moment. 
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B. Dynamically Equivalent Mass System 

Figure L8 shows the location of mass* m • , with respect to 
the coordinate system of Figure 1. Thex"-y"-z" axis 
system is parallel to the beam-chord axis system with its 
origin at the c.g. I^, 1^ t and I z _ are second products 

of the mass with respect to the y"-z", x' f -z", and x'-y" 
planes y respectively. 


The masses * m^, my* and and the offsets* a^* ay* and 
a^y form a dynamically equivalent mass system when: 


2m^ + 2my 


+ 2m Z = m i 



m. 



2a z 2m z 
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In Figure 19 all axis systems shown are parallel and the 
dynamically equivalent mass system is the same as shown in 
Figure 18 . Figure 20 shows the equivalent mass system 
dimensioned relative to the x""-y ,,,, -z M " axis system. The 
accelerations of the x M,, -y" M -z MM axis system are, from 
Figure 19: 

2 

-n R in the z direction 

-(n 2 + in the x direction 

2 

- w Ay in the y direction 

If bx» by, and b 2 are the x"', y'", z"' coordinates of a 
mass, m, then the forces produced by that mass are: 


F v - n 


Fy = m 


= m|n 2 ^R ♦ b z ) + - b x j + ^ 2 p(- Y - b„)J 

* m| n 2 (- X + b x j ♦ w 2 *(y + b y j + w 2 *(b z ) + (« 2 A w^xJ 


2 <t> ^X - b x j + w 2 p|b z J + w 2 AY 


Where : 

r x = r c COS 0 + TgSIN 9 
ry s r^SIN 0 - TgCOS 0 
X = r x -*r Y 
Y = r y + * r x 


From Figure 20: 


m x 

± a x COSy 

N, 

± a x SINy 

; a K' 

b J. 

|pSIN 0 -^COSoj 

my 

i aySIN 7 

± ayCOSy 

’ “y! 

JpCOS 0 +OSIN 0 J 

m z 

+ 

a z 

1 pa z 

! ®Z 
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Since Ax, Ay, 6, p, and 4* represent small displacements, 
products of these terms may be ignored. The forces are: 


Produced 

3E 


m y \u 


|u 2 |r t &x (3 SIN6 * a x ‘J'COSej +u. 2 4^|: 
* w 2 p|- Y ± a x SIN ej 
[w 2 ? a y p COS 0 ? a,, + SIN + w 2 -it 


w 2 


+ u, 2 p|- Y + a y COS ejj 

*z !b z[» 2 ( r i a z) + -M*) " w M' 7 )] 
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Summing the forces, subtracting the ntrif ugal force teru, 
2 

n Rm, and making the substitution, m = 2 (n» x + m Y + n^): 


F = 2 

F z - O) 


n (^ r x - P r Y ) 

mr x +4>^l 2 + u» 2 jmr Y ♦ 2 + w 2 1 


f x =-n 2 mr x +<^« 2 + w 

2 2 
= w <l>mr x + u> mAY 


The moments about the x" , -y , "-z" , origin, substitutinp I zz , 


X BB and I CC’ are: 


I cc ( +SINY)/iQ 2 COSv+ w 2 ‘t>SINej+ (+C0SY)(+0, 2 «t> co^ejj 
I BB ( ± COSY)|i f 2 SINYf w 2 *COSej+ ( + SINy)(+u > 2 <fr SINe)J 

I zz o 2 * * “ 2 *) + (**)(* ,2p )] 


i cc ^f3SiNe+^cose^+ w 2 ^oosej + (isiNY)(n 2 + w 2 )|±(isiNe+‘icose) 
i BB ^+pcos 0 + 4 >s iNej |+ w 2 4 »sin e j +( +cos y) (a 2 +u, 2 )(tPcos o+^sin ©J 
i zz [cn>(io 2 p)* <ip>(in 2 )] 

21 

I cc[( ± ps ^ 6 *^ cos ej^cos y ? w 2 <i>sin ej 

5 cosy n 2 + w 2 j ^psiNe?4cosejJ 
I BB[( fPCOS0 ^ SINe ( ±il2sINY ±w 2 ^COS«j 

i siNY(n 2 +w 2 ^pcose **siNejj 

I zz[ (±1) ( ±n2 ‘ ,,iw2 ^) + (^)( ±ft2 )j 
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Sxamming and subtracting static moments: 

r = i cc ( -rf * cos 2 0 + o 2 ^sin 2 0 + co 2 0cos^ 

+ I fiB ^n 2 0 C0S2 e + 0 > 2 <t> cos 2 0 + <u 2 0 SIN 2 ©j 

M* = I cc |a'i 8 SIN 2 0 -n 2 0 SIN 0 COS 0 -a> 2 0 SIN 0 COS 0 + 6 > 2 /?SrN 2 0 j 

+ I BB ^ft 2 /? cos 2 0 + ft 2 0 SIN 0 COS 0 + 6> 2 0 SIN 0 COS 0 +Jpcos 2 o^ 

* ^zp fl - fi20 ) 


Q’ = I, 


i 2 0 cos 2 e - g> 2 a sin 0 cos 


+ I BB p 0 SIN 2 0 + a> 2 /3 SIN 0 COS©) 


:(" 2 ' 4 ) 


T’ =<M£l 2 ^I 

’ = p|n 2 (- 


I BB “ 1 


Cc)COS2^c- 2 (l BB - I cc )j 


M' = /3 


I z2 + ^ + i oc sin2( ’) 


I BB cos 2 0 - I cc SIN 2 © 


+ 6,2 ( X ZZ + I BB 

+ 0|^fi 2 + <u 2 )^ 


X BB " I CC 


)siN 0COS0 


<i’ = ^(l B B I Cc) SINl,C0S9t ‘ , ’M I ZZ* I BB S ™ 2 '’* , CC C0S2( ’) 
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Taking the moments about the x"' -y"' -z"' axis due to F„, F , 
and F^ and ignoring the centrifugal force term, fi^Rm: Y 

T* = JF X * xF y 

«' * -VF Z 
5- = xF z 

T' = mfr Y * * r x )|~ft 2 ry*(Q^ + w 2 ) ♦ r v) j 

* '( r X - * r y)(“ 2 )(4y * * r x ) 

M* = -m ( r y + $r^ j(u> 2 ) ^r x 

Q' = m [ r x ~ < * ,r Y )('*' 2 ^(‘ l ' r Y 



Summing moments and subtracting static terms: 

F x = ®(fi 2 ♦ u> 2 )mr y + (q 2 + w 2 )mAx 

2 2 
F y = w * mr x + w mAy 

T ’ Icc)COS2e. m (r„ 2 - r, 2 )] 

+ “ 2 [ ! bb * I cc * "( r x 2 * - 2 )J| 

M ’ = f 2 {- r zz * 'bb 00 ^ 8 * W 1 " 2 ") 

* “ 2 ( J zz * W 0 * 2 6 * i co sin2 e * mr y 2 ; 

+ *[(Q 2 + «> 2 )|i bb - I cc JsiNe COS 6 - Jn 
<?’ - I cc )SINeCOSe - mr x r y J 

*^(hz + 'w*™** + ^C 0052 ® + ar X S 
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Let: I c s I zz + I BB 

I ST + T 

B L ZZ cc 

X R * I BB + I CC + m ( r X 2 + r Y ) 

"mi = “ hz + ^b 00 * 2 ® + I cc sn,2tf 
= I^COS 2 e + I n SIN 2 e + mr 2 

PfScj O 15 Y 

Vir ( r c - i b ) sd,9C0S9 

M w«, = ( X C - I B ) Sil,9COS8 - mr X r Y 

*V^* I C SIN? 9 + V 0 ^ 9 * "r x 2 

T ««n = ( t bb " x cx:) C0S2 * * m ( r Y r x ) 

Changes in forces and moments due to mass system ignoring 
centrifugal force term, £l2 m R^ become: 

2 2 

L* = a> mr^ 0 + cj mAy 

D f = ^Q 2 + a> 2 jmr,.0 + |ft 2 + <u 2 j mAx 

M ’ ■ ( a \ fi a * “%*„)* 

T ’ ' - 2i r)^ 
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C. Action line of centrifugal force, F 

Figures 21 and 22 show schematically the forces and moments 
acting on a section. F is the centrifugal force and F__. 
and F.pY are moments due to the centrifugal forces. 

sta i 

F sta i = Z a2 "“ 

tip 

F TX i = 2- fl 2 mzr x 

— ^ Q mzTy 


299-099-724 


41 





BELL 

HELICOPTER company 




Kx>rr* R « « wvt» 


D. Recurrence Fomulas 

Referring to Figures 10 through 17, 21, and 22: 

Vl * V ^i‘'QV"( L i" F i<»i)“LV- ( D l' F i*t)*'W 

^iH 3 Ay i n - M iV*«i d QY- ( L i ‘ F 1 ) d LY- ( °i ' F i*i ) d DY" T i d n 

*1*1 = WmX-Vqx" ( L i“ F i (, i)*'LX"( D i" F i*i)*'0X 

“in 3 ix i' M i d m" , i d QX~( L i" F i®i) d IX~( D i' F i^i) d IW"^' 4 in" T i d TX 

L in * L i* L ’ 

“in 3 V 

"in * «i- F t (“i- ) - p TXi(V *in)‘ V * 

«in 3 V F i( te i- ta i +l )- F m(*r*in> * V * «’ 

T in 3 T i” f tx (V f m) ■ F Ty(*i~*in)* T ’ 

V 1 * *i* “piA * W FD*i - M L L i - Vi - Vi 

Introduce the following quantities that do not depend on F: 
d *LY d LY“ ZI/ LY 

d *DY d DY~ Zt/ DY 

d 'my d MY~ Zv m 

d QY d QY~ Z|/ QY 

d 'lX d LX‘ Zl/ LX 

d *DX d DX' ZV DX 

d *MX ^X" ZV VL X 

d ? QX d QX‘ ¥l/ QX 
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Introduce the following quantities that do depend on F, 


d FLY 

= F i d> LY 

d FDY 

* Vdy 

d FLX 

= F d' 

*1 LX 

d FDX 

■ F i d ’DX- 

^FLY 

= F i‘ / LY 

"fdy 

= VdY 

^FDX 

= P l“'DX 

l FLX 

II 

H. 

sf 

V FL = 

- F i W L 

o 

II 

' F i W D 

Introduce: 

6/8 = 

U FLY&i 

6y = 

d FLY^i + 

6iA = 

’'FDX^i + 

8x = 

d FL X^i + 

60 = 

+ V 




d TY T i 


d TX T i 
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The recurrence formulas then become: 

0 in = ^ + 


in 


s Ay . * 8 ' 


►in = *i * 8 * 


Ax i+1 

=Ax^ 

* 8x 



*i+l 

= *. - 

8f6 



L in 

= L. + 
1 

"in 

r xi n" 2 *in * m in" 2iy in 


D i+i 

= D. + 

l 

(b 2 * 

" 2 ^"in r yi»i*in ^° 2 * " 2 

)m i 

M i+1 

= F.8y 
l 1 

+ F TX 5 * +M i +zL i + 

%! 


:n 

*i+l 


(%y * w 2 )^ 


Q i+ 1 = F.6x + F ty ^0 + Q. ♦ zD. ♦ o» 2 M ptnPi+i* « M^frin 

T in ■ F ix sfl * f ty s * * T i * (“ 2i r * ° 2 W^in 

+ (fi 2 + <0 ) m i+i r Yi+l Ax i+l + 0,2 ' n i+l r Xi+l Ay iH 


Since the forces and moments are calculated inb’d of the masses 
they will not be zero at the tip. The tip values are: 

^IP = l ’tip 

n ~ D’ 

U TIP TIP 


M, 


TIP 


TIP * M* 


Q n p s Q’ tip 
T ST* 

TIP TIP 


299 - 099-724 


45 



MEUCOPTER OOMBMW 


APPENDIX B 
SAMPLE PROBLEM 
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iNn i unru_ r ncuucnu ■ -Lrn va i u* j 
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SYM MAX AMPLITUDE 
O VERT PLANE 

A HORIZ PLANE 44.0 FT. OM^y 

+ TORSION -10.0 DEG TWIST 


ROOT COLLECTIVE =10.0 .15.0 .20.0 OEG./ 



2.0 BLADES GIM8ALE0 HUB 


VSOFT = 0 
VMASS= 0 
HS0FT= 0 
HMASS= **. 
RS0FT= 0 


K=1 1.561 


K=3.446 

K=6.3t2 


K= 1.056 
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© VERT PLANE 
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STM MAX AMPLITUDE 
VERT PLANE 
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-10.0 OEG TWIST 
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8EUL. 

H6UOOPTER ooMFWtiv 


0S/360 FORTRAN M 


OPTIONS - NAME* MAIN »0PT*02» LINECNT*60,S IZ F*OOOOk , 

S0URCE,F8CDIC,N0LIST,N00ECK,LCA0,MAP,N 'DlT.IP.XRfF 
COMMON /COMA/ JHUB, Nl, ' LOT ,POWT ,ITCF( 191 »NAMF <2 » .NOI2) ,NP6 
*.CDATE(2),JHU81,080M(10),RCOLL(10),Z(41>,INPUN 
COMMON /COMP/ CK, IRCOL * XRCOt(lO). IBOM, RBOM(IO) * 

* SMZ141), ZB AR( 40 i * EYEBI 1201 , 

* EVEC(120i,SBI40) , SC 140) « VMB(40), VMCI40), VF8(40>, VFC (40) , 

* 0F8I40) t DFCI40), THI41), THE 1401 » WTI40), SMI 4?) , ISJ< ,XQSOM| 200 1 
*,AZ BAR, RPWA,RPMB,RPMC,COLLA, COL! B, COLIC .CHORD 

* *RB 1411 *RC (41) 

COMMON /COMC/ N,IER(7), OFFSET 

COMMON /COMO/ CHAT(5»5) »SOMNAT( 200,3) »IPLNI 200.3) « INODE (200,3) * 

1 MM3,MM4,MM5,CT(41), ST(4l), IB ,IST, 

2 IBS! 10,10,3) ,IBE(10,10,3),ISTS(10,3),ISTE(10,3) 

COMMON/M/ VLXI40), VDX(40), VLYI40) , VMX(40), VQXI40), VMV(40), 

* 0PLX(40) ,DPOX( 40 ) ,DPLY (40) ,OPMX( 40) ,OPOX (40) ,0PMV|40), 

* OF L X ( 4 0 ) ,OFOX(40),OFLV(40) ,VFt. X( 40 ) , VFDX (40) »VFL VI 40) * 

* F| 41 ) » BOMS, 0TXI41), 0TYI41), SX(41), SV(4l), EMRX(*1), 

* EMRV (41 ) , EMp«w(41), EMPB0I41 ) « EMBPP(41), EMPPWI41), EMPP0I41), 

* THHO (41 ) , F TX ( 41 ) , FTY(41), WFL(41), WFD(41), FMPPMI41) 

COMMON /COM!/ DE7.MSZ , ICCOFO , SOM, OVRG 

COMMON /COMT/ FYX(41), EY8( 41 ), EVC(41), VB(41), VCI41), X1MK41), 

* XITI41), FYRI41), EMRBI4I), EMRC(41), f MRR (41 ) , FMPS0I41), CC02 
»,OVPLT,rVLTN,SVLIN 

*, PL*rES,MUPTYP 

REAL *3 SPECIFICATION FOR IBM; COMMENTED FOR COC VERSION. 

REAL *B CMAT 

LOGICAL LOT, OfT , LGP 1 , LGP2, OVRG, CC02, SVLTN 
1, FIRST 
1, OFFSET 

DIMENSION PP ( 200,3) , IN(3), PQI5.1), SQ(S), 00(5), TOMNAT (SO) , 

* SM ZR X ( 4 1 ) , SM7RV (41) 

DIMENSION A A ( 17 ) 

EXTERNAL ABOUMP 

CALL FRRSFT ( 207, 0,0, u, ABOUMP ,0 ) 

CALL ERRSET ( 20« , 0, 0, 0, ABOUMP, 0 ) 

CALL ERRSE T ( 209,0,0,0, A80UMP.0) 

CALL ERRSET ( 251 ♦ 0, 0, 0, ABOUMP, 0) 

TYCE * 0.0 
CVR PS *0.1047 108 
REAO (5,H AA 
WRITE (6,?) 

1 FORMAT (17A4) 

2 FORMAT (1M1) 

00 3 1*1,6 

3 WRITE (6,4) AA 

4 FORMAT (17X, 17A4 //////// ) 

LOT « .FALSE. 

10 CALL INPT(Tv C E, FIRST) 

IF (TYCE .N-. 0.0) GO TO 1000 
MM3 *3 
MM4*4 
MM5 *5 

1 FI .NCT.CC07) GC TO 20 
MM3 *2 

MM* *3 
MM5 *4 
) NOB *1 
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HELICOPTER cjompxumv 


INI 1>*C 
INI 2) *0 
INI3)*0 

c cottecTive angle swee? * 

(;******♦****•%♦♦*****♦***♦ 

00 TOO IST*1»IRC0L 

C ***« ****** ************ »***•*•••*•***••*••**• »••*•»**»** 

C CALCULATE COEFFICIENTS DEPENDENT ON COLLECTIVE ANGLE • 

C ******** ************** ****•••**••«**•••**•*• *•**•••**•• 
ISTSIIST»l)*IN(l)4l 
ISTSIIST»2)*INI2)*1 
ISTSIIST,3)*INI3)*1 
DO PO I=N0B.N1 
ZTH*TH(I) 

IFII.GT.JHU"> ZTM*ZTH4XRC0LII$T> 

ST I 1>=SINIZTH> 

CTI I)=COSIZTH> 

SS ST ( I )**2 
CCT. I.-SSTH 
SCTH*STm*CT(l) 

EMRXI I)=EMRC !1>*CT« 1> ♦EMRBIT>*ST< I) 

EMR Yl I )*EMRC (I ) *ST( I ) -gMRB! I )*CTt I ) 

EM80W (I )*vc f I )*CCTH»vp( l )*SST h 

EM8B0 1 I 1 »( E YB «I)-€YXfT))*CCTM**2*|FVC I I >-F YXfl) >*SSTH**2 
ENBPW(I»«-XINim*SCTH 

ENBPOII>*(EVRU)-PYX(m*CCTMY(EVXI I)-FTC<I))*SSTH 
TF(SVLIN) FMBPOt I> =0. 

Fwppwcn* ypfr)*ccTM» Ycm*ssTH 
Ew»pom*xnm*scTH*»2 

TMNOI 1 ) * IEMR S0( I ) ♦X IM I ( I ) )* fSSTH-CCTM 1-FMRR ( I )*SCTH 
IFII.E0.N1) GO TO 80 
Z TH »T ME 1 1 ) 

IF! I.GT.JHUP ) ZTHsZTH*XRCOL«TST> 

STH*SINIZTW) 

CTH»COSI ZTH) 

SSTH«STH**2 
CCTH*1.— SSTM 
SCTH*STH*CTw 

sxt n»o. 

SV» I>sO. 

IF! SVITN) GO TO 30 
SXt I) "SC II) *CTH*SR| 1 ) *STM 
ST! I) “SCI I ) "STH-SBI I ) *CTh 
30 DTXII)*SVll)*WTm 
DTYII>*SXlI>*NTin 
VLXtI>*SCTH*(VFB< u -vfc II)) 

VOX II )*VFBI T )*SSTM*VFC( I)*CCTM 
VLV1I )*VFB( ! l*CCTH*VFCt I)*SSTH 
vNxm*scTn*fv««B! n-vMCfin 
VOX !I )*VMB ( T )*SSTM»VMCI I)*CCTM 
VNY I >*CCT«»VMC! 1 )*SSTM 

DPI XI I)*SCTH*(DrBU )-OFCI I) )*SXin*OTxm-ZBARII)*VtXIII 
DPOXI 1)«DFB! I)*SSTH*OFC!I )*CCTH*SV( I)*OTXII > -Z8AR 1 1 )*VDX 1 1) 
DPI VI I ) «DFB m*CCTV ♦OFC ! I )* SSTHyS xm *oty m -ZBAR 1 1 » *VL v I n 
opmxii>*vixi i wbap u >*vmxi n 

DPOXI I ) "VOX I I )-ZB AR 1 1 >*VQX( I ) 

DPMVII)"VLYC n-ZBAR m*VMY! I) 


MAI NO PTC 
MAIN05BC 
MAIN059C 
MAINCOOO 

HA1N0M0 

MAIN0620 
NAIN0B30 
MAIN0040 
MAIN0P50 
MAINOfcBO 
MAINOETO 
MAIN06P0 
NATN069C 
MAIN0700 
MA1N0TI0 
MAIN0720 
NA1N0730 
MAI NO 740 
MAIN075C 
MAIN0760 
MAI NO 770 
MAIN07P0 
"AIN0790 
MAINOPOO 
MAIN0P10 
MAIN0"20 
NAIN0P30 
MAI NO 840 
MAIN0850 
MAIN08B0 
MAIN0P70 
MAINCPPO 
MAINCP90 
MAI N0°00 
MAI NO«IO 
MAlNC°?0 
MAIN0930 
MAI NO 940 
MAIN09*0 
MAIN09B0 
MAIN0970 
MAIN0OP0 
MAIN0990 
MAIN1000 
MAI Nl 01 0 
MAIN! 020 
MAIN1030 
MAI Nl 040 
MAIN! 050 
MAIN1060 
MA INI 070 
M AT Nl OP 0 
MA INI 090 
MA INI 100 
MAIN1110 
MAIN! 120 
MAIN1130 
MAIM 140 
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80 CONTINUE 

SMZRXIN 1 )* 0 . 

SHZ RY INI )* 0 . 

J«N 1 

DO B 5 1*1 *N 

KcJ 

J*J -1 

$M 2 RXtJ>*SMZRXm*EMRXlK»*ZIKI 
85 SN 2 RY IJ )*SNZRY(K ) ♦EMRV(K)*Z (K| 

N 08 *JHUB *1 

C ROTOR RPM SWEEP * 

C******** ********** 

DO 710 IB* 1 tIBOH 

C **************** ******** ************************ 

C CALCULATE COEFFICIENTS DEPENDENT ON ROTOR RPM * 

C ************************************** ********** 

FIR ST*. TRUE. 

BOMS*RBOM( IB )**2 
tBSnST*IP t l)*IN(l)*l 
IBS (IST t IB,?)*IN( 2 )*l 
IBStIST t ie* 3 >=IN( 3 )*l 
DO 110 1 * 1 »N 
F H )*BOMS*S«Z(I) 

FTXII )*P 0 MS*SM 7 RXfI ) 

FTYin*eOMS*SMZRV(i > 

WFL ( I )*F ( I ) *DTY ( I ) 

WFDII >«Fm*OTxm 

DPI XI I) *F ( I) *DPLX II ) 

DFDXI I) *F II) *DPDX II)- ZB AR II) 

DFL YII) *F (I) *DPL V < I WB AR 1 1 1 
vFLxm*Fin*vixm 
vFoxm*Fm*vDxtn 
110 VFL Yl I ) *F ( 1 1 *VL V ( I ) 

DET * .TRUE . 

M$Z*MM 5 

C ****** ******************* 

C CALCULATE DETERMINANTS * 

c ************************ * 

CALL COEF ( 1 , it. FALSE. tISOM f XOSOM f PP) 

C**** ****************** ******************* 

C MODE (COLL ECTIVE .CYCL TCt SC ISSORS ) SWEEP ♦ 

C **************************** ************* 

DO 320 1 * 1 t 3 

IF IOFFSE T . AND. I . EO . 2 > GO TO 320 
C ******************************** **** 

C CHECK FOR DETERMINANT SIGN CHANGE • 

C******** ************** ************** 

IGP 1 *PP 1 1 . 1 > .GT.O. 

DO 120 J* 2 « I SCM 
LGP 2 *LGP 1 

LGP 1 «PPC JtD.GT.O. 

IFI (LGP 1 .AND .LGP? ) .OR • .NOT • (LGP 1 . 0 R »LGP ?) ) GO TO 120 
CALL ITER(ItXOSOM(j-l),XOSOM|J) t PP| J-ltl) *PP|J t TM 
IF! .NOT.OVRG) On TO 120 
ini n«iNii)4i 
INI * I N 1 1 ) 

SOMNATIINI tD'SOM 


MAIN115C 
NAIM1AC 
MAI N1 170 
MAIN11RC 
MAIM 190 
MAIN1 200 
MAIN! 210 
MAIN122C 
MAIN1 230 
MAIM 240 
MAIN12S0 

-MAIM 260 
MAIN1270 
MAIN1280 
MAIN1 290 
MA1N1 300 
MAI N1 310 
MAIN1320 
MAIN1*?0 
MAIN1340 
MAIN1350 
MAIN1 360 
MAIN1370 
MAIN13B0 
MA INI 390 
MAI N1 400 
MAIN1M0 
MAIN1420 
MAIN1430 
MAIN1440 
MAIN1450 
MAIN1460 
MAINi470 
MAIN14B0 
MAIN1490 
MAIN1ROO 
MAIN1510 
MA INI 520 
MAIM 530 
MAIM 540 
NAIM *50 

>-MAINl*60 

MAIM *70 
MAINISBO 
MAIM *90 
MA INI 600 
MAIN1610 
MAIN1620 
MAIM 630 
MAIM 640 
MAIN1650 
MAIN1660 
MAIM670 
MAIN16B0 
MAIN1690 
MA INI 700 
MAIN1710 
MAI Nl 720 
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120 

CONTINUE 

MAIN1730 

C 

GO TO 320 

MATN17.0 

C CHECK FOR HISSED ROOTS * 

MAIN1750 

C*«***»4n|t««**m*#*«*f«**t 

MAIM 760 


KNfT *0 

MAIN1770 


IFUBSItST.IBtlKGT.INmiGO TO 320 

MAIN17P0 


S«AN*.5*(X0$0Ma )*XQS0M( ISOM I ) 

MAIM 7®0 


Jl*IB$f 1ST * I B*I > 

MAIM «00 


J2»INII) 

mainioio 


00 130 K*1 * I SOM 

MAIM 8?0 


00 130 J*J1*J2 

MAIM630 

130 

PPfK f Il*PP(K*ll/ABS<XOSOMfK> -SOMNAT U *ni*ABSISWEAN-SOHNATt J, 

IUMAIM8.0 


Kl«l 

MAINI 650 


K?* ISOM-1 

MAINI 860 

140 

PF«S1GN( l.«PP(Kl»ll) 

MAINI 670 


Kl*Ki*l 

MAINI 860 


IPIK1.GT.K2) GO TO 270 

MAINI 690 


IFf PF9PP (K 1 * I ) • IT .0 . ) GO TO 140 

MAIN1900 


DO 265 J*K1*K? 

MAIN1910 


IP|Pf*PPt J*1 ? I).GT.0. ) GO TO 150 

MAIM 920 


K1*J*1 

MAIN1930 


GO TO 140 

MAINI 9.0 

150 

IFf ABSfPPf J*I)KGT.AB$fPPfJ-l*I 11.0R.A8Sf PPf J»1 1) .GT.ABSfPP • *!♦! 

• 1 IMA1M950 


• 11 GO TO 265 

MAIM 960 


tl*J-2 

MAIM970 


DO 160 K*3* 5 

MAINI *80 


li*u*i 

MAINI 990 


PQf K « 1) *PP f L 1 * I 1 •PP 

MAIN2000 

160 

SOf K)»X0S0M(L1) 

MAIN2010 


ICK *© 

MAIN2020 

170 

00 180 K*1 * 2 

MAIN2030 

180 

SOf K)«.5*<S0fK42>*S0fK*3) 1 

MAIN2040 


CAtt COEFCI, I*.FAtSE.*2*SQ,PQI 

MAIN2050 


00 190 1*1*2 

MAIN2060 


QQf U*P0<U1!*PF 

MAIN2070 


00 190 K*J1 * J2 

MA1N2080 

190 

00(0*00(1 )/ Aft SC SOf U-SOMNAT(K*I) 1*A8Sf SMEAN-SOMNATf K *! 1 » 

MAIN2090 


00 200 K*l * 2 

MAIN2100 


IFfOOfKt .tT.0.1 GO TO 230 

MAIN2110 

200 

CONTINUE 

MAIN2120 


ICK«ICK4l 

MAIN213C 


IFf ICK. tE. 101 GO TP 205 

MAIN21.0 


WRITE t6 *91 X ) 

MAIN2190 

911 

FORMAT f 34H CONVERGENCE FA ItURE-tOCAt MINIMUM 1 

MA1N2160 


GO TO 265 

MAIN2170 

205 

IPtOOm.LT. 0016 * 1 ) ) GO TO 210 

MAIN2180 


IPC0Cf2).tT.P0(4 f HI GO TO 220 

MAIN? 190 


POf 3 m1)«OOC1> 

MAIN? 200 


POf 5,1>*00<2> 

MAIN? 210 


SOf 3 1 *$0(1 > 

MAJN2220 


$0(5>*S0<2> 

MA I N? ?30 


GO TO 170 

MAIN??A0 

210 

POf 5* 1 ) *POf 4*1) 

MAIN? 250 


POf 4* 1) »00< 1 ) 

MAIN??60 


SOf 5)*S0f4) 

MAIN? 270 


SOf 4) *SOf 1 ) 

MAIN? ?80 


GO TO 170 

MAIN? 290 

220 

POf 3*1) *POf 4 *11 

MAIN? 300 
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PQt4 t l)*0Qf?> 

$0( 3) *SQ(4) 

$Q<4)*SQ!2> 

GO TO 170 
230 00 245 M*3,5 

DO 240 L*Jl*U2 

240 POl Ht 1) *P0t H *1 )*( $0 <M )— SOMNAT IL*I))/t SMEAN— SOMNATlL *T ) I 
245 POX Mf DESIGN (P0fM«l ! f PF| 

Li*2*K 

L2*3*K 

L3*0 

CALL lTER(I«S0(U>»SQ(K),P0(Lltl)«P0(K v 11 > 

250 1H.N0T.0VRG) GO TO 260 
KNT *K NT ♦ 1 
TOMNAT(KNT)«SOM 
260 CONTINUE 

IML3.NE.O) GO TO 265 

CALL ITER<I*SC(K)*SO(L2)fPO(K*lltPOIL2tl) > 

L3»l 

GO TO 250 
265 CONTINUE 

270 IHKNT.eO.O) GO TO 310 
Jl*IN<n*KNT 
IN( I )*J1 
J3«KNT 

00 300 J*ltKNT 

280 lHJ2»LT»18$(lSTf]8tI KOR.TOMNAH J3 ) *GT »$OMNATf J? f I ) > GO TP 290 
SOMNATf J1 ♦! ) cSOMNAT ( J2v I ) 

J2* J2 -1 
Jl*Ji-l 
GO TO 280 

290 SONNATUltl >*T0*NATU3) 

J3* J3-1 
300 J1«J1-1 

310 CONTINUE 
C 

320 IBEUST,I*,1>*IMI> 

C 

OET «.EAl$E. 

MSZ*HM4 

C 

710 CALL AMPLTD 
C 

ISTElIST f 3)*IN<3) 

IS^H 1ST *1 ) * IN ( 1 ) 

700 I ST EX 1ST » 2 ) * IN ( ? ) 

C*tM ****** ****** ************ ** ** ****** 

C CALCULATE AND PRINT OUT MPOE SHAPES * 

£************************************** 

CALL $U*PY 

£******** Mt**M**f«* 

C PLOT NATURAL PRFO. VS ROTOR RPM • 

IEt LOT) CAIL RLOUT 
GO TO 10 
1000 STOP 
ENO 


MAIN2310 
MA1N2320 
MAXN2330 
MAIN? 340 
MAIN2350 
HA1N2360 
MAIN2370 
*AIN2380 
NA1N2390 
MAIN2400 
MAIN2410 
MAIN2420 
MAIN2430 
MAIN2440 
NAIN2450 
MAIN? 460 
MAIN2470 
MAIN2480 
MAI N2 490 
MAIN2500 
MAIN2510 
MAIN2520 
HAIN253C 
MAIN?*40 
HAIN2550 
MAIN2560 
MAIN2570 
MAIN? 580 
MAIN?590 
MAIN2600 
MAIN2610 
MAIN2620 
MAIN2630 
MA1N2640 
MAIN2650 
MAIN2660 
MAIN2670 
MAIN26P0 
MAIN2690 
MAIN? 700 
MAIN27I0 
MAIN2720 
MAINP730 
MAIN2740 
NAIN2750 
MAIN2760 
MAIN2770 
MAIN27P0 
HA IN? 790 
MAIN? 800 
MAIN2P10 
MAIN?P?0 
MAIN2P30 
MAIN? P40 
MAIN?850 
MAIN2860 
MAIN2P70 
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OS/ 360 FORTRAN H 


OPTIONS - NAM** MAIN,0PT«02 t LINECNT*60,SI2E*0000K t 

SOUR CEtCBCOIC, NOt 1ST t N00FCK ♦LOAD f MAP f NOEDIT f IO,XRf F 
BLOCK DATA 

COMMON /COMA/ JHUB, Nl« LOT , POUT ,1 U F f 19) f NAME ( 2 I f N0 1 2 

* t C0ATE(2) ♦JHUPl ,DBOM( 10),RC0tm0)*7f4l> t INPUN 
COMMON /COMC/ N ,1 EP (7 ) ^OFFSET 

COMMON /COMf / *B(205) ♦ ZX(205) t 20(205), 2L(205), ZS(?05), 

♦ ♦ZM (205 ) ♦ 20(205), ZH(205>, ZT(205> 

COMMON /COMTP/ DEG(200,3>,PLNE(2 t 4) f 00E$(2,5) 

LOGICAL IER, OFFSET 

REAL *8 SPECIFICATION FOR IBM; COMMENTED FOR COC VERSION 
REALMS PINE, ODES 

REALMS ZB,2S,ZX,ZY f ZH»ZM,ZO,7L,7D,ZT 


BLKD0010 
BLK00020 
)*NPG BLKP0030 
BIK00040 
BLK 00050 
ZV ( 205 )BLKP0060 
8LK00070 
BLKOOOBO 
B IK DO 09 0 
• BLKP010C 

BtKDOll 0 
BLK00120 


REALMS 

BEBE 

,BFP S 

♦ BEX 

♦ BEPH 

♦BEY 

BLKD0130 


, X8E 

, XPS 

, XX 

♦ XPH 

♦ XY 

BLKD0140 


, OBE 

, OPS 

* OX 

• OPH 

♦ OY 

BLKD0150 


»ELBE 

,ELPS 

♦ ELX 

♦ ELPH 

♦ELY 

BLKPO 160 


• PSBE 

,PSPK 

,PSX 

♦ PSPH 

♦ PSY 

BLKP0170 


, YBE 

, YPS 

♦ YX 

♦ YPH 

♦ YY 

BLKD0180 


tBMBE 

,BMPS 

♦ BMX 

♦ BMPH 

♦ BMY 

BLKD0190 


,DE8E 

,DEPS 

♦ DEX 

♦ DEPH 

♦OE V 

BLK00200 


,PHBE 

,PHPS 

f PHX 

♦ PHPH 

♦ PHY 

BLKD0210 


, TBE 

, TPS 

♦ TX 

♦ TPH 

♦ TY 

BLKD0220 

DIMENSION 


BEBE (41) , 

BEPS(41) 

♦BEY(41> t 

BEPH (41 ) 

♦BEX(4l)BLKD0??0 

♦ 


,P$BE (41 ) , 

PSPS (41) 

♦ P SY ( 4 1 ) # 

PSPH(41) 

♦PSX(4l)PLK00240 


EQUIVALENCE 


EQUIVALENCE 


, VBE ( 41 ) 
,XBE(41 ) 

, P W BE (41 ) 
, BP BE (4 1 ) 
, QBE ( 41 ) 

, FLBE (41 ) 
,0f BE(41) 
t TB F ( 41 ) 

E (ZB ( 1 ) , BEBE 1 1 ) ) 
( ZB ( 1 24 ) , BE PM ( 1 
(ZXm v XBE(l) ) 
( Z X ( 1 24 ) , XPH ( 1 
(ZO(l)v QBE ( 1 ) ) 
'Z0(124>, 0PH(1 
( ZL ( 1 ) ♦ ELBE ( 1 ) ) 
(ZL(124),ELPH(1)) 
(ZS(1 ), PSBEM )> 
(ZS(124),PS» (1 

E (ZV ( 1 ) , YBH1H 
(?vn?4), YPH ( i 

(ZM(1) V B*BEI1 )) 
( ZM ( 1 24 ) , BMPH ( 1 
( ZD ( 1 ) , DE BE ( 1 ) ) 
(ZD ( 1 24 ) » DE PH ( 1 
(ZH(1),PHBE(1 >) 
( ZM ( 1 24 ) » PHPH ( 1 
( Z Y ( l ) , TB E ( 1 ) ) ♦ 


,YP$(4l ) » YY (41 

,XPS(41> ? XY ( 4 1 

»PHPS ( 41 ) , PHY( 4 
tBMPS(4l) ,BMY(4 
tOPSI 41 ) t OY (41 
tELPS (41 ) » EL Y( 4 
tDEPS (41 ) ,DE Y( 4 
,TPS(41) t TY (41 
( ZB (4? ) »BEPS( 1 ) ) 
) , ( ZB ( 165 ) ♦ BE Y( 1 
(ZX (42) ♦ XPS( 1) ) 
) , ( 7X ( 165 ) ♦ XY( 1 
( ZQ (42 ) t OPS( 1) > 
>,(20(165), 0Y(1 
( ZL (42) fELPS(l) > 
( ZL (165 ) , EL Y( 1 ) ) 
(ZS(42),PSPS(l>) 
) , ( 2$ ( 1 65 ) , PSY( 1 
( Z Y (42 ) , YPS(D) 
> , ( ?Y ( 165 ) , YY( 1 
( ZM (4? > ,BMP S( 1 ) ) 
),(ZM(165)*BMY( 1 
( Z0(42) ,DEPS( 1) > 
>,(ZD(165>,DEY(1 
(ZH(42) , PHP S ( 1) ) 
) , f ZH ( 1 65 ) t PHY ( 1 
ZT (42 ) , TP S ( 1 ) ) 
) , ( 2T (165 ) t TY(! 


(ZT ( 1 24 ) , TPH ( 1 ) ) , ( ZT ( 165 ) , TY( 1 )) 
DATA CDATE/4H0?/2,4H5/75/ 

DATA OFFSET /.FALSE./ 

DATA IER/7*. FALSE./ 

DATA BEBE(l), BEPS(l), BEYfl), BEPH(l), BEX(l), 


, YX (41) ,YP 
, XX (41) ,XP 
) , PHPH (41) « PH 

) f BMPH (41 ) ,BM 
,0X (41) ,0P 

> ,ELPH(41) « FL 

> , OE PH (41) t DC 

, TX (41 ) , TP 

(ZB( 83 ) ,BFX ( 1) ), 
), 

( ZX( 83) , XX( 1 )) , 

) , 

(Z0( 83 ) , QX ( 1 ) ) , 
) , 

( ZL ( 83 ) ,ELX ( 1 ) ), 

(ZS(83),PSX( 1) ), 
) 

(Z Y( 83 ) f YX( 1) ), 
), 

( ZM ( 83 ) , BMX ( 1 ) ), 
), 

(Z0(R3),DSX(i; ), 
)• 

(2H( 83) ,PHX ( 1) ), 
) , 

( Z T( 83 ) , TX ( 1) ), 
) 


, vpH( 41 )BLK DO 250 
« XPH( 41 ) BLK DO 260 
, PHX (4l)BLK00270 
,BMX(41)BLKD0?80 
,0PH(4l)BLKD0290 
*FLX(41)BLKD0300 
,DEX(41)BLK00310 
,TPH|41 )BLKD0320 
)), BLK 0033 0 

BLKD0340 
), BLKD0350 

BLKD0360 
)>, BLKD0370 

BLKD0380 
)), BLKD0390 

etKD0400 
)), BLKD041C 

BLKD0420 
)), BLKD0430 

BLK DO 440 
)), BLKD0450 

RLKD0460 
))♦ BLKD0470 

8LKD0480 
)), BLKD0490 

BLKD0500 
)), BLKD0510 

8LKD0520 
BLKD0530 
BLK 00540 
8LKD0550 
BLK 00560 
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PSBEIl) , 

p$psm. 

PSYU)t 

PSPHIDt 

PSXtl), 


BLKD0570 


VBEI1I, 

YPSf I>, 

YYf 1), 

YPH ( 1 ) , 

YX (l)t 


BLKD05P0 


XBEIll , 

xpsm* 

X Y( 1 ) t 

XPHf l)t 

xxm. 


BLKD0590 


PHBE(1 ) * 

PMPSf 1 ) , 

PHY ( 1 ) , 

PHPM f 1 ) t 

phx m, 


BLKD0600 


B»BEU), 

B*pst i) * 

BMY( 1) * 

BHPM(X)t 

B«Xtl>, 


RLKD0610 


OBtm , 

opst n f 

0Y< 1 ) * 

OPHfUt 

OX(l) t 


BLKD0620 


EtBE(l) . 

ELPSt 1) , 

EL Y( 11* 

ELPH ( 1 ) , 

ELXUlt 


B LX DO 630 


oesem. 

DEP$m, 

DEVf !), 

OEPMa ) f 

DEXf IN 


BLKD0640 


TBE(l), 

Tpsm, 

TY( l) f 

TPH(n f 

Txm 


BLKD06^0 

* / 

l.DO, 

O.DO , 

0,D0» 

O.DOt 

O.DO, 


buKD0660 


O.DO « 

1 * DO t 

O.DOt 

O.DOt 

O.DOt 


BLXD067C 


O.DO, 

0.00, 

l.DO, 

O.DOt 

0.00, 


BLKD0680 


0.00. 

0*00* 

O.DOt 

0. DO* 

l.DO, 


8LKD0690 


O.DO, 

O.DO, 

O.DOt 

l.DOt 

O.DO, 25*0. DO/ 


BLKD0700 

DATA 

PLNE/BHVERT PLA.BHNE 

,8MH0R JZ 

PL,BHANE ,8H 

TORSION, 

BLKD0710 

1 

3*8H 

/» 

ODES 

/8MRIGI J 

BO ,8 HOY ,8M 

1 NOD, 

P LX DO 7?0 

2 

3 

END 

8ME 

»8M ? NOD, BUBS 

/ 

,BM3 OR M0R,8Hf 

NODES , 

BLXD073C 

BLXD0740 

BLKD0750 
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3PTI0NS - NAME * HA IN f 0PT«0? f LINFCNT*60 f SI7F*OOOOK f 

SOURCF f FP CO IC ,*01 1 ST f NODFCK - LOAD » HAP , NOFD I T , I D* XR FF 
SUBROUTINE AMPtTD AHPL0010 

£******************** *♦*****♦******##**##****♦#***♦*#**♦*** AHPL00?0 

C THIS SUBROUT tNF CALCULATES AND PRINTS OU* MODE SHAPFS * AHPL0030 

C********** **%•'**♦******** ************************ ********** AH»L0040 

COMMON /COHA/ JHUBf Nl, LOT t POUT , I TL F ' 19 ) , NAME i 2 ) « NO < ? ) t NPG AMPL0050 

*,CDATE(2>,JHU91,DB0M(10>,RC0LL( 10>,Z(4i),INPW AMPL0060 

COMMON /COMB/ CK, IRCOc, XRCOL(IO), IBOM, RPHMUO) , AMPL0070 

• SMZ( 41 ) » Z B AR ( < -0 ) » EYEB( 120) , AMPlOOPO 

* EYEC(120),Se(40) , SC (40) t VMB(40), VMC(40), VFB(40), VFC(40), AMPj.0090 

• DFB(40), OF C (40) , TH(41), THE(40), WT(40), SM( 4/ ) , ISOM ,XOSOM 1 200 ) AMPL0100 

*,AZEtAR,RPMA,RPMB,RPMC,COLLA,COLLB,COLLC, CHORD AMPL0110 

• »RB HI) »RC(41) AMPL0120 

COMMON /COMD/ CMAT(5,5),S0NNAT(200,3),IP» N( ?00, 3 ) .INODE ( 200 ,3) , AMPL0130 

1 MM3,MM4,MM5,CT(41 >, ST(41), IB ,!ST, AMPL0140 

2 IB S( 10, 10,3) » IBE (10*10*3 ) »I STS (10*3) »JSTE(10»3) AMPL0150 

COMMON /COME/ ZBtZOM, ZX(20?(, ZG(?05), ZL ( 205 ) * ZS(205)» ZYI205 )AMPLD160 

*»ZH(205) t ZD (203) , ZHC205), ZT(205) AMPL0170 

COMMON /HINGES/ LCH ,L CHP1 ,LFH,LFHP1 ,CHOFF, FHPFF « FCH, FFH, 1PUNCT AMPLOIB^ 

* ,RPMPUN,COLPUN AMPL019C 

**LPH» (.RHP1 , PMOFF, FPH AMPU0200 

*,BOMM,TWSM , LDVN5 AMPL0210 

•»LOTS AMPlO?2f> 

*»B0M1 ,TWSI,DELBOM,DFLTWS AMPL0230 

COMMON /COMI/ OET.mSZ, IGCOFD , SOM, uVRG AMPL024C 

COMMON /COMT/ EYXU1). EYP(4l), EYC(41), YB(41), YC(41>. XI Ml ( 4 1 ) , AMP 10/50 

* XI 7(41), E YR (41 ) , FMRB (41) , EMRC(41), E MRR ( 4) ) , EMRS0(4l), CC02 AMP LC 260 

*,0YPL7,0VLIN ; rvtIN AMPL0270 

*, BLADES.HUBTYP AMPLO2B0 

COMMON /CO* 1 TP/ DEG(200,3),PLNF(2,4),DDES(2,5) AMPL02«0 

LOGICAL DET AMPL0300 

LOGICAL CONV , AMI, AM ?, POUT, C<-02 AMPL0310 

•tRBTEST AMPL0320 

♦.LOTS AMPLO j30 

•,L0YN5 AMPL0340 

6,0VPLT,0VLIN,SVLIN AMPL0350 

C REAL *8 SPECIFICATION FOR IBM; COMMENTED FOR COC VERSION, AMPL03A0 

REAL*B CMAT,DMAT,VEC,PVEC,DA,A,DUM AMP 0370 

*, PLNE , ODES AMPL03P0 

REAL*8 ZB,ZS,ZX,ZY,ZH,7M,ZO,ZL,ZO,2T AmP!0390 


RFAL*8 

BEBE ,BEP S 

tBFX 

♦6FPN tBFY 

• 400 


, XBE , XPS 

p XX 

* XPM t XY 



, OBE , OPS 

t OX 

* QPH f OY 

AM 2420 


.ELBE ,EIPS 

vFLX 

tFLPH tFl Y 

AMP 1043 0 


,PSBE ,PSPS 

,PSX 

,P«PH fPSY 

AMr^04/-0 


, YBE , YPS 

t VX 

9 YPW t y> 

AMPL04<;0 


,8MBE ,BMPS 

tBHX 

♦BHPM f CHV 

AMPL04<0 


,DEBE ,OEPS 

tDFX 

,DFP* f DF v 

AMPL04T0 


,PHBE ,PHP S 

t PHX 

*PMPM f Pf- * 

AMPL04P0 


, TBE , TPS 

* TX 

» TPM T « 

AMPL04S 0 

DIHFNSION 

BE BE f 41 ) , 

BFPSftn 

,BFY(41» . fi F PM ( 4 ? * 

♦ •r> »4i -Pioroo 


, PS BE (4 1 ) , 

PSPSKI) 

• f*' *PSPHC4;i 

,r J X(4) )AM. 10510 


,YBF(4i) , 

YP$<41> 

YY(4!) » YX ( 4 1 > 

, YPH (41 )AMPL0520 


» XRF ( 41 ) 
,PHBE(41) 
,BMBE (41) 
»Q8F( 41 ) 


,XPS(4I 
, PHP' , 4i i 
,BM’'S (41) 
,. «»S(41) 


V Y (4 1 ) 

. >HV( 41) 
,BMY(41) 
, OY ( 4 1 ) 


, XX 1 41 ) 

, PHPH (41 ) 
,pMPH(.,l ) 
,0X(41) 


, XPH( 41 ) AMPL('*3'J 
, PHX (41 )AMPL0S40 
, BMX ( 41 )AMPL055C 
,OPHf 41 ) AMPL0560 
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EQU 1VALENCE 


EQUIvA.Ef.CE 


* -EL8F(«1» ,ELP$(4l) ,ELV(41) .HPWUll 

* ,PEBE(41) ,DFPS(4l) ,0EV(41) ,0EPM(41) 

* ,TB6(41) ,TPS(4l) * TY 141 ) ,TX(41) 

DIMENSION t'- AT (5 ,5 ) » VEC(5), DVEC(5), A(*>), DOH41, 

DIMENSION B(?,41 I,0UMY(1 ,1),S( I> 

1»BBM( 41 ) »C8M ( 41 ) »TOR( 41) 

EQUIVALENCE (DOM* 1 > ,BEBE fill 

EQUIVALENCE (ZB ( 1 > ,BEBE ( 1 ) ) , ( ZB <42 ) ,BEPS(1) ) . f ZB ( 83 ) ,BEX( ] 

1 (7B(124),BEP«(1 >>»(ZB(165),BEY( l)>, 

2 (ZX(1 ),XBE(1> ) ♦ ( ?X (42) , XP$(1)> (7X(83),XX(11 

3 (ZX(12<), XPH(1 >),(ZX(165), XV(i)). 

4 (7041)* QBEtl)), (70(421* Of E (1) ) , (70( 83 ) , QX( 1 

5 (70*1 24 )» QPH(1 >>,120(165), 0Y(1)>, 

6 (rL(l),ELBE l))»(7L(42),ELPS(l>»,(7t(B31 f ElX(I 

7 (ZL(124),ELP'-Ml)) ,(2L(165),ELY(1) I, 

8 i75 f 1 ) • PSBE ( 1 ) ) , ( ?S (4?) ,PSPS( I) 1 ,*7S(83),PSX( ] 

9 (ZS(12‘.),PSPMm>,(2S(165»,PSY) l)> 

EQUI v A.Ef.CE (2 V (1 >, YBE(l)) ,( 7Y (42) , YPS( 1) ) ,(7Y(83), YX(1 

1 ( Z v ( 1 24 ) , y?m(1)1,(7V(165), VY(1)1, 

2 (ZM(t),BMBE(l)),(7M(42»,BMPS(in,(2M(8?),p»X(] 

3 (ZM« 124>,BM<»N(1 )) , (ZM(165),BMV( 1)), 

4 (ZCO ),CEBE(1 M ,(70(42), DEPS(l)).(ZD(83),DExr 

5 (Z0U24',DEPM(1>),(ZD(165),PFY(1)), 

6 (7H(1. . PWBSd ) > ,(7M42) »PMPS( 1) ) ,(ZM(83) ,PHX( : 

7 (7H(124),PMPH(1 )),(ZM(165), PMY( 1 ) ) , 

•* (71 (1 ),TBF(1) ),(ZT(42), TPS(1)),(7T(83), TXC 

( ZT (124), TPM(l)), (77(165), TV( 1 ) ) 

DATA CVCPN/9 ,5402966/ , CV0T/5. 729578/ 

C****** ******************** *********** 

C MODF LOOP M*l POP COU'CTIVE MODE * 

C M*2 FOP CYCLIC MODES * 

C M*3 FOP SCISSORS MOOES • 

C ******** ******** ******** ****** ******* 

T =0 .0 
NCOLMeO 
NCYCM*0 
DO 227 Mel, 3 

C BYPASS C01LECT1VE ANP CYCLIC MOOES FOR t -.RTJCUIATFD ROTOR. 
1F( M,NE .3. AND.FHOFF ,NE »C » GO TO 227 
IF(M.NE.7.AND.CH0FE .NE.O) GO TO 227 
MOOENOO 

IF( IBE4TST,: J,m) .LE. IBS(IST,IB,M) )G0 TO ??7 
c******** ********* *** ***r- ****** ******** ******* 

C SWEEP NATURAL FREQUENCES SCORED IN SOMNAT * 

C ***» ******************** ********************* 

NPT • 1+ IPF(IS’,I8,M) -IBS (1ST, IB, M) 

OC 223 NP*1,NPT 

nps * np ♦iescsT,ip,M) -i 

SM )»SOMNAT ( NPS ,M ) 

CALL COEF(M,M,.TPbE.,l,S , DUMY) 

PNA7 * CVCPM *SOPT;S(D) 

SOrtNAT(NPS,M) e FNAT 
OFT*. TRUE. 

I********************************* 

C S OMATY eVEC FOR A 

C n.e> peTa(TIP)*.*( 1)*?(TIP) • 

C VOIP) * A( 3) *7 (TIP ) • 

C PST(T’» *A( 2) *7 (TIP) * 


,ELX(41)AMPL0f7O 
,DEX(4I)AMPL05P0 
,TPH( 41 )AMPL0590 
,5,10) AMPL0600 


»rt;culatfd rotor. 
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t M* 1*4 lUU *’»» tt»»\ JV4; 

t.' thr »r\|t ♦* Ihni ,«n Itir fit 




C PHJITIP) *A(4)*2(TIP) * 

€•♦**•****♦ •*••*•»*•******* ******* 

DO 250 1*1, MM5 
Alt 1*0*00 
•>0 2*0 J*1,MM5 
OMATII,J)*CM/T»I,j) 

2*0 AID* OMAXlt A(I) ,0ABS(0MATt 
250 Atn*1.00/A(l> 

DO 255 1*1, MM5 
DO 255 J*I,MMS 

255 OMAT(l v j)*DMAT(I»j|*Am 
P2*0. 

00 280 1*1, MM5 
IRW*0 

DO 270 J*1,MH5 
IF(J.EQ.T) GO TO 270 
iew*lftw*l 
DO 260 K«1,MM4 
260 CHAT* IRW,K)*0MAT) j,K) 

270 CONTINUE 

CALL INVDET(Pl) 

IFI ABS(Plt.L T .P2> GO TO 280 
JRW*I 

P2*A8S(P15 
280 CONTINUE 
IRW*0 

DO 300 I*l,MM5 
Afi »*0.00 

IFII.EQ.JRW) GO TO 300 
IRW*IRW*1 

VEC(IRW)«-0MAT(T,MM5) 

0V€C< IRW1*VEC(IRW) 

DO 290 J*1,MM4 
CMATI TRW, J ) *DMAT< I,J) 

290 DMATURW,J)=DMAT1I,J) 

300 CONTINUE 

DF7*. FALSE* 

CALI TNVDFTtOUNVV) 

KLM*0 

50 CONV*,TRUt . 

DO 70 1*1, MM A 

0A«0*D0 

00 60 J*1 ,MM4 

60 OA*OA*CMAT{ I,J)*OVEC( J> 

!F( DA .£0*0*00) GO Tr 70 
A (I )*AtT)«OA 

CONV*CONV. ANC.DA/ At I ) .LE . .000001 
70 CONTINUE 

IFICONV) GO TO 100 
KLM*KLM*1 

IFOCIM.GT.25) GO TO 70 
DO 80 1*1, MM* 

OVECf TUVECf II 
DO 80 J*1,MM4 

80 OVECm DVECU)-OMAT( 1,J)*A(J) 

GO TO 50 

90 WRITE 16,907) PNAT , DBOMfIBI 

07 FORMAT <2*M CONVERGENCE FAILURE AT ,F9.2,18M CRM, ROTOR RRM * 


AMFL1150 
AMPH160 
AMPLII70 
AMPLl 180 
AMP LI 190 
AMPLl 200 
ANPL1210 
AMPLl 220 
AMPL1230 
AMPLl 2*0 
AMPH250 
AMPLl 260 
AMP 11 ?70 
AMP LI 280 
AMPH790 
AMPL1300 
AMPL1310 
AMPL1320 
AMPLl ?30 
AMPL1340 
AMPL135C 
AMPLl 360 
AMPLl 370 
AMPL13P0 
AMPLl 390 
AMPL1400 
AMP 11410 
AMP LI *20 
APPL1430 
AMPL1440 
AMPL14S0 
AMPL1460 
AMPL1470 
AMPL1480 
AMPL1490 
AMPLl 500 
AMPLl 510 
AMPL1520 
AMPLl 530 
AMPLl 540 
AMP1155C 
AMPLl 560 
AMP11570 
*>F 11580 
AMPLl r 90 
AMPLl* 00 
AM*LltlO 
AMPH620 
AMPLl 630 
AMP L 1 64 0 
AMPL1650 
AMPU 660 
AMPL1670 
AMPL1680 
AMPLl 6«0 
amp LI 700 
AMPL17I0 
AMPLl 720 
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* F9.2) 

INODE (NP$,M) « 5 
IPLNINPE,M> * 4 
DEGfN?S»M> * 180. 

GO !j 223 

******** *•*****.« 

C CALCUtATF MODE SHAPES * 

£*********. ************** 

100 IF! CC02 ) A(4)*0.D0 
L*N1*1 

DO 115 1*1 1 Ml 
l»t-l 

Bfl *L )*YBE I I >*A 1 1 )♦ YPSI I 1*A I21^VYI I )*VX fl )*A|31 ♦VPHI I )*A( 41 
6f2»L>~XBE<I >«AUHXPS<I>*A<2l*XYm + XXm*A<3)4XPH<It*A<4l 
IFICC02) B I 3*L )*0.0 

IF! M .E0.3. ANO.FHPFF .NE .0.0 .AND. I »GE . INI— LFHl iBf I «L ) *0.0 
IFI M .E0.3. ANO.CHOFF .NE *0 .0 .AND .1 .GE. INI— LCH) IB ( 2 f L 1*0*0 
110 IF< .N0Y.CC0?)B<3,l)*<PMBEm*AI 1? *PHPSI1) *A I 2>*PHV| I )*PHX (I :*A. 

* *PMPH< I }*Af4))*tVDT 

115 IF I PMOFF.NE .0.0 .AND • I .GE. INI— L PH) )B (3 .L )*B| 3 *1* 1 ) 

ABSCL*0. 

DO 120 I *1 * MM3 
DO 170 J*1 ♦ N1 
ABSP* ABS f B f T ,J) > 

IFUBi^L.GT.ABSB) GO TO 120 
IPLNlNP$*“) = ! 

scale =en-j) 

ABSCt *A6$B 
T 20 CONTINUE 

SCALE*! ./SCALE 
ABSB* SCALE 
00 150 1*1*MM3 
IFII.E0.3) ABSB*10.*AB$B 
DO 150 J*I *N 1 
150 B(IfJ)*B<I, J)*AB$B 

IFI .NOT. POUT) GO TO 10 

1FIRC0LL! 1ST) . NE . COLPUN .OP. DBOMtlBl .NE . RPMPUN) GO TO 10 
C BYPASS MODE PLOTS IF THESE MODES ARE NOT TO BE PRINTED. 

CALL MDPIQ T f 8 *NP»MfNPT# FNAT .DB0MI1B) ) 

10 CONTINUE 
NNOOE * 2 
K*IPLN(NPS,M) 

AM2*BlK t 3).GT.B(K*2) 

DO 160 !*4,N1 

AMI *B IK? I ) .G T *8 IK tl-l I 

IFI .NOT. (AMI AND.AM21 .AND. I AMI .0R.AM2) t NNOOE *NNODE M 
IFINNOOC .vE.4) GO TO 180 
160 AM2*AM1 

IFINNODE .NE .2 .OR.K.EQ .3 .OR «K .EO. M) GO TO 180 

C*T 

IFIK.EQ-2) 1*3 
DO 170 I*l?Nl 

IFISCALE*fDUH(ia t L)PAm*DUMlI f 2 t L)PAt?)*0UMII ?3tL)*AI3>* 
0UMUf4 f L)*A<4)^DUM|I t 5 f L) I.GT.1.E4) GO TO 180 
IV CONTINUE 

NNOOE *1 
180 1 BSCL*0. 

CO 3i0 I *1 ?N 1 


AMPL1730 
AMPL1740 
AMPU750 
AMPL1760 
AMPL1770 
AMPL1780 
AMPL1790 
AMPL1800 
AHPL1810 
AMPL1820 
AMPL1830 
AMPL1840 
AMP LI 850 
AMPL1860 
AMPL1870 
AMPL1880 
AMPLl P Q 0 
•.MPL1900 
AMPL1°10 
AMP L 1920 
AMPL1°30 
AMPL1°60 
AMPL1950 
AMPL1960 
AMPLl 970 
AMP LI 980 
AMPL1°90 
AMPL200C 
AMPL2010 
AMPL2020 
AMPL2030 
AMPL2040 
AMPL2050 
AMPL2O60 
AMPL2070 
AMPL20PC 
AMP L? 090 
AMPL2100 
AMPL2110 
AMPL2120 
AMPL2130 
AMPL2160 
AMPL2150 
AMPL2160 
APPL2170 
AMPL2180 
AMPL2190 
A"P12?00 
A**P l? 21 0 
AMP12220 
AMPL2230 
AMPL2240 
AMP L? 250 
AMPL2260 
AMPL2270 
AMPL2280 
AMPL229C 
AMPL2300 
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IKi* til »li\» klMit i* «!f (LIj lilt tins JkX** 1 is 
suti^il tu thi* i»*stnt I hmi mi th«* til Ir* jkkft* 


F8M*Et*CTU>+DE*ST( J) 

FtD*0E*CT(J>-£L*STO> 

WRITE (6*911 1 J» ZtJ), B(?*J>* OBM* 0C0. FBM t FCO 

IF( .NOT.CC02) WRITE (6*912 ) P(3*J>* T 

tct-1 

IF(J.NE.Nl) GO TO 230 
220 CONTINUE 

IF (K AT. 3) WRITE (6*914) 

IF (K.EO. 3) WRITE (6*915> 

9|4 FORMAT (1WO*49X*34MNOTE ll> PER INCH MAX DEFL F CTION } 

915 FORMAT (1HO*51X*?1WNOTE (1) PER 10 DEG MAX TWIST ) 

GINT*0.0 
00 410 N*lt2l 

GINT*GINT*>$M!N>*<B( 1, N) **2*B( 2tN>**2>*EYR (N)*(B 1 3 *N >/57 # 3 )**2 
Iff *NOT *L0 V N5) CC TO 410 
IF(ABS(B f ~*N)> .LE. 0.01) 8(3*N)*0.0 
410 CONTINUE 
412 CONTINUE 

WRITE (6*407) GINT 

40? c OR MAT (//29H THE GENERALIZED INERTIA X* *F10.5* 

1 16M IN— L3F -SEC**2 ) 

225 CONTINUE 

CALL CARDS CM t MODENO*NPS# ICOL*IRPM, P I 
223 CONTINUE 
227 CONTINUE 
RETURN 

230 BMsfBMBE<n*A< 1)*BMPS(L >* At 2) *BMV IL >*BMX ! L>* A(3l*BMPHIl l* At 4) I 

• ♦SCALE 

0«C0BEm*Am*0*>S( L>*At2K0YtU*0XtL)*At 3)^0PW<l!*A<4> )*SCALE 
EL*tELBEU)*Am*ELPS!L)*At 2imv CL HELXt L)*At 3UELPHU >*At 4M 

* 9SCALE 

DE* tDE8E(l)*A(l )*DEPSfL )*Af 2)^DEY (L )*OEXtL )♦ At3 )*DEPH tL l*At 4) ) 

♦ *SCALE 

IF( .NOT.CCO?) TMTBEI L) FA t 1 )*TPS( L I *A t 2 1 ♦TY (L >*TX (l)*At 31 *TPH tL ) 

* *At4))*SCALE 

I F ( PHOF F .NE *0 • • AND* J *LE * LPH) T * 0 *0 
GO TO MBR1* ( 200* 220) 

END 


AMPL2890 
AMPL2900 
AMPL2®10 
AMPL2920 
AMPt2930 
AWPL2940 
AMPL2950 
AMPL2960 
AMPLZ^TO 
AMPL29P0 
AMPL2990 
AMPL3000 
AMPL3010 
AMPL3020 
AMPL3030 
ANPL3040 
AMPL3050 
AMP L3060 
AMPL3070 
AMP 13060 
AMRL3090 
AMPL3100 
AMP13U0 
AMPL3120 
AMPL3130 
AMPL3140 
AMPL3150 
AMPL3160 
AMPL3170 
AMPL3180 
AMPL3190 
AMPt 3200 
AMPt3?10 
AMP L3 220 
AMPL3230 
AMP13240 
AMPt3?*0 
AMP13260 
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0S/360 FORTRAN H 

OPTIONS - NAME * MfiTN.0P*-O2,lINECNT*60,SlZE*0000K, 

SOURCE »t BCD IC , NCI 1ST , NODECK »LOAD»MAP»TiOED IT , IO,XREF 
SUBROUTINE CARDStM,MPDENO»NPS* ICOL.IRPM, B » 

C ******** **************** A******************* ************** 

C THIS SUBROUTINE PUNCHES OUT NODE SHAPES • 

C******** **************** **************** ****** ************ 

COMMON /COMA/ JMUB, Nl, LOT, POUT ,ITLE 1191 ,NAME I 2 ) »ND I? > .NPG 

* , CO ATE f 2) ,JHU81,DBOM( 10 ) ,RCOLL t 10), ? (41 )• INPUN 

COMMON /COMO/ CMAT«5,5),Sn*«AT<200,3),IPLNI200,3),INODE»200,3), 

1 HM3,MM4,MM5,CTI41), ST(41 ) , IB ,IST, 

2 IBS( 10,10,3) , I BE ( 10, 10,3 ), I STS ( 10,3) ,ISTEtlO,3) 

COMMON /H1NGFS/ LCH ,L CHP1 ,L FH.LEHPl ,CHOEF,EHOFF,FCH,FFH, IPUNCT 

* ,R PM PUN, COL PUN 
*.LPW,LPHP1,PH0FF,FPH 
*,BOMN,TWSr .L0YN5 
*,LCTS 

*,BOMI ,TWS1,DF'_B0M,0ELTWS 

COMMON /COMT/ F VX (41), FYB t 41 ), EYCI41), YB 141) , YCU1), XIMII41), 

* X* T 1 41 ) , E YP (41 ) , EMRB(41), EMRCI41), EMPRI41), EMPS0I41 ), CC02 
*,OVPLT,OVLIN,-VLIN 

*, 8L AOE S ,HLfBTV D 

LOGICAL OFT 

LOGICAL CONV , AMI, AM?, POUT, CC02 
*,RBTEST 
*,LOTS 
*,LDYN5 

6,0VPLT,0VL!N,SVlTN 

C REAL *8 SPECIFICATION FOR IBM; COMMENTED FOR COC VERSION. 

REAL*8 CMAT 
DIMENSION 8(3,415 
2 ,PUN!2?,3> 

3»D( 6,3,6) 

DIMENSION YENtlC,3,T,3,A,3) ,YENFROI 3 ,3,4,3 ) 

C«*** ************************ ********* 

C MODE LOOP M*1 FOR COLLECTIVE MODE * 

C M*2 FOR CYCLIC MODES * 

C M*3 FOR SCISSORS MODES * 

C******** ********************** ******* 

MOO EN0*M00EN0*1 

IFfRCOI LUST).NE.COLPUN.OR.DBOM(IB) .NF.RPMPUN) GO TO 702 


DO 25 KK*1 ,21 
K*KK 

PUN(KK,1)*B( l.K) 

PUN <KK,?)«-B(?,K) 

PUN(KK,3)*B(3.K»* 12. 

25 CONTINUE 

DO 200 J*2 ,21 
RNcW»(j-15*2(?l>/20. 

FRACTeO.O 
DO 205 1*2,21 

IFf RNEM.GT.2 (T-l) .AND.RNEW.LE.2«I )) 
1FRACT=(RNEW-Z(I-1 )>/( 2 II) *2 (1-1 )) 

IF( FRACT.EO.O.O) GO TC 205 

PUN ( J , 1 ) *B ( 1 ,1-1)*FRACT*(BC1.I)-B(1,I-1 )) 

PUN (J»2)*B(?,I-1)*FRACT*(B( 2, I >-B ( 2 , 1 -1 ) > 
PUN(J,3)*(B(3,I-1 )*FR ACT* (6(3, I )-B( 3,1-1) ))*!?. 
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IF! FR ACT ,NF .0.0) GO TO 20« 

205 CONTINUE 
209 CONTINUE 
200 CONTINUE 

PUN(22,1)*S0MNAT'NPS.M)/DB0M!IG) 

IFfN.FO. 1>PUNI??,2>* I. 

IF!M.EQ. 2>PUN!?2,?>*-1. 

I r IM.E0. 3)PLN122,2>*0. 

PUNI22»3“=.02 

ANGLLH*0.0 

IFIZUl ).FO.O.O) GO TO 208 
IHNG«Z!41) 

ANGLLH*12.*57.3*!B! 2,IHNG42)-6!2,lMN64l ))/! 2UHNG-»2)-2IIHNG4l» I 
1 - IB 12 ,IHNG)— B 12 »IHNG-1 I )/ (2( IHNGI—2 1IHNG— 1 1 1 

WRITE!6,73)ANGLLH 

73 FORMAT! 73H TmE LEAD-LAG ANGLE FOR THIS HOOF .SCALED TO 1 FOOT MAX 
1 DISPLACEMENT,^, F10.S.10H DEGREFS* 

208 CONTINUE 

IF! INPUN .NE.l) GO TO 7 
IFIMOOENO.GT.fe) GO TO 7 
RBTEST*. FALSE. 

IF! ABS!PUNI2?,1)-1.).LE..05) RBTEST*.TRUE . 

00 450 KN*1 . 21 

IFIM.E0.2.AND.RBTEST) PUN !KN,3) *0.0 
450 CONTINUE 

IFIM.EQ.2.AND.RBTEST) PUN!??,3)*0.0 
IFIRC0LL(IST).F0.1.) GO TO 7 
IDBOM=OBOM( IB) 

ISW*Pl*M22,2> 

DO 600 KKIC*l»?l .2 

IFIRCOLLIIST) .NE.CPLPUN.OR.OBOMIIR) .NF.RPMPUN) GO TO 600 
KRRP1 *KKK*1 

WRITEI7.2T) IPUNIKKK.I) ,1*1 ,3), PUNIKKKP1 ,1), 1*1, 3). NAME, 
1MOOENO,ISK,RCOLL( 1ST) ,IDBOM 
600 CONTINUE 

WRITE !7, 30OANGLLW, NAME ,MOOENO,ISW,RCOLL! 1ST), IOBOM 
300 FORMAT! F10.6,50X, A4 , A 2, I 2,1 3, F4 .0,1 5) 

TO? CONTINUE 

IF! HODENO.GT .6) RETURN 
SORT CYCLIC DETUNING OATA FOR CB1 

IF! IB .EO. 1 . AND. 1ST. FO. 1)011 ,M,MOOENO) *SOMNAT!NPS,N) 

IF! IB.E0.1.ANC.IS T .E0.3)0(? ,M ,MOOENO) *SOMNAT|NPS.M ) 

IF! IB .EO.3.AN0.ISY.E0.1>0!3,M,M00ENO)*SOMNAT!NPS,M) 
IF(1^.E0.3.AN0.IST.E0.3)0!4,M,MP0EN0)*S0MNAY!NPS.m) 


27 F0RMAT(6F1C.6,A4,A2,I?,13,F4.0,I5) 

7 CONTINUE 

FOLLOWING 36 CAROS ARE COMMFNTED TO SAVF CORE SPACE. 

IF MOOES TO PIT BMC PROGRAM DYN5 ARE DESIRED, REMOVE C»S HERF 
AND ON DIMENSION S T ATEMENT FOR YEN ANO YENFRO. 

IF! LDYN5 ) GO TO 460 
IF IMOOENO .GT. 4) GO TO 469 
00 437 IS=G*1,10 
DO 437 ICPMP=1,3 

IF( IC0MP.P0.3)B!3,ISEG»24l)*B!3,ISEG*2*l)/57.3 
VFN ! I SEG , I ST , IB , 1 COMP ,MOOENO,M ) *B ! 1C0MP , I SE G*24 1 ) 

437 CONTINUE 

mENFRO!IST,IB,MOOENO,M)*SOMNATINPS,M) *6.28/60. 
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tK* .•» ili\i .'I *Utj n tint js**r i\ 

I t * 1 till' t l'\ll u (ltM| >>ll |ll|* ttfcft* 


469 IFIM .NE. 3) 60 TC 460 

IFIRCOLLdSTKNE.COLPUN) GO TO 460 
IF( OBPM ( IB ) »NE .RPMPUN ) GO TO 460 
IF (NFS ,NE. IBE( IST,IB,M))G0 TO 460 
OO 459 I TYPE *1,3 
00 459 !H0DEN*i,4 
WRITE (7,465>NAME 
WRITE (6,465) NAME 

465 FORMAT! 30X,A4, A?) 

WRITE ( 7 , 464 ) ( (YFNFRO! ICOL ,IRPM, IMOGEN,! TYPE » ,ICOl *1 ,3 )« !RPM*1 ,3) 
WRITE (6,467) ( ( YENFRO (ICOL, I RRM, IMOGEN, I TYPE ) »IC0L=1 ,3)»IRPM*1»31 

467 FORMAT (• • 1?F6.0) 

IVEN-’ITYPE-l 
WRITE (7,466) IYEN 
WRITE (6,466) IYEN 

466 FORMAT U5H 111 ,15 ) 

00 459 IC0MP«=1,3 

WRITE (7,463) ( ? ( YEN( TSEG ,ICOl , IRPM , I COMP , IMOGEN, TTVPE ) 

1 »ISE6*1»10) , ICOL *-1 ,3 ) , TRPMr 1 ,3 ) 

WRITE (6,468) ( ( (VEN( ISEG . ICOL , IRPM , 1 COMP ,1 MODEN, ITV?E ) , 

• I SEG* 1 , 10 ) , ICOL =1 ,3>.IRPMel,3) 

468 FORMAT (• • 13F6.3) 

459 CONTINUE 

460 CONTINUE 

464 FORMAT) 13^6.0) 

463 FORMAT! 13F6.3) 

IF( INPUN, NE.l) RETURN 
C PUNCH CYCLIC DETUNING CARDS FOR C81 
883 FOR MAT ( 71 1 0 ) 

ID80M*OBOM(2) 

ICOLL =* COLL (2) 

DO B IMel,3 
DO B IJ*1»6 

D(5,IM,IJ)*( RCOLL (3 )-RCOLL( 1 ) >*»5 
D(6,IM,IJ)*(D90M(3»-DE0M) 1) )P.5 
IF(IM.EO.l) ISW=1 
IFtTM.Eo,?' ’SW*— 1 
IF( IM .EC.' SW -0 

IF( IB.EO. ' NP. IS T . EQ ,3 • AN0.MPDENC.E0.6 . ANO,M,EO,T ) 

1WR I TE (7,28) ( 'Oil ,IM,IJ),I«1 ,6),RC0Ll )2),DB0M(2),NAME,IJ,ISW) 

28 FOR MAT(6F10.0,2F5.0,A4,A2,2I2) 

8 CONTINUE 
RETURN 
END 
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OS/ 360 FORTRAN » 

OPTIONS - NAPE* PAIN ,0PT*02,LlNECNT=60,SlZF*000OK, 

SOURCE *FPC0 IC *NOL 1ST *NOOECK .LOAD, NAP ,NOEDIT , I D, XRFE 
SUBROUTINE COEFIKIN, KAX, DET, IMAX, SQSOM, PP> 

C ******** ************ a***************** ****** ********** ••***••***• 
THIS SUBROUTINE CALCULATES THE DEFLECTION OF EACH STATION 
AS A FUNCTION OF DEFLECTIONS AT THE ROTOR BLADE TIP FOR 


NODES I A*K I N TO KAX 


I A* 1 FOR COLLECTIVE MODES 
I A* 2 FOR CYCLIC MODES 
IA=3 FOR SCISSORS MODES 


IMAX*NO OF FREQUENCIES TO BE CALCULATED * 

SQSOMO TO IMAX) C0NTA1NFS SOUARES OF FREQUENCIES * 

• 

DET*.TRUE. USED to find MODE SHAPE FOR KNOWN NATURAL FREQ. * 

• 

OET*.FALSE. used TO find NATURAL FREQ. -the determinates of * 
the POINDRY CONDITION MATR1CIES ARE CALCULATED * 
AND STORED IN PPfl TO IMAX, KIN TO KAX). • 

**************** ****** ********** ****** ***•**•••*****•**»•**• ****** 

COMMON /COMA/ JHUB, Nl, LOT, POUT, ITLEU9), NAME <2 ),ND»2 ) ,Nt 

*, CD ATE I ?: » JM'JBl ,DBOM( 10 ) ,RCOL L t 10 ) , Z f 41 ) , INPUN 
COMMON /COMP/ CK, I RCOl » XRCOL < 10) » IBOM, RB0MI10) , 

* SMZ< 41 ) * ZBARI40), EYEBU20), 

* EYEC ( 120 ) , SP f 40 ) , SC(40>. VMBJ40), VMC(40), VFP<<.0), VFCIAO), 


COFFOOIO 
COE F0020 
COEF0030 
COE F0040 
COFF0050 
COE F0060 
COE FOOT 0 
COEFOOPO 
COE F0090 
C0EF0100 
COE FOIl 0 
COEF0120 
COE FO 130 
C0EF0140 
C0EF0150 
C0EF0160 
C0EF0170 
C0EF01R0 
COE FO 190 
COE FO 200 
COE F0210 
COEFO 22(> 
C0EF0230 


* 0FB«40), 0FC(40), TH (41), THE ( 40 ) , WT!<,0), Sm< 4? ) , I SOM ,X0SOM ( 200 )COE F0240 

*,AZBAR,RPMA,RPME.,RPMC,COLlA,COLLB,COLLC,CHORD COEF0250 

*,RB (41) .PC I 41 ) C0EF0260 

COMMON /COMO/ CMAT(5,5),S0MNATC200,3),IPLN)200,3) , INODE ( 200 ,3 ) , C0EF0270 

1 HM3,MMa,mm5,C t (41), ST ( 41 ) , IP , I ST , COEFO2P0 

2 IBS( 10,10,3) » I BE (10,10,3),ISTS( 10,3) ,ISTFI10,3) C0EF0290 

COMMON/H/ VLX(40), V0XJ40), VLY (40 ) , VMXI40), V0XI40), VMYI40), COEF0300 

* D®LX(40) ,DP0X(40),DPLY(^0) *DPMX( 40 ) ,DPOX (40) ,0PM VI 40) , C0EF0310 

* DFLX»40) ,DFDX(40),DFLY(40),VFLX»40),VFDXC40),VFLV)40), C0EF0320 

* F|41>, POMS, 0TX)41), 0TYI41), SX»41), SYI41), EMRX(Al), COEF0330 

* EMRY (41), EM88W( 41 ), EMBB0I41), FMBPP(41), EMPPWI41), FMPP0I41), C0EF0340 


* TH HO ,41), CTXKI), FTYI41), WFLC41), WFD»4l), FMBPW(41) C0EF03S0 

COMMON /COMM/ SPRIP .FLPSPR COEF0360 

*,VSOF ,VMAS ,HSCF ,HMAS ,RSOF C0EF0370 

l.SOFI COEF03PO 

*, TORSO COEFOT-O 

1,ANC,LE,STR C0EF043O 

1,IL0C ,TANAL' COEF0410 

COMMON /HINGES/ LCH,LCMP1,LFH,LFHP1 ,CH0FF,FH0FF,FCH,FFH,IPUNCT COEF0420 

* ,R PM PUN ,CCLPUN C0EF0430 

*,LPH,LPMP1 , PHO e F , FPH C0FF0440 

*,BOmm,TwSm ,L0YN5 C0EF0<50 

*,L0TS COEF0460 

*,BOMI ,TvSI,DELPOM,DELTWS COEF0470 

COMMON /COMT/ E YX (41), EVp(41), EYC(41), VB(41), YCU1), XIMII41 ) ,COE F04P0 

* XI T t 41 ) , E YR (41 1 , EMR0U1), FMRC(41), EMRR 14' ) , EMRS0I41), CCO? COEF0490 

* ,0V PL T ,OVL I N ,SVL IN COFFO'OC 

PLAOES,HUBTYP C0EF051O 

COMMON /CO-E/ 78(?OS>, 2X(205>, 70t20*>>. 7L(204>, 7S(20*), 7V I 205 )COF F0S20 
*»ZM (205 ) , Z0I205), ZMI205), ZT(205) C0EF0530 


COF F0 K 0C 
C0EF051O 

7L ( 20*> ) , 7 S ( 20 * ) , 7V I 205 )COF F0520 

C0EF0530 


LOGICAL "ET, fCO? CrEFD*40 

REAL *P SPECIFICATION FOR IBM; COMMENTED FOR COC VERSION. C0EF0*50 

RFAL*R 2E,ZS,ZX,7Y,2M,2m,20,2L,2D,2T C0EF0S60 
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*, CHAT, OB, OV, OS, OX, OH 

DIMENSION SOSCMU), PPI20G,3) 

C ******** **************** ************ 

C CALCULATE DEFLECTION COEFFICIENTS * 

C******** ******** ************** ****** 

DO 220 I I»1 , IMAX 
S0MS*S0S0M(IT) 

SOM SY*SOMS* BOMS 

ZMC l ) *SOMS*EMPBMl N1 )*B0MS*E MBB0IN1 ) 

ZOI l)»SOMS*EMBPw(Nl)*BOMS*EMBPO(Nl) 

ZMI 42 )*ZQI 1 ) 

ZOI 42 >«SOMS*EMPPW INI ) *BCMS*EMPPOtNl ) 

ZOI 83)«S0MSV*SM(N1) 

ZTI 83)*S0MSV*EMRV»N1> 

ZLI 124)*S0MS*FMRXIN1) 

ZDI 124)=S0MSY*EMRV|N1 ) 

ZTI 124) *SOMS*EYR INI >*B0MS*THH0IN1 ) 

ZLI 165) =SMI N 1)*S0MS 
ZTI 165)*SnMS*EMRXINl) 

DO 135 J*1 , 5 

IF( CC02.AN0. J.E0.4) GO TO 135 
M*N 1 

Ll«J**»l-40 
DO 130 I*2,N1 
12*11 
ll*Ll*l 
M*M— 1 

DB v VF LY I M)*ZB ( 12) ♦VFL XI M) *ZSI 12)— VMY I M) *ZMI 
*— VI Yt M) *Zl I L ?)-VLX I M) *ZOI 12 ) 
DS*VFlXIM)*ZfML2)*VF0XIM)*ZSIl2)-VMX|M)*2MI 
•-VL XI M)*ZL I l?)-VOX I M) »ZOI L2 ) 
0Y«DFLYIM)*ZB'L2)*DFLXt M)*ZSIL2)-0PMVtM )*ZM 
*— DP LYtM)*ZL IL2) — D PLXt M)*?DI L2 )-DTV| M) *ZTI 1 2 
DX«0FLXIM)*ZBIL2)*0 K 0XIM)*ZSIL2)-0PMXIH)*ZM 
*-D01XIM)*ZHL2)-0K>X(M)*ZDIL2)-DTX|N)*ZTIL2 
DH*Wf LI M)*ZBIL2)*WFDIM)*ZSIL2 l-DTYIM )*ZL I L2 
•ZTI L2) 

ZBI LI )*ZBIL2)*0B 
?Slll)*ZSIL?)*DS 
ZYIll )«ZYIL2)*0Y 
ZXI LI )*ZX( L2 )*DX 
ZMI M )*ZHI L2 )*DH 

ZLI LI )*ZltL2 > .SC'MSMEMRXIMWMtll )*SM|M)*ZY 
ZDI LI >*ZDU2 /*SOMSV*EMRYIM)*ZMIH >*SOMSY*SM 
1000 CONTINUE 

ZMI LI )«FIM)*DY*FTX(M)»0H*ZMIL2)*ZBARIM)*ZH 
•♦ ISOMS*EMBBW IM) ♦BUMS*EM8B0I M ))*?BIL1)+I SOMS 
•)*2Slll) 

Z0IL1 )*FIM)*DX*FTY( M)*DH*ZOIL2)*ZBARIM)*ZOI 
•♦ISOMS*EMBPWIM)«.ftOMS*cMBPOIM) )*ZB(L1 ) 

•♦ I SOMS*EMPPw IM) ♦BOM $*EM*PC! M))*ZSIL1) 

ZTI LI )*FTXIM)*0B*FTYtM)*0S*ZTIL2)*IS0MS*FYR 
•)*SOMSY*EMRV (M) *ZX( LI )*S0MS*EMRX|M)*ZYIL1 ) 


L2)-VMXIM)*Z0IL2) 

L2)-VOXIM)*ZOI12) 

IL2)-0RMXIM)*Z0IL2I 

) 

IL2)-OPOXIM)*ZOU2) 

) 

)-DTX(M)*ZDIL2)-WTtM)* 


111 )) 

IM)*ZXIL1) 


* FMB f 'W IN )*BOMS*E MBPO I M ] 


IM)*B0MS*THHP! M) )*ZHIL1 


130 CONTINUE 
135 CONTINUE 

C******** ********* *«**»*••• ****** 

** CALCULATE BOUNDRV DEFLECTIONS • 
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I to tin* f«*\tcntn»n on t; • title p**)* 


£************************ *** 4*4 ** 

K*N1-41 
150 KK*0 

00 220 I A*KIN*KAX 

KK*KK *1 

L*K 

00 200 J«1,MM5 
L*t*4l 

IF<CC021 GO TO 193 

C TORSION COMPONENT BOUNDARY CONDITION 

CMAT(5 t J)*7T(L-LPM)-PPM*| Zi a-LPH)-7Ta-LPHPl)> 
1 -CK* (ZMfL-LPH)-FPH»fZHU-LPH>-ZM(L-LPHPl )) 

GO TO 197 

193 IF'J.EQ.4) L*L*41 
197 CONTINUE 

GO TO t 400 t 500* 600) *1A 
400 CONTINUE 

C BOUNDARY CONDITIONS FOR COLLECTIVE MOOES 

CMAT(l f J) * 7UI) — ZY f L ) *f VSOE 4S0MS YVMAS ) 
CMAT(2 f J) * ZX(L> 

CMAT(3tj> * zeu.) 

CMATf 4, J ) « ZOft) -7 SU) * TORSO *• .Efc 
GO TO 200 

500 CONTINUE 

C BOUNDARY CONDITIONS FOR CYCLIC MOOES 
CMATfl.J) * ZY ( L ) 

CMATf 2» J ) * ZOCL) - Z Xf L ) *fHSOF 4S0MS *HMA$> 
CMATf 3* J ) * 7M(l) - ZBCL) *R$OF 
CMATf 4* J ) * ZStL) 

GO TO 200 

600 CONTINUE 

CMATf 1 * J ) * ZYfL) 

CMATf 2, J) * ZX(L) 

CMATf 3* o f « ZB f L) 

CHAT f 4* . ) * ZSf L) 

I F f CHO r F * E 0 ♦ 0. 0 ) GO TO 300 

CMATf r t J)«r?x fL-LCM)-FCH*( 7X f L -LCM )-ZX (L -L CHP 1 >) 
CMATf 4* J)*ZO(L-LCW)-FCHPf ZOfL-tCM>-ZOfL-LCMPl>) 
1 - SOFI*(ZSf l-LCH)-FCMtfZS(L-LCM)-ZS(L-LCHPn> 
ZSf LCH) *0.0 
300 CONTINUE 

IFf FHOFF .FO.O.O) GO TO 705 

CMATf 1, J)=?Y<L-LFH)-FFM*f?YfL-LFH)-ZVfl-LFMPl>) 
CMATf 3 f J)*ZMCL-LFM)-FFH#f ZM f L-CFW WM fL H. FHP 1) > 
1-RSOF •f?B<L-CFH)-FFM*f7BfLH.FH)-79fL-CFHPl )) 
305 CONTINUE 
200 CONTINUE 

IFfOFT) RE’JRN 
220 CALL INVDETf PPfll *KK) ) 

RETURN 

END 


COE FI 150 
COE FI 16^ 
C0FFH70 
COEF110O 
C0EF1X90 
C°F FI 200 
COE F 1 21 0 
COE r l 220 
COEF1230 
C0EF1240 
C0EF1250 
COE FI 260 
C0EF1270 
COE FI 280 
COE FI 290 
C0€ c 1300 
CPEF1310 
CPEF1320 
C0EF1330 
CPEF134C 
CPEF1350 
CPF FI 360 
CPEFn70 
CPEF1380 
CPE FI 390 
COE FI 400 
CPEF1410 
CPFF1420 
CPEF1430 
COE FI 440 
COE *1450 
CPF F 1460 
CPE F1470 
COEF14BO 
CPE FI 490 
COE FI 500 
C OE f 1 5 1 0 
C OF FI **? 0 
COE *1*30 
CPFF1540 
CPE Fi *50 
COE F 1 56 0 
CPF F 1 c 70 
CPE Fl 58 0 
C OF F l 590 
CPF F 1 600 
CPE F 1 61 0 
C Of F 1 *?o 
CPEF1630 
C^E FI 6^ * 
COE F 1 
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OS/ 360 FORTRAN M 


OPTIONS - NAMF* MAIN »0 p T *02, LlNECNT*feO »S IZE=0000K, 

SOURCE t EBCDIC ,N0L I ST .NODECK , LOAD, MAP, NOEOIT , IP.XREF 
SUBROUTINE INPT(TYCE, FIRST) INPTOOIO 

C***9 ****** ************** ******** **************** ****** ******** INPT0020 

C THIS SUBROUTINE RFAOS AND PRINTS OUT INPUT DATA • 1NPT0030 

£*************•*****•*****************•****«*****•*******•****• INPTOOIO 

COMMON /COMA/ JHUB, Nl, LOT, POUT, I TLE ( 19) ,NAMF ( 2 ) ,N0 (2 ) ,NPG INPT0050 

*,CDATE( 2) .JHl'Bl ,090M( 10), RC 01 1(10), 71*1), INPUN INPT0060 

COMMON /COMB/ CK, IRCOL, XRCOL(IO), I BOM, RBOM(lO) , INPT0070 

* SM7(41), ZP AR ( 40 ) , EYEBU70), INPTOOBO 

* EYEC(120),SB(40) , SC (40), VMB(40), VMC 140) , VFBI40), VFC (40) , INPT0090 

* OF B ( 40 ) , DFC (40) , Th(41), THE (40), WT(40>, SM(42 ) , ISOM.XQSOM (200)INPT0100 

*,AZBAR,RPMA,RRMB.RPMC ,COLLA,COLLB, COLIC, CHORO INPT0110 

*,RB(41),RC(41) INPT0120 

COMMON /BIRO/ DUMMV(&3),M(6), INPT0130 

1 WTPK21) ,Eie(?0),FlC(?0),GA(?0),GI(?0), INPT0140 

* THD(?l) INPT0150 

COMMON /COMC/ N,IER (7), OFFSET INPT0160 

COMMON /COMF/ SOMM, TWIST, OIA, SOMI, DELSOM INPT0170 

COMMON /COMM/ SPRIP ,FIPSPR INPT01P0 

*,VSOF ,VMAS ,M$OF , MM AS , RSPF INPT01O0 

1 ,SOFI 1NPT0200 

*, TORSO INPT0210 

1,ANGLE,STR INPT0220 

1 , IL OC ,T ANAL F INPT0230 

COMMON /HINGES/ LCM,LCMP1,LFH,LFHP1 ,CHOFF,FHOFF,FCm,FFH, IPUNCT 1NPT0240 

* ,RP*»PUN,COLPUN INPT02S0 

*»LPM,LPHP1,RH0FF,EPm INPT0260 

*»BOMM ,TWSM , LDVN5 INPT0270 

•.LOTS INPT02B0 

*,BOMI ,TWSI,0E*.BO4,DELTwS INPT0290 

COMMON /COMJ/ VSOFT ,HMASS,VSOFT,VMASS,RSOFT INPT0300 

l,SPRLG INPT031 0 

COMMON /COMT/ EVX(41), EV9C41), EVC (41 ) , VB(41), VC ( 4 1 ) , V)MI(4l) .1NPT0320 

* XIT(41), E YR (41), FMRB (41) , EMRC(41), EMRR 1 41 ) , EMRS0(4l), CC02 1NPT0330 

* ,0V PL T,PVL IN ,SVL IN INPT0340 

*, BLADES, HUBTYP INPT0350 

LOGICAL LTWS, 12 BAR, LEIB, LEIC, LGA, LGAMp, LG I , LGAMC, INPT0360 

♦FIRST, LOT, POUT, LEVP, LEVC, LGAM, CC02 INPT0370 

*,0V PL T,OVL IN ,SVL I N INPT03P0 

*, SOUTH INPT0390 

* ,L0VN5 1NPT040E 

* .LOTS INPT041C 

DIMENSION W(?l ) ,RH01 XX(21) ,RH01VV(2I ) INPT0420 

DATA lASK.IPLS/'**, •/ INPT0430 

I,IVM/*E*/ INPT0440 

*, SOUTH/. FALSE./ INPT045P 

DATA NDECK/ * DECK* /, NNAME /• NAME •/ , NPLOV ‘PLOT* /, NMODE/' MOOE •/ INPT04B0 

* , NPUNCH/* P'JNC •/ , NT OR/* TORS* /, NT Wl ST/» TWI S •/ INPT0470 

*,NDVN5/*DYN5V INPT0480 

*,NALLMD/»ALLMV INPT0490 

NAMELIST /INPUT/ SOMI, SOMM, DELSOM, BOMI , BOMM, DELBOM, R, JHUPINPT0500 

*, N, LTWS, TWSI. TWSM, DElTwS, TWIST, THD , INPT0510 

* WTPL,7BAR,EIB,E1C ,GA,EVEB,EVEC,G1 .SB, SC ,RB,RC INPT0B20 

•»CC 0? ,OVPL T , OVL IN ,S VL IN INPTDF30 

*, SPRIP .FLPSPR INPTO«- 

*, TORSO INPTOF 7 

*, VSOFT, VMASS ,HSOF T, HM ASS ,R SOFT INPTCS60 
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* *CMOPF,FHOF* 
l*$PRLG 9 CYCLE 9 AZBAR 9 BL AOES 

♦ 9 PHOFF 

2*RPMA tRPWB*PPHC .COLLA 9 C0LLB .COUC 

itCK 9 ANGLE 9 SIR 

**P$QR 9 PLAST f OP 

UMUBTVP 

1 flTLF 9 NAME 

**Z 

1 CONTINUE 

OVPLT*. FALSE . 

OVL IN*.FALSE. 

SVLIN * .FALSE* 

L2BAR=. false . 

LETB *. FALSE. 

LEIC *• FALSE • 

LGA *. FALSE. 

LOTS*. FALSE . 

LGAMB*.TRU r . 

LGA MC*. TRUE. 

LEYB *♦ F ALS C » 

LEYC *. FALSE. 

LGI *. FALSE. 

£****************** 

N * 20 
1 f41 ) *0.0 

C READ TITLE CARD * 

C******** ********** 

IFf LOT. AND. SOUTH) 00 TO *00 
SOUTH*. TRUE . 

10 READI5 f °01tEND*?20)M 
901 FORMAT( A4.6X t A4 f 6X. A4 f 6X v A4«6X 9 A4 9 6r, A4 > 
K * I YM 
LOT * .FALSE. 

POUT # .FALSE. 

INPUN * 0 

cco? * .true. 

I TM$ * .FALSE. 

LDYN5* TRUE. 

DO 96F 1*1 1 f 
I Ff M( I ) .EC .NOEC A > K» IASK 
IFf Hf I ) .EO.NNAHE ) y« IPLS 
IFf MC I ) .EO.NMODE ) POUT *.TRUE . 

IF(M( I).EO.NALLMD)LOTS *.TRUE. 

IFf MU) .EO.NPLOTUOT«.TPUF. 

IFfMfD.EC .N PUNCH ) I NPUN*1 
IFf Mf I).E0.NTOR)CC02«.FALSE . 

IFf Ml T) .EO.NTVIST )L TWS* .TRUE. 

IFf Mf I) .EO.NOYNM LOYNS*.FAL SE . 

965 CONTINUE 

IFfK.EO.IYM) GO TO 320 
NPG*0 

IFf K.EO. TASK ) GO TO 20 
IF IK .NE. iPLS) GO TO 10 

C******** ** ******* *9* ** •* ** **** ** 


C READ CHANGES TO PREVIOUS CASE * 

£ ************************** ****** 
TWIST * TwSAVc 






INPT0570 
INPT05B0 
INPT0E90 
INPT0600 
INPT06I0 
INPT0620 
INPT.I630 
INPT0640 
INPT0650 
INPT0660 
INPT0670 
1NPT06B0 
1NPT0690 
INPT0700 
INPT0710 
I NP TO 720 
INPT0730 
INPT0740 
INPT0750 
INPT0760 
INPT0770 
TNPT07P0 

inrtotro 

INI -OBOO 
INPT0P10 
INPT0F20 
TNPT0B3C 
INPT0B4 ? 
INPT0P50 
INPT0B60 
tNPT0870 
iNPTOPBO 
INPT0B90 
INPT0900 
INRT09I0 
INPT0920 
INPT0O30 
INPT0940 
INPT0OS0 
INPT0960 
INPT0970 
INPT09B0 
INPT0990 
INPT1000 
TNPT1010 
INPT102C 
INPT1030 
INPT1040 
TNPTI050 
INPT1060 
INOT1070 
INPT10B0 
INPT1090 
TNPT1 100 
INPTU10 
1NPT) 1?0 
INPT1130 
1NPTH40 
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REAP * 5, INPUT } 

*MOVLIN> PQMl*AM AX 1 ( 1. ,BOMI ) 

^0 TO 50? 

C SET UR FOR ^COUPLED MOOES FOR FAN PLOT 

500 BOM 1*1 

0€LB0M«RPMA 
SOUTH*. FALSE . 

SVL1N* .TRUE . 

TW$ 1*1 • 

TWS M* 1 . 

TWSAVF * TWIST 
TWT S T *0. 

OYPtT*.TRUf . 

OVL IN*. T RUF • 

SOM I* .2 

OF 1 . S0M*0.05*RPma 
GO TO 50 

C ******** *********** 

C READ IN NEW CA^F * 

20 READ! 5,90?) NAME ,( IT IE (I) ,1*1,10) 

902 F0RMAT(4X,A4,A2,20X,10A4> 

iF f INPUN. EQ* 1) WRITE <7 ,600 )N AME , I ITL E tl ) ,1*1,10) 

600 FOR MA T f 10X,*4,A?,9X,10A4) 

REA Df 5,904) CYCLE ,T PR $0 , VMA SS ,HMA$S ♦ VSOFT ,H$OFT , RSOFT 
SPRIP * HSPFT 
FLPSPR * RSPFT 
JHUB = IFIXfCYCLF) 

REAP (5, 2) AZBAR, RPMA ,RPM8 ,RPMC tCOLLA *LOLLB ,COL L C ,TWI ST , 

* BLAOES,CHPRO,PSOR ,OP,PLAST,HUBTVP 

♦ ,C K)FF ,FHOF F 
1, SPRIG 

3,PH0FF, ANGLE ,STR 
TWSAVF*TWIST 
2 FORMAT (14E5.0) 

00 400 1*11,10 

400 IT\E(I)*I#SK 
502 IF( AZBAR.NE.O.O) R*N*AZ3AR 
BOM I * RPMA 
BOMM * RPMC 

TFfRPMC .FO. 0.0) POMM * RP^B 
I F ( RPMp .FO. 0.0) 50MM s RPMA 
DFL90M * 100. 

IF( (BOMM-BOMI) .NF. 0.0) OELBPM * ! BOMM P0Ml)/2. 

TWS I = COLL A 
TW$M * COLIC 

IFfCOLLC .80 • 0.0) TWSM * CPLLB 
IF ( CPLlB .80. 0.0) TWSM * COllA 
DEL * ; * 10. 

TF ••VT'.'SI) .NF. 0.0 OELTWS * < TWSM-TWS I )/? . 

IF T P •EO.O.OPSOR*.! 

: ' ** ' OP *BPMI 

SOMM * -VAST *BOMM 

DEL COM - OP *BOMI 

IF( SOMM.LE.O.O) SOM m* 10.4 BOMM 

2F(DELS0M.LE .0.0) D Ft SO M* .? 5*B0MM 

IF(K.FO.IPLS) GO TO 50 

JMUB*M* XOfJHUBfO) 


1NPT1 150 
INPT1160 
INPTH70 
INPTH80 
INOTH90 
INPT1200 
INPT1210 
InPT1??0 
INPT123C 
TNPT1 ?40 
INPT1250 
INPT1260 
INPT1?70 
JNPT12P0 

inpti?qo 
INPT1300 
INP T 1 ?] 0 
INPT132C 
1NPT1330 
INPT1?40 
TNPTl 350 
INPT1360 
INPT) 370 
INPT13P0 
INPT1390 
INPT1400 
INPT1410 
T NPT) 420 
TNPTl 430 
INPT1440 
INPT1450 
INPT*460 
I N°T 1 470 
TNPT* 4F0 
INPT1490 
1NP T 1 500 
INV <1510 
TNPTl *?0 

INPTH30 
INPT1540 
INPTlS^O 
;NPT15so 
INPT1 «70 
INPT1 SPO 
TNPTl 590 
1 NP T 1 600 
TNPT161C 
TNPTl A?P 
INPT1630 
TNPTl 640 
1NPT1 650 
!N paw n 
TNPT1670 
« NPT 1 6P0 
l NPT1 690 
TNPTl 700 
NPT1T10 
!NPT17?C 
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Use or disclosure ol data on this page is 
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CALI DATE(ND) 

Ni*N*l 

IM AZBAR*GT .0.0) CO TO 145 
Zfl)*0.0 

REA*}! 5.904) (7fI).Is2*?l) 

R*0 *0 

00 147 1*1 *20 
ZBARm*Z<I*l)-Zm 
Rtft* ZB AMI) 

147 CONTINUE 
145 CONTINUE 

REA Of 5t 904 | OUMMV 
602 FOR NAT( 7F10 *4) 

READ (5.901) OUM 
904 F0RMATI7F1C.0) 

FtPlNT*0*0 

totmas*o*o 

Zfl )«0.0 

00 960 ISEG*lt?0 

WTPK I$EG)*OUPPY< ISBG*42) 

EI8 (I $EG)*DUMMY ( I SEG+1 !*10**6 
EIC (1SEG)*DUMMY< TSEG+22)*10**6 
GAf I$EG)*1. *10**6 
I*ISEG 

IF( AZBAR.NE.0.0)Z(I*1 >«Zm*AZBAR 
TOTMAS*TOTMAS«*TPl( I) 

FtPlNT*FtPINT4WTPtm4(R*«2«l-n/40.l^9? 
960 CONTI WE 

WTPL(21)*0UMMY<631 
T0T«AS*TCTPA$*R/20. ♦ WTPt(21l 
FLP1NT*FLPINT4WTPU 21 l*R*R 
FLPINT*FLPINT/<3? .2*1 44* > 


TlPWT*WTPt(21l 

C 

IF(lTWS) REA0(5.904) (THO (I) .1*1 .N1 ) 

IFC CCO?) GO TO 45 
READ (5.9041 OOP eCK 
READ! 5.^04 1 DUMMY 
00 962 I$EG«1. 20 

EYEB( lSEG>*OUMMYf ISEG) 

EYECf ISEG)*OUMMY< TSEG+21) 

GI (I SEG )*DUMMY f I SEG+4? )*lO**t 
962 CONTINUE 

READ! 5.904) ( SB ( I ) . I*l.N) 

REAO( 5.904 ) (SC(I). I*1.N) 

READ! 5.904) (RBC I)>T*1.N1 I 
READ! 5.904) ( RC( I) .1*1 .Nl) 

00 3B6 1*1.20 
DUMPY (1 >*EYEBf I ) 

DUMPY (1*21 ) *EYfC ( I I ♦WTPL II 9*RC If 1**2/386 *4 
386 CONTINUE 

DUMMY! 21 )*0.0 
DUMMY!42>*-1.0 

C » » » » » » »» »»♦♦»» »» »»»»♦♦» 


1NPT1730 
1NPT1740 
1WT1750 
INPT1760 
1NPT1770 
INPT1780 
INPT1790 
INPT1600 
INPT1P10 
INPT1R20 
TNPT1P30 
INPT1P40 
INPTIB50 
INPT1B70 
INPT18B0 
INPT1 P90 
INPT1900 
INPT1910 
INPT1920 
INPT1930 
INPT1940 
INPT1950 
INPT1°60 
INPT1970 
INPT1980 
INPT1990 
INPY2000 
INPT2010 
INPT2020 
INPT2030 
INPT2040 
INPT2050 
INPT2060 
INPT2070 
INPT2080 
INPT2090 
INPT2100 
INPT2110 
INPT2120 
INPT2130 
INPT2140 
INPT2150 
INPT2160 
INPT2170 
INPT2180 
INPT2190 
INPT2200 
INPT2210 
INPT2220 
INPT2230 
1NPT2240 
INPT2250 
INPT2260 
INPT2270 
INPT2280 
INPI? >90 
INPT??00 
INP T ?310 


£ ^ 5 ^ 
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C CALCULATE PARAMETERS FOR PRINT OUT * 

50 IFI.N0t.CC02) GO TO 55 
45 DO 40 I*1,N1 
SBf I>*0. 

SC) I)«0. 

RBf f)*0. 

RCf 1)*0. 

EYESf 1>~0. 

EYECfI>*0. 

WTfI>«0. 

40 EYX fl 1*0 . 

EYXfNl>*0. 

55 JMUB1*JHUB«-1 

IF f INPUN. NE.l) GO TO 57 
IFf AZBAR.EQ.0.0) MRITEt6,202) 

202 FORMAT) 1M . 5 8HTHE INPUT MASS AND INERTIA DATA FOR UNEQUAL SEGMENT 
1LENGTM//41HMAVE BEEN RECAST BEFORE BEING PUNCHED OUT I 
WTT *0 .0 
DO 205 1*1,20 
Mil >*0.0 
RMOIXXt I )*0.0 
RHOIVY»I)*0.0 
X I* 1 1 —1 > *R/2 0. 

XIPl*I*R/20. 

FRAC*0.0 
DO 206 J*l,20 
VJ*2f J> 

YJP1*Z) J*l) 

IFfXI.GT.VjPl.OR.XIPl.LE.VJ) GO TO 206 

IFf XI.LF.VJP1.AND.VJP1.LE.XIP1.AN0.VJ.LE.XI) FRAC*YJP1-XI 
IFf XI .LE.YJP1.AND .V JP 1. LE . X IP1. AND. XI. LE. YJ. AND .YJ.LE.XIP1) 
l FRAC*V JP1— VJ 

IFf VJ.LE.XIP1.AND.X1PI.LE.VJP1) FRAC*XIP1-VJ 
IFtYJ.LE.XI .ANO.XIP' . LE.YJP1) FRAC*XIP1-XI 
Wfl >*Wf mwTPLf J)*FPLC 
RHOIXXt I )*RHOIXX( I ) ♦EVEBf J)*FRAC 

RHO IYVf I )*RHOIVYf I) ♦EVEC t J) *FRAC«Hf 1 >*RC) 1 >PP2/3B6.4 
206 CONTINUE 

Mfl >*WfI)*20./R 
PHOIVVf I )*PHOIVVf I) *20./R 
RHOIXX) I)*RH0IXXfl)*20./P 
WTT*WTT*Wf I ) 

205 CONTINUE 

W f 2 1 ) *T IPWT 

WTT *WTT*R/?0 . ♦Mf?l> 

RMOI XX 1 21 ) =0.0 
RHOI YV t21)*0.0 
WRI TEf 7, 210 MW(1> ,1*1,21) 

WRITE l7,210>fRHOIXX«I),I*l,21) 

WRITE <7, 210MRH0IVVU) ,1*1,21) 

210 FORMAT) 7F10. 5) 

57 CONTINUE 

CALL START) TOTMAS »FLP INT,TIPwT,R, LTWS ,F !R ST ) 

RET UR) 

320 CALL PLTIME 

TALL PLOTfO. ,0. ,499 ) 

TVC5 * 1.0 

RETURN 

ENC 
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HEUCX>PT6R cxiMnuvv 


0S/360 FORTRAN H 


OPTIONS - NAME* MAIN,0PT*0?,LINECNT*60,SIZF*0000K, 

SOURCE, EBCDIC.NOLIST.NOOECK, LOAD, MAP, N0ED1T, ID, XREF 
SUBROUTINE INVOET(O) 

THIS SUBROUTINE INVERTS ANO FINOS THE DETERMINANT OF A SOU ARE MATR] 
STIFF INPUT MATRIX — INVERSE UPON RETURN 

N ORDER OF STIFF IN BV Nl 

0 DETERMINANT UPON RETURN 

OET .FALSE. — INVERT STIFF AND FIND P 

.TRUE. — FINO P ONLY ISTIFF IS DESTROYED) 

NS I SIZE OF ARRAY STIFF IN THE CALLING PROGRAM INSZ BV 

IGOOFED RETURNS O FOR NO ERROR CONDITION, 1 IF OVERFLOW OR 

DIVIDE CHECK OCCURS, AND 2 IF MATRIX IS SINGULAR 
CHECK OCCURED. INOT USED IF DETb.TRUE.) 

COMMON /COMD/ CHAT (5,5) .SOMNATI 200,3 ) ,1 PLNI 200,3) .INODE 1 200,3) , 

1 MM3,MM4,MM5,CT(41), STI4I), IB .1ST, 

2 IBSI 10,10,3) ,IBE <10, 10,3 ),IST$( 10,3), ISTE (10,3) 

COMMON /COMI/ OET.MSZ, IGCOFD • SOM, OVRG 
LOGICAL OFT 

INTEGERP2 NDEX 

REAL *8 SPECIFICATION FOR IBM; COMMENTED FOR CDC VERSION. 
REAL*8 BIGA, HOLD, STIFF, P.CMAT 
DIMENSION NDEX (50) 

DIMENSION STIF F ( 5, 5) 

EQUIVALENCE (MSZ, N) , ICMATll.l) , STIFF) |,1) ) 

DATA NSZ/5/ 

L*1 

I GO OF D =0 
Pel .00 

SFARCH FOR LARGEST ELEMENT 

DO 80 K*1,N 

IF(DET) L«K 

IFIK.FQ.N) GO TO 45 

RIGAeO.DO 

DO 20 J«K,N 

DO 20 I*K,N 

HOLD e DABSISTIFFII ,J)) 

I FI BIGA.GE .HOLD ) GO TO 20 
BIGA * HOLD 
IRON el 
JCOL ej 
D CONTINUE 

INTERCHANGE ROWS 

NDEX(K) * JCOL • NS2 - NSZ • IROW 
IF ( IROW.LE .K ) GO TO 35 
DO 30 IeL,N 
HOLD e— STIFF (K, I ) 

STI FF IK »I ) e STIFF) IRON, I) 

0 STIFF (IROW, I) e HOLD 
INTERCHANGE COLUMNS 
5 IF! JCOL.LE.K) GO TO 45 
DO 40 J»L »N 
HOLD e -STIFFIJ.K) 

STI FF I J »K ) e STIFF! J, JCOL) 

0 STIFF (J, JCOL) e HOLD 

DIVIDE COLUMN BY MINUS PIVOT 
5 BIGA e-sTIFF(K.K) 

IFIBIGA.EO.O .DO ) GO TO 160 
DO 55 1 C*L »N 
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Pm- 01 >t»\, lOM'tl* t»l il.tt.l lh|\ ivhio tv 
Mi[«NS( t«* fft«* <|*vll «* !h‘»l **** ilU' |VN)«' 


55 IFIIC.Ne.KI STIFF*IC,K» * STIFF! IC#K) /RIGA 

c Reouce matrix 

PP 65 I*L *N 
IFII.FO.K ) GO TO 65 
DO 60 J*1,N 

60 IFIJ.NE.K) STIFF 1 I , J) * STlFFn,K|FSTIFF(K,J|4STIFFII.J» 
65 CONTINUE 

IFIOET) GO TO 7T 
C DIVIDE ROW BV PIVOT 

00 75 JR *1 » N 

75 IFCJR.NE.K > STIFFIK,JR) * STIFF! K , JR >/ST IFF < K ,K> 

C REPLACE PIVOT BY RECIPROCAL 

77 P«P*STI FFIK » K ) 

Q*P 

ifioet.ano.k.eq.ni return 

80 STIFF IK »K) * l.DO/STIFFCK.KI 
C FINAL ROW ANO COLUMN INTERCHANGE 
K*N 

100 K<=K-1 

IFIK.LE.O) GO TO 150 
J * I NDEXIK) - 1) / NSZ 
IROW * NDEXtKl - J * NSZ 
IFIIROW.LE.K) GO TO 120 
DO 130 1*1, N 
HOLD * STIFF ( I, K) 

STIFF(I,K> * -STIFF (I ,IROW) 

130 STI FF (I , IROW J * HOLD 
120 JCOl *J*1 

IFI JCOL .LE.K ) GO TO 100 
DO 110 J*1 ,N 
HOLD* STIFFCK.J) 

STI FF IK , J) « -STIFFIJCOL, J) 

110 STIFFIJCOL, J1 * HOLD 
GO TO 100 
150 CONTINUE 
RETURN 
160 CONTINUE 
0 * 0.0 
return 

END 


INVP0570 
INV DO 580 
TNVP0590 
I NV 00 600 
INVP0610 
INVP0620 
INV00630 
INV00640 
INV 00650 
I NV 00 660 
INV00670 
INV006B0 
INVD0690 
INV 00700 
INV00710 
INV00720 
INVD0730 
INV007A0 
INV 00 750 
INV00760 
INV00770 
INVD0780 
INV DO 790 
INV 00800 
INV00810 
INVD0820 
INV 00830 
INV008A0 
INVr0850 
INV DO 860 
INV00870 
INV00880 
INVD0890 
INVOOOOO 
INV 00910 
INV00°20 
INV 00930 
INV009AP 
INV00950 
INV 00 960 
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Uni* «it th\i liiMiir (4 <UIj uti this ****** is 
st'lip'it lit Hit* M*sln<tuMi nti Hit* lillr |wm* 


OPTIONS - NAME * MAIN ,0PT*02* LINECNT*60 .S 1ZF*0000K f 

SOURCE « EBCD IC *NOL 1ST ♦ NODECK tLOAD*MAP f NOFDIT * ID*XR EF 
SUBROUTINE I TER I IM* SOMA *S0M5* PA f P5) 

£*#********************** ******************************** 


c this subroutine iterates to THE NATURAL FREQUENCY I 
c 

C IM*MQDE USED AS ARGUMENT IN CALL TO COSE 

C S0M4 and SOM5* SOUARES OF FREQUENCIES THAT 

C BRACKET NATURAL FREQUENCY 

C PA AND P5* DETERMINANTS OF BOUNDARY CONDITION 

' MATRIC1ES ASSOCIATED WITH SOMA AND 

C SOM5 

(>♦**♦*** ********* *********** A**# *********************** 


* 


COMMON /COM I / D ET *M SZ » IGOOFD , $OM t OVRG 
LOGICAL XF * OVRG 
DIMENSION Sll)* Pll*l> 

OVRG* .TRUE. 

XF*P4.LT.O. 

P1*ABSIP4) 

P2*-ABSIP5> 

SOM 1* SOMA 
SOM 2* SO MS 

Sin*.5*tS0Ml4$0M2> 

ICOUNT*0 

350 CALL COEFIIM.IM,. FALSE. *1*S*PI 
IF(P(ia>.FO.O.) GO TO 380 
IFIXF) Pfl, 1>*-P(im 
SOMTP*Sm 

IFtPt 1,1).IT.P1) GO TO AO 
35 Sin* .5* I S0MTP4S0M2 > 

GO TO 30 

AO IFfPiltl) .GT .P2 ) GO TO 50 
A5 S(ll*.5*IS0Ml#S0MTP) 

GO TO 30 

50 0P*AB$<PX/P2) 

IFIDP.LT.1.E5.AND.DP.GT.1.E-5) GO TO 60 
IFIPf l,l).GT.O.) GO TO 35 
GO TO A 5 

60 X1*S0M2-S0M1 

X2*ISU )-SOMD/xi 
XA*pn f n-pi 
X5*X2*( P2-P1 )-X4 

IF! A8SIX5).GT. 0.001 •ABStXAl ) GO TO 10 
X2*P1/(P1-P2) 

GO TO 20 

10 X3*X2*( X2-1 • )/X5 

D*»P1*X3 

C«.5*l ♦P2*X3*1 • ) 

E*$ORT<C*C-D) 

X2*C-E 

IFIX2.LT.0.) X2*C*E 
20 Sin«SONl*X?*Xl 
30 I COUNT* I COUNT* 1 

IF( ABSI ISm-SOMTP)/SOMTP).LE..O02 ) GO TO 380 

IFf ICOUNT.GT *20) GO TO 370 

IFISOMTP.lT.Sin) GO TO 360 

SOM 2*S0MTP 

P?*Pf 1,1) 


1TFR0010 

ITER0020 

ITER0030 

ITEROOAO 

ITER0050 

ITER0060 

ITER0070 

ITFR0080 

ITER0090 

ITER0100 

ITER0110 

ITER0120 

ITER0130 

itefoiao 

ITER0150 
ITER0160 
ITER0170 
ITER0180 
ITEPOIAO 
ITER0200 
ITER0210 
ITER0220 
ITER0230 
ITER02A0 
ITFP0250 
ITER0260 
I T ER0270 
ITER0280 
ITFR02P0 
I T EP0300 
ITFR0310 
ITEP032C 
ITER0330 
ITER03A0 
ITER0350 
ITER0360 
ITER0370 
I T ER0380 
TTEP0390 
ITEROAOO 
ITER0A10 
HFP0A?0 
ITER0A30 
I TpR0A40 
ITEPOAfO 
ITFR0A60 
ITFP0A70 
ITEPOABr 
ITFROAQO 
ITER0500 
ITER0510 
ITFR0520 
ITEP0530 
ITFR05A0 
JTER0550 
IIFR0560 


raaSRAL p 
gtJQOIl RBAUW 
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it>i. .»t *UiJ »*n |hi\ ( Mti« >*. 

I»* tt»r i r\(i i» IIm ttllr 


GO TO 350 
360 S0H1«S0MTP 
Pl*Pt 1*1) 

GO TO 350 

370 WR1TH6»901) S0M1 t S0M? # $0MTP, PI . P2* P<l*H 

901 FORMAT < 20H CONVERGENCE FAIIURF ♦ 13X *1H1 t ?0X f 1H2 ,20X * 1H3 / 
1 T21 * 3HS0M , 3E20.8 / T22»1HP * 3E20.8 > 

QVRG*»FALSE * 

380 $0M*$(l) 

RETURN 

ENO 


ITER0570 
IT* *0580 
ITER0590 
ITFR0600 
1TER0610 
ITER0620 
ITER0630 
ITER06A0 
ITER0650 
ITER0660 
ITFP0670 
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OS/36© FORTRAN H 

OPTIONS - NAME* MAIN*OPT=02*LINECN7*60.$IZS*OOOOK* 

SOURCE*EBCD IC *NOL IST*NOOECK *LOAD*MAP*NOEDIT* IP* XREF 
SUB ROUTINE MDPLOT (B *K * I M006 *K ENO * F * RP* ) 

01MFNSI0N X ( 63) *v<6 3) *B (3 *11 *TITLF<6 f *WCR0S C 3*3 ) 

OATA X /3*0 .00*3*0.06 *3*0*1 2*3*0.18 *3*0.24* 

% 3*0 . 30*3*0. 36 *3*0.42*3*0.48 *3*0.54* 

$ 3*0.60*3*0.66 t*'*© .72,3*0. 78.3*0.84* 

S 3*0.90*3*0.96*3*1.02*3*1.08*3*1.14,3*1.20 / 

OATA N /63/ * NR/3/ *1 SC ALE/1/ 

DATA TITLE /• •.• /R • * * EV *,3* # •/* 

* WORDS /*COLL • * *ECT I* * *VE •**CVCL»,MC •*• •* 

* *SC IS * * f SOR •** • / 

00 100 I*l*?l 

IT * 3*1 

VCIT-2) « -P(1*I> 

V(IT-l) * ~M2*I> 

YtIT) * -B13.IV10. 

100 CONTINUE 

OMEGA « F/RPM 
CALL CORE !A*4> 

WRITE (6*3) OMEGA 
3 FORMAT IF4.2) 

TITLE ( 1 ) * A 

00 200 •t«l*3 

TITLE ( J*3 ) * WOROSf J*!MO0E) 

200 CONTINUE 

KB * M0D(K t 4) 

IF(KB .EO. 0) KB « 4 

1 PR NT * 0 

IF(KB .EO. 4 .OR. K .€0. KENOI IPRNT * 1 
CALL XYPLOT( Y*X*N .NR* 1 SCALE *KB * IPRNT* TI TL E ) 

RETURN 

END 


ORIGINAL PAGE IS 
OB POOR QUALITY} 


** ' l *« t ‘ I t 


MDPL0010 
MDPL0020 
MDEL0030 
M0PL0040 
MDPL0050 
M0PL0060 
M0PL0070 
MDPL0080 
MrPL0090 
M0PL0100 
WOPLOIIO 
MDPL0120 
MOPL0130 
MOP LO 140 
MDPL0150 
MDPL0160 
MOP LO 170 
MDPL01B0 
MOPL0190 
M0PL020C 
MDPL0210 
MDPL0220 
M0PL0230 
M0PL0240 
MOP LO 250 
M0PL0260 
MDPt0270 
MOP 10280 
MOPL0290 
M0PL0300 
M0PL0310 
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0S/360 FORTRAN H 

OPTIONS - NAME * WAIN *0PT*02 t LINECNT*60.$IZF*OOOOK * 


SOURCE# EBCDIC *NOLIST*NOOECK *10 AD* MAP «NOEO IT* ID*XR EF 



SUBROUTINE M INMAX f X*XMAX *XMIN,N* PR ANG»PRINCH* J*K I 

MNMX0010 


DIMENSION XI 1 ) 

HNMX0020 


XMAX *= xm 

*N«*X0030 


XMIN « XU) 

MN«*X0040 


00 100 I*2*N 

HNHX0050 


IFfXMAX— XU) )3*4*4 

MNMX0060 

3 

XMAX e X(I) 

MNKX0070 


GO TO 100 

MNMX0080 

4 

IM XMIN—Xf I ) )100* 100 * 5 

HNMXC090 

5 

XMIN * XU) 

MNMX0100 

100 

CONTINUE 

MNMX0110 


IFf XMAX ♦ XMIN) 1*1*2 

WNMX0120 

2 

IF(XMAX) 6*1*7 

MNMX0130 

6 

XMAX * 0*0 

MNWX01A0 


GO TO l 

MNMX0150 

7 

XMIN * 0.0 

MWX0160 

1 

CONTINUE 

■WX0170 


XR * XMAX -XMIN 

MNMX018C 


IF(XR .EQ. 0.) XR si. 

HNMX0190 


UPIM * XR / PRANG 

MNMX0200 


JSW « l 

MNMX021C 


IF1UPIM .IT. 1) JSH * 0 

MNf X0220 


J * IFIXULOGIO(UPIM) ) ♦JSH 

MNMX0230 


PWR * 10.**J 

MNNX0240 


PRINCM * * 1 * PWR 

MNMXC250 


X * 1 

WWX026C 


IF! PR INCH .Gf. UPIM) GO TO 10 

HNHX027C' 


PRINCM « .2 *PWR 

MNMX0280 


K « 2 

*N**X0?90 


IFf PR INCH .GE. UPIM) GO TO 10 

MNMX0300 


PRINCM * .5 *PWR 

MNHX0310 


K * 3 

MNMX0320 


IFt PRINCM .GE. UPIM) GO TO 10 

KNHX0330 


PR1NCH « PWR 

MNMX 03*0 


K * 1 

MNMXO3S0 


RETURN 

KNHX0360 

10 

CONTINUE 

WNMX037C 


J • J-l 

MNMX0380 


RETURN 

MNHX0390 


END 

MNHX0400 
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OS/360 FORTRAN » 

OPTIONS - NAME * MA IN »OPT *02, LINECNT*60,SI 2 F*OOOOK , 

SOURCE, EBCDIC ,N0L 1ST, NOOECK , LOAD, MAP, NOEDIT, ID, XREF 
SUBROUTINE PIDUT PL0U0010 

C******** ****************** *********** PLOU0020 

C THIS SUBROUTINE PRODUCES FAN PLOTS * PL0U0030 

C******** **************** ************* PL0U0040 

COMMON /COMA/ JHUB, Nl, DUMB, POUT, IUEU9), NAME I 2) » NDI2 )PLPU0050 

*,NPG,CDATE»2),JHUB1,DBOMI10),RCPLLI10).7I41 ) PIOU0C6O 

*. INPUN MLOUOOTO 

COMMON /COMB/ CK, IRCPt, XRCOLUO), IBOM, RBPMIIO) , PLOUOCBO 

* SHK4U, ZBARI40), EVEBI120), PL0U0C90 

* EYEC 1 120) » SBI40) » SCI40), VMBI40) » VMCI40), VFB140), VFCI40), PLOUOIOO 

* OFBI 40) , DFCI40), TH(4l), THEI40), WTI40), SMI42 ) , 1SOM,XOSPMI200)PIOU0110 


*,A78AR,RPMA,RPMB,RPMC,COLLA,COLLB,COLLC,CHORD 
*»RB 141) ,RC I 41 ) 

COMMON /CPMD/ CMATI5,5),SOMNATl200,3),IPLNt?00,3),INODEI200,31, 

1 MM3,MM4,MM*>,CT|4I ), ST| 41 ) , IB ,IST, 

2 IBS! 10, 10,3), IBB » 10,10, 3), ISTS1 10.3) ,ISTE 110,3) 

COMMON /COMF/ SOMM, TWIST, DI A, SOMI, DELSOM 

COMMON /HINGES/ LCH,LCHP1,LFH,LFHP1 ,CHOFF ,FHOF* ,FCM,FFH, 1PUNCT 
* ,R PM PUN, COL RUN 
•»LPH.lPMPl,PHOFF,FPH 
*»BOMM,TWSM ,LDVN5 


PLPU0120 

*LPU0130 

PL0U0140 

PLPU0150 

PLPU0160 

PL0UC170 

PIPU0180 

PL0U0190 

PLOU0200 

PLPU0210 

PLPU0220 

PLOU0730 

PLOU0240 

PLPU0250 


*»80MI ,TWS1,DELB0M,DELTWS PLOU0730 

COMMON /CPMJ/ HSCFT ,HMASS,VS0FT,VMASS,RS0FT PLOU0240 

1, SPRIG PLPU0250 

COMMON /COMT/ EVXI41), EVBI41), EYCI41), YBI41), YCI41), X)MI |41 ) ,PL0U0?60 
* XITI41 ) , EYRI41), EMRB 141) , EMRCI41), EHRRI41), EMRS0I41), CC02 PLOU0270 
*,QVPLT, OVLIN, SVLIN PLOU0280 

*, BLADES, HUBTYP PL0U0290 

REAL *P SPECIFICATION FOR IBM; COMMENTED FOR CPC VERSION, *LPU0?00 

REAL *8 CMAT PLOU0310 

LOGICAL OPEN, CC02, OVPLT, OVLIN PL0U0320 

DIMENSION IBUFF I 4094) , XMI4), YM)4), XI200), YI200) »LOU0330 

DIMENSION UI10), VI10), JQl(lO), JQ21 10 ) , J03I10) PLOU0340 

DIMENSION ITLEH10) ,ITIE2I9> PLPL'0350 

EQUIVALENCE UTLE <1 ),1TLE1I1) ) *I1TLE?I 1) ,ITLEI 11 ) ) PLOU0360 

DATA OPEN /.FALSE./ PLOU0370 

SUMH1N*FH0FF*CH0FF PLPUC*BO 

IFI SUMH1N.NE .0*0) GO TO 100 PL0U0790 

IFIOVPLT.OR. OVLIN) CALL PLOTI-27.,0. ,-3) PLOU0400 

GO TO 105 PLPU0410 

100 IFIOVPLT .OR. OVLIN) CALL PLOT 1-9. ,0., -3) »L0U04?0 

105 CONTINUE PLPU0430 

IFI OVLIN) GO TO 450 PLPU0440 

IF( OVPLT) GO to 3 PL0U0450 

IFI OPEN) GO TO 2 PL0U0460 

CALL PLOTS 1 1 BUF e » 4098 ) PLPU0*70 

OPEN*. TRUE. PLPU0480 

CALL PLPTI0.5,0.,-3) RLPU0490 

? XMI 1) *0. PLOU0500 

XM<?)*0B0MI IBOM)* 100. PLPUP510 

MAXX*MIM0U5,2*tIFlXIXMI2)*.01*.9)-IFlXtXMt 1 )* .01 ♦ . 01 ) ) ) PLPU0520 

XMAX « 4. PL0U0530 

XMIN*.54IFIX1I8.1“XMAX)*,5) PLPU0540 

XMI 3) *0.0 PL0U0550 

XMI 4) *100. PL0U0560 
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IFIDBOMIIBOM) .GT. 400.)XM|4) « 200. 

IFIDBOMIIBOM) .GT. 800.)XM|4) * 800. 

XMX*XM( 3) ♦XMAX*XM |4 ) 

YM|1)*SPMI 

YM!2)«SPMM4DELS0M 

MAXY*MINOI7,IFIX(YM!2)*.01*.9)-IFIX|VM!1)*.01».01> > 

YMAXsMAXY 
YM| 3) *0.0 
VM(4)*4.*XM!4) 

VMX *YM| 3 )4YMAX*YM 14 ) 

$FcXM 14 )/YM 14) 

3 00 435 1*1,3 

IF< I.NE.3.AND.FH0FF.NE.0) GO TO 435 
IF! I.NE.3.AND.CH0FF.NE.0) GO TO 435 
IF(OVPLT) GO TP 7 
CALL PLOT!-. 5, 0.0.3 ) 

NBLO* BLADES 

CALL PLOT I 0...5.2) 

CALL AX IS IXMIN,1., ’NATURAL FREQUENCY-CPM' ,2l,YM AX,90.,YM| 3) »YM| 4) , 

*10. ) 

IF(I.EO.l) CALL SYMBOL! 3.7, 1.20+YMAX ..1 25. ’COLLECTIVE MODF',0.,151 
IF! I.E0.2) CALL SYMBOL ! 3.9, 1.204YMAX , .125 ,'CYCL TC M0PF',0.,11> 

IF! 1.E0.3) CALL SYMBOL < 3.8 .1.204VMAX ,. 125. »$C I SSORS MODE*. 0., 13) 
CET«3 ,66-.258*FL0AT ( 1 RCPL ) 

CALL SYMBOL I CET , 8.5, .1, ‘ROOT COLLECTIVE *', 0., 17> 

CET-CET41.457 

DO 5 J*l » IRCOL 

CALL NUMBER I CET , 8.5, .1, RCOLLIJ), 0., 1) 

CET »CET4.429 

IFI J.Gf. IRCOL) GO TO 5 

CALL SYMBOL * CET , 8.5, .1* ',', 0., 1) 

CET ■CET*. 086 
CONTINUE 

CALL SYMBOLICET, 8.5, .1, • DEG.' , 0., 5) 

IFI .NOT.CC02) CALL SYMB0H2.36, 8.85, .1 ,3,0., -1) 

IFI .NOT ,CC02)CALL SYMBOL! 2.93,8.80, .1,'T0RSI0N',0.,7) 

CALL NUMBER I 5.5 ,8 .8 ,. 1 .TWIST ,0. ,1 ) 

CALL SYMBOL 16.01,8.8, .l.'OEG TWIST',0.,9) 

CALC SYMBOL 1 7.0.9.0 ,* l , 'V$0 B T*' ,0. ,6) 

CALL NUMBER 17.75,9.0, .1 .VSOFT ,0.0,2 ) 

CALL SYMBOL * 7.0,8 .8 ,. 1, *VMA SS*' ,0. ,6) 

CALL NUMBER 17.75,8.8, .l,VMASS,0.0,2) 

CALL SYMBOL I 7. 0,8. 6, .1, 'MSOFT*' ,0.,6) 

CALL NUMBER 17.75,8.6, .1 ,MSOFT, 0.0,2 ) 

CALL SYMBOL 17. 0,8 .4, .1, *HMA SS*' ,0. ,6) 

CALL NUMBER 17.75,8.4, .1 ,MMA SS ,0.0,2 ) 

CALL SYMBOL 17. 0,8 .2, .1, 'RS0FT**,0. ,6) 

CALL NUMBER I 7.75, 8. 2,. 1,R SOFT, 0.0, 2) 

CALL SYMBOL I 2.36,9.05,.1,2,0.,— 1) 

CALL SYMC-'LI2.68,9.0,.1,»HORI2 PLANE • ,0., 11 ) 

CALL NUMBER I 5. 58, 9. 0, .1,01 A, 0.,1) 

CALL SYMBOL 16.01,9.0, .1, 'FT. DIA',0.,7) 

CALL SYMBOL I 2. 36, 9. 25,. 1,1,0. ,-l) 

CALL SYMBOL! 2.76,9.2, .1,' VERT PLANE • ,0. ,10) 

CALL SYMBOL 12.25,9.4, .1,' SYM MAX AMPLITUDE '. O., 18) 

CALL PLOT! 2 • 25,9.31 5,3) 

CALL PLOT 12 .48,9.375, 2) 

CALL PLOT 1 2. 68 ,9. 37 5, 3) 
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CALI PL0TI3.77,9.375,2> 

CALI P10T(3.02,9.90,3) 

CALL SYMBOL ( 3.02,9.6, .1,ITLE2 ,0.,35) 

CALL SYMBOL ( 2.93, 9.8,.1,ITLE1 , 0.t37) 

CALL PLOT I 0. ,9.5,3) 

CALL PLOT! 0.,10.,2) 

CALL SYMBOL! .5,10. 00,.1»*BHC PROGRAM 0F1758 *,0.,18) 

CALL SYMBOL (2.38, 10.00, .15, 'COUPLED ROTOR NATURAL FREQUENCIES*, 
*0.»33) 

CALL SYMBOL (6. 85, 10.00, .1, NAME ,0.,6) 

CALL SYMBOL 1 7.43, 10. 00,. 1,*( )*,0.,10) 

CALL SYMBOL (7, 516,10. 00,. 1,ND ,0.,8> 

CALL PLOT (8. 5, 10., 3) 

CALL PLOT (8.5, 9. 5,?) 

7 CALL PL0T(XMIN,l.,-3) 

IF(OVPLT) GO r 0 60 

CALL AXIS(0.,0.,»ROTO. RPM* ,— 9,XMAX ,0 ., XM (3 ) «XM( 4) ,10.) 

IN01*0 

IN02*0 

IF(HUBTYP.FC .0.0. AND.NBLD .E0.4) IN0I«2 
IF! HUBTVP.EC .0.0. AND.NBLD .E 0.4) IN02«6 
IF! HUBTYP.EO .0.0. AND.NBLD ,E0. 6) IN01-3 
00 200 1FF«=1,8 
ICMOD«NBLD* ( 2PIFF-1 )/2 
YSPOT«I FF 

IF(l.NE.l) GO TO 210 
C THIS PATH FOR COLLECTIVE MOOES 

IF! MOO! IFF ,NBLO) ,NE .0) GO TO 210 
0ELY«IFF/40. 

DEL X*XM AX/40. 

XL*0.0 

YL*0.0 

00 238 M*1 , 40 

MC«MOO(K *2) -*2 

XL*XL*uELX 

YL«YLtOELY 

CALL PLOT! XL ,YL ,MC ) 

238 CONTINUE 

CALL NUMBER (XMAX,YSPOT-. 05, .1»Y$POT,0.,-1 ) 

CALL SYMBOL! XMAX* .0857.YSPOT-.05, . 1 , */REV • ,0. , 4) 

CALL PLOT 1 0,0,3) 

210 IF! I.NE .2) GO TO 220 
C THIS PATH FOR CYCLIC MOOES 

IF! HUBTYP.EO .1 ) PC TO 220 
IFIMOO! IFF»NBLD).EQ.0> GO TO 220 
IF! IFF.EQ.IN01.0R.IFF.E0.IN02) GO TO 220 
DEL Y« IFF/40. 

0ELX«XMAX/40. 

XL«0.0 

YL*0.0 

00 237 M»l,40 
M"«MOO( M ,2 ) ♦ 2 
XL*XL ♦OELX 
Yl»YL*OELY 
CALL PLOTIXl.Vl.MC) 

237 CONTINUE 

CALL NUMBER (XMAX, Yf'OT— .05, .1 ,YSPOT ,0.*-l ) 

CALL SYMBOL ! XMAX* “C7,YSP0T-.08, . 1 , • /R EV • ,0. ,4 ) 


°l 0U1 i50 
PL0U1 160 
PL0U1I70 
PLOL' • IP V 
PL0U1190 
PL0U1200 
PL0U1210 
PL0U1220 
PL0U1230 
PL0U1240 
PL0U1250 
PL0U1260 
PL0U1270 
PL0U1280 
PL0U1290 
PLOU1300 
PL0U1310 
PL0U1320 
PLOL'1330 
PLOU1340 
PL0U1 350 
PL0U1360 
PL0U1370 
PL0U1380 
PL0U1 390 
PL0U1400 
PL0U1410 
PL0U1420 
PL0U1430 
PL0U1440 
PL0U1450 
PL0U1460 
PLOU1470 
PL0U1480 
PL0U1490 
PL0U1500 
PLOU1510 
PL0U1520 
PL0U1530 
PL0U1540 
PL0U1550 
PL0U1560 
PLOU1570 
PLOU15PO 
PL0U1 490 
PL0U1600 
PL0U1 610 
PL0U1620 
PL0U1630 
PL0U1640 
PLPU1650 
PL0U1660 
PL0U1670 
PL0U1680 
PLOU1690 
PLOl'1700 
PLOU1710 
PL0U1720 
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CALL PLOT !0» 0,3 ) 

220 IFfI.NE.3) GO To 200 
C THIS PATH FOR SCISSOR HOOFS 

IFtHUBTYP.EO.O.O.AND.IFF.EQ.INOl) GO TO 230 
IFfHUBTYP.EO.O.O.AND. IFF.EQ.IN02> GO TO 230 
IFtHUBTYP.EO.O.O) 60 TO 200 

IFtHUBTYP.EO.l.O.AND.MOOtTFF NBLDI.EO.O) GO TO 200 
230 CONTINUE 

0ElVrlPF/40. 

DEL XrXM AX/40. 

XL“0.0 

Vl»0,0 

DO 236 Hr 1,40 

MC«MODtM,2>*2 

Xl*XL*DELX 

VL«YL*DELY 

CALL PLOTtXL »YL »HC) 

236 CONTINUE 

CAL l NUMBER ( XMAX * YS POT-. 05* .1 ,Y$P0T,0.,-1 > 

CALL SVMBOLIXMAS*.085T,YSPOT-.05,.1,'/REV',0.,4) 

CALL' PLOT! 0*0 *3) 

200 CONTINUE 

CALL NUMBER I 0.»-.6, .1 25, BLADES, 0.,1I 
CALL SYMBOL! .5, -.6, .125, 'BLADES', 0., 6) 

IFtHUBTYP.EO.O.O) CALL SYM80Lt2.,-.6».125,'GlMBALE0 HUB', 0. ,121 
IFtHUBTVP.EO.l.ICALL SYMBOL !2., -.6,. 125, 'HINGELFSS HUB' ,0., 13) 
50 CONTINUE 
60 Kl»3 

IFtCC02* Rl«2 
DO 415 K«l,Kl 
K ?*0 

DO 410 IB«1 , IBOM 
DO 410 ISTrl.IRCOL 
J1»IBS1IST,IB,I) 

J2rIBEt I ST, IB, I ) 

IF! J2 .LT.Jl ) GO TO 410 
DO 400 JrJl , J2 

IFtIPIN »J,1).NE.K) GO TO 400 
K2r*?*l 

XtK2)rOBOM| IB) 

Y !K 2) «SOMNAT ! J , 1 ) 

400 CONTINUE 
410 CONTINUE 

IFIK2.E0.0) GO TO 415 
XIK2*l)»XMf 3> 

XtK2*2)*XMt4) 

Y!X2*U»YM! 3 I 
VIK2*2)*YM!4) 

CALL LINEtX,V,X2,l,-l,K> 

415 CONTINUE 
420 CONTINUE 

IF(OVPLT) GO TO 430 
CALL PL0T!8.5-XMIN,-.5,3» 

CALL PLOT!8.5-XMIN,-i.,2) 

430 CALL PL0T(9.-XMIN,-l.,-?> 

435 CONTINUE 

440 CALL TIMEXITU,TT,TL ) 

ET«60.*TT 


PL0U1730 
PLOU1740 
PL0U1750 
PL0U1760 
PLOUI 770 
PL0U1780 
PLOUI 790 
PLOUI POO 
PLOL'IPIO 
PLOUI P20 
PLOUI 830 
PLOUI P40 
PL0UIP50 
PLOUI P60 
PLOUI P70 
P10U18P0 
PLOUI 890 
PLOUI 900 
PLOU1910 
PL0U1920 
PLOUI 930 
PLOUI 940 
PLOL'1950 
PLOU1960 
PLOUI 970 
PLOUI 930 
PLOUI f 
PLO* 

PLOL 

Pl~*7i,20 
Plot .CzO 
PL0U2040 
PL0U2050 
PLOL'2060 
PLOU2070 
PLOU2080 
PL0U2090 
PLOU2100 
PLOU2110 
PL0U2120 
PLOU2130 
PL0U2140 
PL0U2150 
PL0U2160 
PL0U2I70 
PLOU21PO 
PLOU? 190 
PLOU2200 
PL0U2210 
PL0U2720 
PL0U2230 
PL0U7240 
PL0U2250 
PL0U226O 
PL0U7270 
PL0U2280 
PL0U2290 
PLOU? 300 
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WRITE 16 .901) ET 

PLOU2310 

901 

FORMAT 1 38HOPLOT REQUESTED -PROGRAM EXECUTION TIME f FR.2, 

PL0U2320 


1 8H SECONDS ) 

PL0U2330 


return 

PL0U?3*C 

450 

xu >«xMm 

PL0U235C 


OX« XMAX*XM 141/199. 

PL0U2360 


Ytl)*IXtl)-XM13))/XM14) 

PL0U2370 


DO 455 1*2,200 

PL0U23B0 


x»n«xn-iuox 

PL0U2390 

455 

Y1I)*1X1I )-XM13))/XM14| 

PLOU2400 


YR?*YM13)*YMAX*YM4) 

PL0U2410 


Kl*3 

PL0U2420 


IF1CC02) Kl*2 

PL0U2430 


TALC PLOT! XN IN, 1. »— 3) 

PL0U2**0 


DO 630 1*1,3 

PL0U2450 


IF! SU*9UN.NE .O.O.AND. l.NE.J) GO TO 630 

9101*2*60 


00 *60 18*1 »IBOM 

PL0U2470 


JQ11IB)*IBS1 1,18.1) 

PL0U24P0 

460 

J0?tIB)*I8E(l,IB,I) 

PL0U2490 


00 620 K*1,K1 

PLPU?«00 


00 470 IB*1,IB0M 

PI ''U? 510 

470 

J031IB)«J01»IB)-1 

Pi0U2520 

450 

KNT*0. 

PLPU2530 


00 510 JB*1 , IBON 

PL0U2540 

490 

IF1JQ31IB>.6E.J0?1IB)) TO 510 

PLOU2550 


J031IB)«JQ3IIB)*1 

PL0U2560 


JI8 * J03( I B ) 

PL0U2570 


IFIIPIN 1 JIB , I) *NE*K) GO TO 490 

PL0U25B0 


knt*knt«i 

PL0U2590 


V1KNT)«S0MNAT{ JIB ,1) 

PL0U2600 


UIKNT )*DBPM( ?3) 

PL0U2610 

510 

CONTINUE 

PLOU2620 


IF! KNT.LE .1 ) 0 TO 620 

PL0U2630 


C1*0. 

PLOU26*0 


C2*0. 

PLPU2650 


C3*0. 

PLOW? 660 


C**0. 

PL0U2670 


00 520 IB«1,KNT 

PLPU26B0 


C5*U1 IB )**2 

PLOP? 690 


C6*V1 IB )**2 

PLPU2700 


Cl*Cl*C5 

PLPU2710 


C2*C2*C6 

PLPU2720 


C3*C3*C5*C6 

PLPU2730 

520 

C4*C44C5**2 

PL0U?7*0> 


C7*FL0AT1KNT) 

PLOU2750 


C6*IC1*C2-C3*C7)/1C1P*2-€.4*C?> 

PL0U2760 


C5*1C2-C1*C6)/C7 

PLOU2770 


00 530 L*l» 199 

PLPU2780 


C1»SQRT ( AB$t C5»C6*X It )P*2)1 

PL0U2790 


IF1C1 .GE.YMt 31. AND. Cl . IE. VR 2) GO TO 5*0 

PL0U2P00 

530 

CONTINUE 

PLPU2810 


60 TO 570 

PL0U2820 

540 

CALI PLOT f Y 1 L ), (Cl— YM 13 ))/YMI4),3) 

PLPU2P30 


J*L*1 

PLOU28*0 


00 550 l*J*200 

PL0U2P50 


Cl « SORT t AB S1C5*C6*X (L 1**2 ) ) 

PtPU2B60 


IF1C1.LT.YMI3K0R.C1.GT.YR?) GO TO 560 

PL0U2870 

550 

CALL PLOT! Y ( L ) , 1C 1-YM 13 ) )/YM!4) ,2) 

PLPU2R80 
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560 

Cl*fCl-YMm )/VW(*> 

PL0U2890 


CALI SYMBOL? XMAX*.66,C1#. 2 » 

PLOU2900 


CALL NUMBER ( XMAX+.83-C1 ».1»C6»0«»3) 

PL0U2910 


GO TO ABO 

PL0U2920 

570 

IMKNT.EQ.O) GO TO 620 

PtOU?930 


00 580 1*1, KNT 

PtDU2940 

580 

CALL SYMBOL 1 (UU)-XM<3>»/XM|4»,IVlLI-VMm»/YM|4»,.05.K,0..-l» 

PL0U2950 


GO TO 480 

P10U?*60 

620 

CONTINUE 

Ptf>U?970 


CALL PL"T(9.,0.,-3) 

PLOW 980 

630 

CONTINUE 

P10U2990 


CALL PlOTI-XMIN,-l.,-3> 

PL0U3000 


GO TO 440 

PL0U3010 


ENO 

PL0U3020 
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OS/360 FORTRAN M 

OPTIONS - NAHE * PAIN »0PT*02»LINECNT*60»SI?E*0000K» 

SOURCE, EBCDIC »NOL 1ST, NODECK, LOAD, HAP, NOEDTT ,ID*XRE C 
SUBROUTINE START* T0TMAS,FLP1NT,TIPWT,R,LTWS .FIRST) 

C ******** •»#***»***** *•** ************** ****** ******••»••*•**••* 

C CALCULATES COEFFICIENTS WHICH ARE INDEPENDENT OF SWEEPS * 
C******** ************************************ ****************** 

COM HON /COMA/ JHUB. Nl, LOT* POUT, I TIE I 19I»NAHEI2)»NDI2I»NPG 

* tCOATEt 2} » JHUB1 ,DBOH( 10I,RC0LL( 10 >,2141 ), INPUN 
COHMON /COHB/ CK, IRCOL, XRCOLIIO), IBOM, R BOHI 10) , 

* SN2I41), ZBARf 40 ), EYE6I120), 

* EYEC I120> »SBI40) , SCI40), VHBI40), VHCI40), VFM40I, VFCI40), 

* 0FBI40), DFCI40), THI41), THEI40), WTI40), SHl42l,IS0M,XQSnHt20CI 
*, A2BAR.RPNA ,RPNB ,RPHC ,COLLA ,COLLB ,COLLC ,CHORO 
*,RRI41),RC(41) 

COHNON /BIRO/ DUHMY 463) ,M(6 1 , 
l WTPLI21 I, EI8I?01,GICI 201 ,6AI20I,6II 20), 

* TH0I21) 

COHNON /CPMC/ N, IERI7) 

1, OFFSET 

COHNON /COHF/ SOMM, TWIST, DIA, SOHI, OEISOH 
COHNON /COHH/ SPRIP ,FLPSPR 
*,VSOF , VHAS *HSO e .WHAS ,RSOF 
1,S0FI 
*, TORSO 
1,AN6LE,STR 
1,IL0C,TANALE 

COHNON /HINGES/ LCH.LCHPl »t FH,LFHP1 ,CHOFF,FHOFF ,FCH,FFH, IPUNCT 

* ,RPHPUN,COLPUN 
*,LPH,LPHP! ,PHOFF,FPH 
*,BOHN,TWSH ,LDYN5 
*,LOTS 

*,BOHI ,TWSI ,DELBOH,OFLTWS 
COHNON /COHJ/ HSOFT ,MNASS,VS0FT,VMASS,RS0FT 
1,SPRLG 

COHNON /COHT/ EYXI41), EYBI41), EYCI41), VBI41), YCI41), XIM1I41), 

* XITI41), EYRC4I), EHRBI41), EHRCI41), ENRRI41), CHRSOIAl), CC02 

*»OVPLT,OVLIN,SVL1N ,BLAOES,HUBTYP 

LOGICAL LTWS, L2BAR, LEIB, LEIC, LEA, LGAM6, LEI , LGANC, 

♦FIRST, LOT, POUT, l EVP, LEYC, L6AH, CC02 
*,OVPLT,OVLIN,SVLIN 

* .SOUTH 
*,L0YN5 
*,LOTS 
1, OFFSET 

C REAL *R SPECIFICATION FOR IBHt COHMENTFD FOR COC VERSION. 

REAL*8 HNGE.MNLS 
OINENSION CF (40 1 

DATA CVR,CVRPS,HCVH / 0.0174533, 

*0.1047108, 0.00124539 / 

*, SOUTH/. FALSE./ 

1 , HNGE/8HGIHBALED/, HNLS/8M RIGID / 

RPHPUN*BOHI 
COL PUN*TWSI 

IFI BOHI .NE.BONN .AND .TWSI »NE .TWSN I R P NPUN *BOH 1 *DE LBPH 
I FI BOHI .NE . BOHN .AND .TWS I .NE .TWSH ) COLPUN*TwS I *DElTwS 
IFI .N0T.L0VN5) RPMPUN*RPHC 
IF( .N0T.L0YN5) COLPUN*COLLC 
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c 


$TAP0*70 


IFIDEtTWS.lE.O.I DEITWS*10. 

$7680580 


lRC0l*MAX0lltIFIXUTW$N-TWSIl/DFI.TW$4l.0111 

STAP059C 

700 

RC01LU1-TWSI 

STAR060C 


XRC0LI11*TW$I*CVR 

STAR0610 


IFdRCOl.EQ.ll GO TO 720 

STAP0620 


X2*DElTWS*CVR 

STAR063O 


00 710 I*2 t IRC0l 

STAP0640 


RCOLL 1 1 1*RC0U 1 1-1 > *OEt TVS 

SI A RO 650 

710 

XRCOLtl >*XRC0M!-1)«X2 

STAP0660 

720 

BOM I* AM AX1 ( 0 . *BOM 1 ) 

STAR0670 


IFtDELBOM.lE.O.I DELB0M*100. 

STAR0680 


1B0M*MAX0(1.IFIX< fB OMM-BOMT I/OEIBOM+1 .01) ) 

$TAR0690 

730 

1FI 1B0M.LE.10) GO TO 740 

STAR0700 


1ER (71*1 

STAR0710 


IB0H*1 

S7AP0720 

740 

DBOMC11*AHAXII0..60MI1 

STAR0730 


RBOM(11*DBOM(1|*CVRPS 

STAR0740 


IFdBOM.lE.d GO TO 760 

S7AP0750 


X2»DEIB0M*CVRPS 

STAR0760 


00 750 I*2,It>C~ 

S7AP07-*0 


OBOHm*OeOM(l-l»*OELBOM 

S7AR0780 

750 

RBOMI 1 ) -ROOM ( 1-1 ) +% 2 

STAR0790 

760 

IFIR.LE.O.ANO..NOT.LZBAR) IER(2I*1 

STAR080C 


N1-N41 

S7AR0810 


CALL SET1HE ( 550.1 

STAR0820 

60 

CONTI NUF 

STAR0830 


IF! AZ BAR. LE. 0.0 » GO TO 110 

STAP0860 


ZBARI 11*R/FL0ATIN1 

STAR0850 


Zll 1*0.0 

STAR0860 


DO 70 I*1.N 

STARC870 


Z If ♦! )*Z 1 1 ) ♦ AZB AR 

S7AR0880 

70 

ZBAR(I1*ZBAR (11 

STAR0890 

110 

!F( JHUB.EQ.O 1 GO TO 125 

STAROOOO 


DO 120 I*l» JMUB 

S7AP0O10 

120 

THDIIIbO. 

STAR0920 

125 

CONTINUE 

S7AR0930 


FtPlNT«TIPWT*R*R 

S7AR0940 


TOTHAS*TIPWT 

S7AP0950 

80 

00 40 I*1,N 

$TAR0 C 60 


TOTMAS«TOTMAS4WTPU I1*7BAR( 11 

STAR0O70 


FIR INT*FLRINTMWTRC <Il*ZBARtI l*(Z(Il4Z( 1*11 l**21/4. 

S7AR0980 


VMBd l«EI8t I 1 

STARC990 


vMcm»Eicm 

STAP1000 


0UMMV»I)«61Bm/10**6 

STAR1010 


DUMMY (I 4201 *EIC (11/ 10**6 

STAR1020 

90 

CONTINUE 

STAR1030 


00 85 I *1 »N 

STAP1040 


CF( 11*0.0 

STAR1050 


00 87 IJ*I,N 

STAR1060 


CFflKFlIMWTPlf 1J)*ZBARIIJ1*(ZIIJ14Z( Ij4lll*.5 

S7AR1070 

87 

CONTINUE 

STAR 1 080 


CFC I!*ICF(I >4TI PWT*R1 *CVRPS*CVRPS/386.4 

ST API 090 

85 

CONTINUE 

STAR1 100 

C 


STAP1110 


FtPINT«FLPINT/( 32 .2*144.1 

STAPU20 


IFUTWS1 60 TP 180 

STAP1130 


TWO 11 1*0. 

STAR11A0 
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DUM1«TWIST/R STARH50 

D0100 I*1»N STAR116 0 

100 TH>II-H)*TH0U)*ZBAR(1)*0UM1 STAR1170 

180 CONTINUE STAR1180 

GAMMA *0. 002378*5. 73*CH0RD*R**4/IFLPINT912**5) STAR1190 

C**M**M«*M***«** STAR1700 

C PRINT OUT INPUT * STAP1210 

C«»«t********M**** STAR1220 

NP6*NPG*1 STAR1230 

WRITE (6,905) NPG, COATE* NO, NAME , ITLE STAP1240 

005 FORMAT ! 1HI , 27X ,4HPAGE,I3,12X,2'HBMC PROGRAM 0F1T58 -COMPILED • STAR1250 

1 2A4.11X.2A4 /28X.A4.A2.24X, 19 'NATURAL BLADE MOD^S //48X.10A4/49XSTAR1760 

2 *8A4,A3 ) STAR1270 

NSWFEP * IFIXMSOMM-SOMI) /OELSOM) *1 STAR 1280 

1FINSWEEP .LE. 200) GO TO 20 STAR1290 

OELSOM s (SOMM-SOMT )/199. STAR1300 

N SWEEP * 200 ST API 310 

WRITE (6.1 > OELSOM STAR1320 

1 FORMAT I80H0 MOLE THAN 200 POINTS REQUESTED ON FREQUENCY SWEEP. DESTAR1330 


•LTA HAS BEEN CHANGED TO F10.4 1 
20 CONTINUE 

WRI TE 16,904 > 1 1 , ZB AR tl 1 .DUMMY! 1 1 ,OUMMY (1*20), WTPt I 1 1 ,THOI II, 
1CF( 1) ,I*1,N) 

906 FORMAT! 36X. 

173H SEGMENT FI UB-IN**2' WT/.» TWI 

2 CF AT./35X, 

374H LENGTH BEAM ;H0R0 !LB/1NI INB 

4 1NBD FND»/35 X, 

574H CIN) (E-6) IE-61 fO 


;horo 


WT/.-» 


ILB/IN1 


574H IIN) (E-6) IE-61 

6 ILB/RPM) ,/, I 29X, 1 2, F 11.2,46 .4.3,G14.5) 
IFIN.LE.33) GO TO 190 
NPG*NPG*1 

WRITE 16,905) NPG, CDATE , NO, NAME, ITLE 


STAR1349 
STAR 1350 
II), STAR1360 

STAR1370 
STAR1380 
TWIST AT STAR1390 
ST API 400 
INBD EN0 STAR141C 
STAR1420 
10 tv- 1 STAR1430 

STAR1440 
STAR1450 
STAR1460 
STAR! 470 


190 WRITE (6,907) R, BCIMI, BOMM, OELBOM, JHUB, TWSI, TWSM, DFITWS, S0M1STAR1480 


*, SOMM, OELSOM, TH0IN1) STAR1490 

907 FORMAT (1H0,27X,7HRA0IUS*»F7.2,3H IN, T8I.24H INITIAL FINAL 0ELTSTAR1500 

1A / T70.9HR0T0R RPM ,3F9.2 /26X,I2,13H HUB SEGMENTS ,T64, STAR1510 

2 15HR00T COLL 10EG) ,3F9.2 /T*3 ,16HFRE0 SWEEP 1CPM) ,3F9.2 /28X,STAR152C 

3 13HTWIST AT TIP*,F8.3,4H DEG / I $7AR1530 

WRITE ( 6 , 882 ) TIPwT,TORSO,VSOFT,VMASS»HSOFT,HMASS.RSOFT,SPRLG, STAR1540 

1EH0FE,CMCFF, BLADES STAR1550 

882 FORMAT l/29X,9HTIPWEIGHT,6X,G10.3,4H LBM.7X , 13HMAST TOR ST1F , STAR1560 
1 610.3 ,1 1H IN-LB F/ DEG //30X,5MVS0FT,I0X,G10.3,5H /LBF,6X ,5HVMASS,STAR1570 


8X,G10.3,10H LBM/BLAOE //30X, 5HHS0FT,I0X,G10.3,5H /L8F,6X, 


STAR 1580 


3 5WMASS,8X,G10.3, 10M LBM /BLADE //30X,15HFLP SPRIN6/BL0 ,G10.3, STAR1590 

4 24HFT— LBF/OFG INPL SPRG/BLO ,G10.3,UH FT-LBF/DEG //30X, STAR1600 

5 15HFLP HNG OFFSET ,G10.3,5H INCM,6X ,15HINPL HNG OFFSET ,G9.3, STAR1610 

6 5H INCH //30X ,1 5HNUMBER OF BLOS ,G10.3 ) STAR1620 

WRITE (6,873) PHOFF STAR1630 

IFfHUBTYP.EO.l. )WPITE (6,883) HNLS.CHCRO STAR1640 

:FIHUBTVP.NE.1.)WRITE (6,883) HNGE,CH0R0 STAR1650 

STAR1660 

873 FORMAT 1 1H, , 66X »17H PITCH HORN OFFSET ,F10.4,5H INCH ) STAR1670 

LPHePHOFF/ZBARID STAR1680 

FPHe( PHOFF-L PH*ZBAR 1 1 ) )/2BAR ( 1 ) STAR1690 

LPHP1*LPH STAR1700 

LFHbFHOFF/ZBARI 1) STAR1710 

FFH*(FH0FF-LFH*ZBARI1))/ZBAR(1) STAP1720 
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LCH*CHOFF/ZBAR ( 1 ) STAP1730 

FCH*(CH0FF- LCM*ZBAR( 1) )/2BARf 1) STAP1740 

LFMP1*LFH*1 STAR1750 

LCHP1*LCH*1 STAR176P 

HMAS*-MMASS*BLADES/386.4 STAP1770 

VMAS*-VMASS*BLADES/386.4 STAP1780 

1F(HS0FT.E0.0.0)HS0F *1.E20 STAR1790 

IF(V$OFT .EO.O.O IVSO* *1.E20 STAP1B00 

IFf MSOFT.NE .0.000 )HSOF *2.E7/(R*HS0FT ) STAR1B10 

IFIVSOFT.NE. 0.000 )VS0F *2.E7/|R*VS0FT ) STAR1820 

RS0F *RS0FT*12*57.3 STA»1830 

SOFJbS^RLGPI 2*57.3 STAR1840 

883 FORMAT! JM0»29X»15HHUB TYPE ,AR,14X,13HCH0R0 ,F10.3 ,6HSTARIB50 

1 INCHES) STAR? 860 

WRITE (6»802)TDTMAS,FLPINT STAR1870 

802 FORMAT (/29X.10MBLADE MASS,5X ,G10.3,4H L8M,9X,12HFLAP INERTIA, STAR1880 

1 G10.3,17H SLUG-FT**2/BLADE ) STAP1890 

WRITE 16,913) GAMMA STAR1900 

913 FORMAT ( /29X , 17HBLA0E LOCK NUMBER ,010.3 ) STAR1910 

210 IFICC02) GO TO 240 STAR1920 

NPG*NPG*1 STAR1°30 

WRITE <6,905) NPG, COATE, NO, NAME, ITLE STAP1940 

WRITE (6,935) STAR1950 

935 FORMAT <iOX, STAR1960 

1 23MIBB ( IN-LB-SE C**2/IN) 29M BEAM RAO. OF GYRATION UN) , STARl^O 

2 23HICC I IN- LB- SEC* *2 /IN) 29HCH0R0 RAO. OF GYRATION UN) ) STAR1980 

DO 215 1*1, N STAP19O0 

RAD8B * 0. STAR2000 

RAOCC « O. STAR2010 

i c <WTPL!I) .EO. 0.0) GO TO 440 STAR2020 

RADBB * SORT (EYES ( I ) *386.4/WTPLlI) ) $TAR?030 

RAOCC * SORT (EYEC ( I ) *386.4/WTPll I ) ) STAR2040 

440 CONTINUE STAR2050 

WRITE (6,938) 1 ,E YEB( I ) ,RAOBB,EYEC( I ) .RAOCC STAR2060 

930 FORMAT <8X, 1 2 ,6X,E1 5. 5, BX ,G14.4,15X,E15 .5,8X ,G14.4 ) STAR2070 

215 CONTINUE STAR20P0 

9.0 FORMAT ( 1H* ,82X,2F15.2 ) STAR2090 

NPG*NPG*1 STAP2100 

WRITE 16,905) NPG, COATE, NO, NAME, ITLE STAR2110 

WRTTE (6,919) CX STAR2120 

919 ' 1RMAT ( 46X , 25MC0NT ROL SYSTFM STIFFNESS* ,E15.5.6H IN-IB // 33X . STAR2130 

1 l ?W SHEAR CENTER ,8X,11HC.G. OFFSET, 6X ,1?HGJ)LB-1N**2) / 34X, STAR2140 

2 1 IHOFFSET ( IN) ,12X ,4M( IN) /33X,4HBEAM,4X,5HCH0R0,6X,4HBE AM,4X, STAR2150 

3 5HCH0RD ) STAR2I60 

00 230 1*1, N STAR2170 

WRITE (6,931 ) I, SB( I ) , SC(I), RBII), RC(I), GI ( I ) STAR2180 

931 FORMAT! 28X12 ,2F8.3,3X ,2F8 .3 ,3X,E12.5) STAR2190 

230 CONTINUE STAR2200 

WRITE (6,933 ) RB(21) ,RC( 21) STAP2210 

933 FORMAT) 28X.2HTw.19X ,2F8.3) STAP2220 

C****** ***************** ****************** **••**•**• **************** ***STAR2 230 


C CALCULATE COEFFICIENTS INDEPENDENT OF COLLECTIVE ANGLE AND ROTOR RPM STAR2240 
C**»**« ************************** ******************************** *******STAR2250 
240 2(1 )*0. STAR2260 

00 245 1*1 *N1 STAR2270 

245 TM( I)*THO(I )*CVR STAR2280 

DO 300 1*1, N STAP2290 

2(1*1 )*Z(I)*7PARCI) STAR2300 
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IF(I.LT.2*0R.I.GT.19) GO TO 352 
IF IOFFSET) GO TO 352 

tvi*cvmcii-i>4Vhcii41))/ioo. 

TV2*I VMBI 1-1 l+VMB a ♦! ) 1/100. 

IF! VMCI11.LT.TV1.0R.VMB III. IT. TV2 IOFFSET*. TRUE. 

IF ( OF FSET > Z 1 41 ) *F LO AT 1 1 ) 

352 CONTINUE 

IFIVMBII) .GE .5. EO 2) VMBI I) *2 BARI I 1 /VMBI I ) 

IFI VMCI I > *GE .5.E03) VMCI I)=ZBARII> /VMC II 1 
DUM1*0. 

VFBII l*0.5*ZBARm*VMBIll 

DFBII) * ZBAR(I) *0.6666667 *VFBt!) 

VFC II 1*0 . *?*? BAR 1 1 )*VHC 1 1 ) 

OFCII) * 2BARII) *VFCII> *0.666667 
IC0T*ZI41 

TANALF *TANI ANGLE/ 5 7. 3) 

ILOC*ICOT 
ANGLE *0.0 

280 IFICC02) GO TO 290 
MTI 11*1 ./GI 1 1) 

WTII)*2BAR(I)*WTII) 

290 THE II 1*0. 

ifivmbid.eo.0.0) vmbiii«i*o 

IFIVMCm.EO.O 0) V«CII)*1.0 
IF fVMB 1 1 ) • L T.0 .0) VMBII)*-1./IVMBI 1 1*57. 3*12*1 
IFI VMCII1.LT. 0.01 VMC 1 1 1*-1 .0/1 VMC I 1 1*57.3* 12. 1 
IF! I.GT.JHUB) THE 1 1 Is .5*1 TM * t )*TH| 1*1 1 I 
SMI !♦ 11«HCVM*ZBAR II )*WTPLI1 1 

IFI .NOT.CC02) EYXd^llxrSMI 1*11*ZBARI I) **2/12. 
EMRBI M)=0. 

EMRCC I* i ; *0. 

IFISVLIN) GO TO 292 
EMRBI I*l)*SMlT*l)*RBm 
EmCII*l)*SMll#l)*RCII) 

292 EYR I 1*1 1*EHR8II ♦1)*RBII)*EMRCII*1 )*RCII 1 
EMRRI 1*1 1*4 .*SMf I*11*RBII 1*RC< 1 1 
EMRSQ 11*1 >*SH| 1*1 1* IRC 1 I )**?— RBI 1 1**21 
EVB I I *1 )«.5*EYEBI I) *2BAR II) 

300 EYCII*l)*.5*fYECm*ZBARin 
SMI1)*SM|2) 

EYXI1)*EYXI2> 

EYB II )*EYB I 2 ) 

EYCI1)*EYCI?> 

EYRI1)*EYRI2) 

EMRBI 1)*EMRB<2) 

EMRCI1)*EMRCI?) 

EMRRI 1 )*EMRR 12 ) 

BMR $011 1*EMR SOI 2) 

EYR 121 1*EYR I 21 ) ♦TIPWT*I RCI2il**?*RBl2l)**2) ^386 .4 

00 301 1*2, N 

EYX II )*FYX 1 1 )*EYX I I *1 ) 

EYBII)«EYB( I)*EYBtl4l ) 

EYC II)*EYC II)*E YC 1 1 ♦! ) 

EYR II)*EYR 1 1 ) *E YR (I ♦!) 

EMR B 1 1 1 *EMR B f I ) *E MR B I l+ll 
EMRCII)*EMRCm*EMRCl 14 1) 

EMRRI I1«EMRRII)4EMRR| Ul) 

EMR SO 1 1 )*EMRSOI I ) 4EMRS0I I *1 ) 


STAR2310 

STAR2320 

STAR2330 

STAR2340 

STAR2350 

STAR2360 

STAR2370 

STAR2380 

STAR2390 

STAR2400 

STAR2410 

STAR2420 

STAR2430 

STAP2440 

STAR2450 

STAR2460 

STAP2470 

STAR2480 

STAP2490 

STAR2500 

STAR2510 

STAP2520 

STAR2530 

STAP2540 

STAR2550 

STAP2560 

StAR2570 

STAR25P0 

STAP2590 

STAR2600 

STAR2610 

STAP2620 

STAR2630 

STAR2640 

STAR2650 

5TAR2660 

STAR2670 

STAP2680 

STAR2690 

STAR2700 

STAR2710 

S T AR?720 

STAP2730 

STAR2740 

S T AR?750 

STAR2760 

STAR2770 

STAP2780 

STAP2790 

STAR2800 

STAR2810 

STAR2820 

STAP2P30 

STAR2840 

STAP2P50 

STAP2860 

STAP2P70 

STAR28BO 
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301 SMtI!*SM(I)*9MI»lt 

SMf 21»*SM»21 )-»TIPWT/386.4 

00 302 1*1, Ml 

EYR (I )*EYB( 1 l*EYC IT > ♦EVRU) 

YBii»*EYxm*EYcm 

Ycm*fvxm»Eveii) 

X1MH I)*EVCm-EYB» 1) 

302 XlTtl >*-2.4EVxm*EYBm^EYCII» 
ISOM * NSWEEP 

XSOM *S0M1 
SM2IN1>*0. 

J 2^) J 

00 31 5 1*1 tN 
* J 
J*J-1 

315 $MZ t J ) * SMZfJR) ♦SWfJP) ♦ZtJPl 
00 330 1*1* I SOM 
XOSOH U )*<XS0M*CVRPS>9*2 
330 XSOM *XSOM ♦DEC SOM 
FIRST*. FALSE . 

0IA*R/6. 

RETURN 

END 


STAR2B90 

STAP2900 

STAR2910 

STAR2920 

STAR2930 

STAR2940 

STAR2950 

STAR2960 

STAR29T0 

$TAR?9*0 

STAR2990 

STAR3000 

STA»3010 

STAR3020 

STAP3030 

STAR3040 

STAR3050 

STAP3060 

STAR3070 

S7AR30R0 

STAR3090 

STAR3100 

STAR3110 
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OS/ 360 FORTRAN H 


OPTIONS - NAME * MA IN ,0PT*02» LINECNT*60»SIZF*0000K , 

SOURCE .EBCDIC .NOLI ST .NODE CX .LOAD, MAP .NOEDIT ,10, XR FF 
SUBROUTINE SUMMY SUMM0010 

C * ******* **************** ***•*•*•*••**•**•****• '•******••**••** SUMM0020 

C THIS SUBROUTINE PRINTS OUT A SUMMARY OF NATURAL FREOUFNCIES • SUMM0030 

C***************«*********»************************************ SUM HO 0*0 

COMMON /COMA/ JWUBt Nl, LOT, POUT ,1 TLE( 19 > .NAME (2 ) ,N0 (2 ) ,NPG SUMM0050 

*»CDATE( 2) .JMUBI .OBOM| 10 ), RCOLL t 10), 2(41 ), INPUN 6UMM0060 

COMMON /COMB/ CK. I RCOt * XRCOL(IO), IBOM, RBOM(IO) » SUMM0070 

* SM2(41), ZB AR( 40 ) . EVEB(120), SUM HO 080 

* EV EC 1120) ,SB(40) . SC (40), VMBC40), VMCI40), VFBI40I, VFC(40), SUMM0090 

* 0FBI40), DFC(40), TW(4l), THE 1 40), WT(40), SM(42),ISPM,XQSf)MI200)SUNM0100 


*»A2BAR,RPMA,RPMB,RPMC,C0LLA,C0LLB,C0LLC, CHORD SUMM0110 

*,RB 141) ,RC( 41) SUM MO 120 

COMMON /COMD/ CMATt 5,5) , SOM NAT (200, 3 ) ,1PLN( 200,3) .INODE 1 200,3), SUMH0130 

1 MM3 ,MM4,MM5,CT I 41 ) , STI41), IB »IST, SUMM0140 

2 1BS(10,10,3),IBE (10, 10, 3), I STS (10, 3) ,T ST El 10,3) SUMM0150 

COMMON /HINGES/ LCH,LCHP1,LFH,LFHP1 ,CHOFf ,FHf)FF,FCE FFH, IPUNCT SUMM01&0 

* ,R PM PUN, COL PUN SUM MO 170 

*,LPH,LPHP1,PH0FF,FPH SUMM0190 

*,BOMM,TWSM »LDYN5 SUMM0190 

•♦LOTS SUM MO 200 

*,BOMI,TWSI,DElBOM,OELTWS SUMM0210 

COMMON /COM TP/ OEGI 200,3) »PINE (?,4) ,ODES( 2, 5 ) SUMM0220 

REAL *8 SPECIFICATION FOR IBM; COMMENTED FOR COC VERSION. SUMM0230 

REAL *8 CMAT SUMM0240 

REAL *8 PINE, ODES SUMM0250 

DIMENSION S0MNA1 (200) * S0MNA2(200), SOMNA3(200), SUMM0260 

•IPK200), IP2(200), IP3(200), IN0K200), )N02(200), 1N03(200), SUMM0270 

* 181(10.10), 183(10,10), 186(10,10). 182(10,10), SUMM0280 

*184(10,10), IB6 < 1 0, 10 ) * ISTl(lO), 1ST3(10), IST5(10), IST2(10), SUMM0290 
•1ST 4(10), IST6(10),ITLE1(10),ITLE2(9) SUMM0300 

* ,DEG1(200), 0FG2( 200) , 0EG3(200) SUMM031O 

EQUIVALENCE ( 1TLE (l ) , ITLE1 ( 1 ) ) ,(ITLE2( 1),ITIE(11)) SUMM03V0 

EQUIVALENCE ( SOMNAT (1 ,1 ) .SOMNAl ( 1 ) ) ,( S0MNAT1 1 ,2 ) ,S0MNA2 ( 1 ) ) , SUMM0330 

1 (S0MNAT(1,3),S0MNA3(1)),(IPLN(1,1),1P1(1)) , SUM MO 340 

2 ( IPLN ( 1 ,2 ) , IP2( 1) ) , (IPLNI 1,3) ,IP3(1 )) , SUMM0350 

3 (INPOff 1,1),1N01(1)) , ( INODE (1,2) , IN02 ( l ) ) , SUMM03&0 

4 ( INODE ( 1,3) ,IN03( 1)),(IBS(1 ,1 ,1 ), IB 1(1,1)), SUMM0370 

5 (IBS(1,1,2),IB3(1,1)), (IBS(1, 1,3), 186(1.1)), SUMM0380 

6 ( I B F ( 1 ♦ 1 , 1), 1 62(1,1)), (IBE(1, 1,2), 184(1,))), SUMM0390 

7 (I8E( 1.1.3) ,166(1,1)), (ISTS(1,1),IST1(D) , SUMM0400 

8 (ISTS(1,2),IST3(1)) , ( ISTS(1 ,3) ,IST5tl) ) , SUMM0410 

9 (I5 T E (1 ,1 ),IST2(1)) , (IS T E(l,2),IST4d)) , SUMM0420 

A ( 1STE (1,3 )» IST6 ( 1 ) ) SUMM0430 

EQUIVALENCE (DEGt 1,1) ,0EG1( 1)) , (OEGd ,2),DEG2(1)> , SUMM0440 

1 (DEG( 1,3) »0EG 3(1)) SUMM0450 

IF( FHO* : F.NE .O.OR .CHOFF .NE .0 ) GO TO 3 SUMM0460 

NPG*NPG*1 SUMM0470 

WRITE (6,901) NPG,CDATE.N0,NAME.ITLE1.NAME,ITLE1.ITLE2,ITIE2 SUMM0480 

901 FORMAT (1H1 ,27X,4HPAGE, 13,12*, 29HBMC PROGRAM 0F1768 -COMPILEO , SUMM04O0 

2A4.11X.2A4 /57X,19HNATURAL BLADE MODES //IX, 2( °Y ,A4,A2 ,4X,0A4, SUMM0600 
A 1, 10X ) /1X,2(?0X,8A4»A3,11X) //19X,30HC OLLECTIVE M 0SUMM0510 
0 E .40X.22HC v C L I C MODE/ 2( 5X ,20HNATURAL ROOT ROTOR ,SUMM0520 
4X, 19HMAX 1 mum NUMBER OF ,5X ,3HMAX , 10X ) / SUMM0S30 

2 (6X,4HFPEQ,4X,4HC0LL,3X,3HRPM,4X,9MAMPLITUPE,4X,5MN0DES,4X, SUM MO 640 

10HDEFLECTI0N ,6X) / 2( 5X ,4H/REV, 5X ,3HDE6 ,33X , 10H ANGLE — OEG »6X) 1SUMM0650 


LINES*0 


SUMM0660 
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IN1*1 

IN2*1 

00 60 I$T*1 *IRCOL 

IFt 1ST21 1ST > .LT^ISTlf 1ST) .AND. I STM 1ST) *LT* IST3IIST) ) 

* GO TO 60 
UNE$«LINES+1 
WRITE <6 *90?) 

902 FORMAT 1 1H ) 

00 SO IB*leIBOM 

IF(IB2CIST,ie).LT.IBl(IST»IBI.AN0.IB6CIST,IB).LT.I83(IST,IB)l 

* GO TO 50 
WRITE (6 *902 1 
L1NES*LINES*1 

MB*IB2( IST»1B)-IB1( 1ST* IB 1^1 
MC*IB6(IST,IB)-IB3( IS t ,IB)M 
MA*MAXOIMB,MC) 

IF(LINES.EQ.?.0R.LINES*MA.LE.50) GO TO 20 
NPG*NPG*1 

WRITE 16,901) NPG, CDATE, ND , NAME, ITLE 1, NAME, 1UE1,ITLE?,1TLF2 
WRITE 16, 903) 

903 FORMAT) 1H0) 

LINES*2 

20 LINE$*L INES^MA 
DO *0 1*1, HA 

WRITE 16, 902) 

IF( I.GT.MR) GO TO 30 

SOMNA 1 1 INI )* SOMNA KIND /DBOM( IB ) 

IF * IFl(INl) 

IN * INOK INI) 

WRITE 16,904) S0MNA1IIN1), RCOLLIIST), OSOMdB), (PLNE)J»IP), 
•J*l ,2 ), (ODFS(K,IN) , K*l,2), OEGl(TNl) 

SOMNA 1 ( INI ) * SOMNA K IN 1) *OBOM( IB ) 

906 FORMAT (1H*,F10.5,2F7.1,2X,A8,A3»2X,A8,A2,F7.1) 

INI *IN1+1 

30 1FII.GT.MC) GO TO 60 

SOMNA?) IN2)*S0MNA?t IN?) /OBOMI IB) 

IP * IP2IIN2) 

IN * IN02IIN2) 

WRITE (6,905) S0MNA2IIN2), RCOiLIIST), DBOM(IB), tPlNEIJ,TP), 
•J“l ,2 ), (PDF S(K ,1 N) , K«l,2), 0EG2IIN?) 

S0MNA2I IN?)*SOMNA?( 1N2)*0B0M(1B) 

90S FORMAT (1H«,F76.5,?F7.1,2X,A8,A3,2X,A8,A2»F7.1) 

IN2*IN2*i 
60 CONTINUE 
50 CONTINUE 
60 CONTINUE 
3 CONTINUE 

WRITE (6,907) NPG,CDATE,ND,NAME,1TLE1, ITLE2 

907 FORMAT (1H1 , 27X ,6HPAGE, 13 ,1 2X,29MBHC PROGRAM DF1758 -COMPILED , 

1 2A6,11X*2A6 /57X,19HNATURAL BLADE MODES //10X,A6, A2,6X,9A6, A1 

2 ?1X,8A6,A3//19X,17H SCISSOR S,6X,7HM ODE/ 5X, 

3 2 OH NATURAL ROOT ROTOR ,6X,19HMAXIMUM NUMBER 0F,5X,3MMAX / 
6 6X,6HFREO,6X,6HCOLL,3X,3HRPM,6X,9HAMPIITVDE,6X,5HNODES,6X, 

5 10HDEFLECTI0N /5X,6H/REV,5X,3H0EG,33X,10HANGLF — DEG ) 

LINES*0 

DO 160 1ST*1,IRC0L 

IF I IST6( 1ST ) .LT .1 ST 51 1ST) ) GOTO 160 
LINE$*L1NES*1 


SUMM0570 
SUMM05B0 
SUM MO 590 
SUM MO 600 
SUMM0610 
SUMM0620 
SUMM0630 
SUM MO 660 
SUMM0650 
SUMM0660 
SUM MO 670 
SUMM06P0 
SUMM0690 
SUM MO 700 
SUM MO 710 
SUM MO 720 
SUM MO 730 
SUM MO 760 
SUMM0750 
SUMM0760 
SUM MO 770 
SUMM07B0 
SUMM0790 
SUM MO 800 
SUMM0810 
SUMM0E20 
SUM MO 830 
SUM MO 860 
SUMM086O 
SUM MO 860 
SUM MO 870 
SUM MO 880 
SUM MO 890 
SUM MO 900 
SUMM0910 
SUM MO 920 
SUMM0930 
SUM MO 960 
SUMM0960 
SUM MO 960 
SUM MO 970 
SUM MO 980 
SUM MO 990 
SUM Ml 00 C 

suMMioie 
SUMM1020 
SUMM1030 
SUM Ml 060 
SUMM1050 
/ SUMM1060 
SUMM1070 
SUMM1080 
SUM Ml 090 
SUMM1 100 
SUMM1H0 
SUMMH20 
SUMM1130 
SUM Ml 160 
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WRITE 16 *902» 

DO 150 

lF(IB6tIST,JB).lT.IB5lIST,ieil 60 TO 150 
WRITE (6*002) 

UNES-tINES*l 

WA*IB61 1ST »IB1— TB5f 1ST, 18’ 

IFILINES.E0.2.0R.I.INES4MA.LE.501 60 TO 120 
NR6«NPG*1 

WRITE (6,00?) NPC. COATE, ND,NAME,1TLE1, IUF? 

WRITE (6 ,9031 
UNFS«2 

120 LINES*CINES*HA 
HA«IB5UST,IB) 

WB*IB6( 1ST, T8 I 
00 145 I*MA,MB 
WRITE 16,9021 

S0MNA3U 1*S0MNA3(I I/DBOMIIBI 
IP * IP3( T 1 
IN • INP3( I ) 

140 WRITE16.9C4) S0MNA3U1, RCOLLUST), DBOMtIB), (PINMJ.1P1, 

*J«1,2>» (PDESOUINl ,K «1,21, DFG3U) 

145 S0MNA3I I )rS0MNA3(I 1*DB0MI1B1 
150 CONTINUE 
160 CONTINUE 

CALI TIMEXC TU»TT»TL 1 
ET*60.*TT 

WRITE (6,9061 ET 

906 FORMAT I <>H0RIW TIME ,F8.2 ,8W SECONDS 1 
RFTURN 
END 


SUWM1I5C 
SUM Ml 160 
SUMM1170 
SUMM11R0 
SUMM1 190 

sum mi 200 

SUM Ml ?10 
SUMW1220 
SUMM1??0 
SUMM1240 
StfMM125P 
Sum mi 260 
SUMM127C 
SUMM1 280 
SUM»»1290 
SUMM1300 
SUMM1310 
SUMMI320 
SUM;*1330 
TUMM1340 
$UMM13*0 
SUMM1 360 
SUMM1370 
SUMM1380 
SUMM1390 
SUMM1400 
SUMM1410 
SUMM1420 
SUMMJ430 
SUM-1440 
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0 S/360 FORTRAN H 


OPTIONS - NAME * MAIN,0PT^02,L1NECNT*60,$1ZF*0000K, 

SPURCE»EBCD IC »NOL 1ST, NODECK, LOAD, MAP,NOEO!T, 10, XREF 
SUBROUTINE XYPLOT ( V,X ,N,NR ,ISC ALE ,M, IPRNT, TITLE » 

X , Y INPUT ARRAYS TO CROSS PLOT 
N NUMBFR OF POINTS TO PLOT 

NR ALPHABET CYCLE KEY: CYCLES THRU FIRST NR LETTERS 

ISCALE « 0 MEANS PLOT X AND Y ON SAME SCALE 

NOT * 0 MEANS USE BEST SCALE FOR EACH VARIABLE 
M BLOCK NUMBER TO BE FILLED BY THIS CALL I A BLOCKS TOTAL I 

IPRNT * 1 MEANS TO PRINT THE WHOLE PAGE OF 4 PLOTS 

NOT * 1 MEANS TO RETURN WITHOUT PRINTING 
TITLE 24 ALPHAMERIC CHAR ACTFPS DESCRIBING PLOT VARIABLES 
L0GICAL»1 TITLE,HEAD.NSCL,NEXP,NSIGN, BLOCK, ALP, BLANK, ONE, MINUS 
DIMENSION TITLE! 1) .HEAD! 36 1 ,NSCL 1 3 ) ,NE XP II 01 ,NS!GN( 2» 

DIMENSION BLOCK I 60 * 131 ),ALP! 27 ) 

DIMENSION XI 1 ) »Y ! 1) ,SCALE!2),IX7!4),IVZ!4) 

DATA HEAD /• *,» •, *S • , »C • , »A », »L • , *E # , »S • , • t •, • •» , V»,*E*, 

1 •R»,»T*.».»,2*» ♦ ,».»,»E»,4*» »,*H»,»0*,»R», •I»,»2%».»* 

2 2 *» *, •.»,•£* , 2 *» • / 

DATA NSCL /M','2V5V, NEXP /»0» , • 1 • , *2 * * • 3* , »4* , • 5 • » »6 • , »7» , 

1 • 8*, *9 » /.NSIGN /•♦»,»-»/ 

C DATA ALP/*A , ,’B*, , C , , , 0 , ,*E*,*F*,*G , ,*m*,»I»,»j , ,*k , , , L , , , M*, 

DATA AL P/*B •, , C , »»T».»0 , ,*E** , F*,»G*,»M»,»I*,»J*,*K , ,»L , » , M», 

* •N*,*O»,*P»,»0 , . , R*,»S*,*T»,»U , .»V*, , W*,*X»,*Y*,*Z*,* • / 

* ,MINUS /•-•/, ONE /•»•/ 

* ,IY7/ 2,34,2,34/ ,1X7/ 1, 1,69,6** / 

IF! IPRNT ,LT, 0) GO TO 40 
I C !M.NE.1> GO TO 30 

BLANK * ALP I 27) 

DO 300 I *1,60 
DO 300 J *1,131 
300 BLOCK (I, J) * BLANK 
30 CONTINUE 

CALL MINMAX (X,XMAX,XMIN,N,6.1,SX,JX,KX ) 

CALL MINMAX ( v , YMAX,YMIN,N»4. 33333, SY,JY,KY) 

IF! ISCALE ,NE • 0) GO TO 20 
IFISY .GT. SX) GO TO 10 
SY * SX 
JY * JX 
KY * KX 
GO TO 20 
10 SX * SY 
JX * JY 
KX * KV 
20 CONTINUE 

SCALEX * 10. / SX 
SCALEY * 6. / SY 

I ZERO « -XMIN *SC AL EX ♦O.* 

J7ER0 « -YMIN * SC AL EY *0.5 
IZERO * I2ER0 * IX 2 ( m ) 

JZERO « J7ERP *IY7(M) 

DO 100 K*1 ,62 
K1XZ * K»IXZ(M) 

BLOCK !J2ER0,KIXZ> * MINUS 
IF1K.GT.27) GO TO 100 
KIY7 « K*IYZ(M) -2 
BLOCK IKIY7, 1 ZERO) « ONE 
100 CONTINUE 


XVPL0010 

XYPL0020 

XVPL0030 

XYPL0040 

XVPL0050 

XYPL0060 

XYPL0070 

XYPL0080 

XYPL009C 

XYPL0100 

XVPL0110 

XYPL0120 

XYP10130 

XYPL0140 

XYPL0150 

XYPL0160 

XYPL0170 

XVPL01P0 

XYPL0190 

XYPL0200 

XYPL0210 

XYPL0220 

XYPL0230 

XYPL0240 

XVPL0250 

XYPL0260 

XYPL0270 

XYPL0280 

XYPL0290 

XYPL0300 

XYPL0310 

XYPL0320 

XYPL0330 

XVPL0340 

XYPL0350 

XYPL036C 

XVPL037C 

XYPL03P0 

XYPL0390 

XYPL0400 

XYPL0410 

XYPL0420 

XYPL0430 

XYPL0440 

XYPL0450 

XYPL0460 

XYPL0470 

XYPL04B0 

XYPL0490 

XYPIOKOO 

XYPLOBIO 

XYPL0S20 

XYPL0B30 

XYPL0S40 

XYPL0550 

XVPLOBfcO 
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l « 1 

00 200 K«1,N 

1 * I2ER0 4XIK)*SCALEX ♦0.5 
J * JZERO ♦VIK) *SCALEY 40.5 
BL0CKIJ.lt * ALPIL) 

l « l»l 

ifil.gt.nr) l*i 

200 CONTINUE 

C KILL IN NUMBERS IN HE A01N6 
HEAD! IB) * NSCL1KY) 

HEAOI 32 ) = NSCL(KX) 

NS * I 

IFIJV .LT, 0) NS *2 
HEAOI 21 > * NSIGN<NS> 

NS * 1 

IFtJX .LT. 0) NS *2 
HEAOI 35) * NSIGNINS) 

NS <*»ABSIJV> +1 
HEAOI 22) » NEXPINS) 

NS *IABSIJX) 41 
HEAOI 36) - NEXPINS) 

C PUT HEADING INTO BLOCK 
J « IVZIM>-1 
IS * IX2IM) +1 
DO 400 K«l,24 
BLOCKIJ.IS^K I * TITLE IX) 

400 CONTINUE 

IS « IS 424 

DO 500 K«l,36 

BLOCK I J. IS^K ) *> HEAOI K) 

500 CONTINUE 

40 1FI IPRNT ,E0. 0) RETURN 

WRITE 16,1) 1 1 BLOCK I J, I) ,1*1,130 

1 FORMAT l*l*/l» » I30A1 ) ) 

WRITE 16.2) 

2 format i»o b and c are scaled to l 

RETURN 

END 


XYPL0570 

XYPL05B0 

XYPL0590 

XYPL0600 

XYPL0610 

XYPL0620 

XYPL0630 

XVPL0E40 

XYPL0650 

XVPL0660 

XYPL0670 

XYPL06B0 

XYPL0690 

XYPL0700 

XYPL0710 

XVPL0720 

XYPL0730 

XYPL0740 

XVPL0750 

XYPL0760 

XVPL0770 

XYPL0780 

XYPL0740 

XYPLOPOO 

XYPLOBIO 

XVPL0P20 

XYPL0P30 

XVPL0840 

XYPL0B50 

XVPL0B60 

XVPL0R70 

XYPLOBBO 

t J* 1 » 60 ) XYPL0P90 

XYPLOOOO 
XVP 10910 

INCH? T IS SCALED TO 10 DF6.» )XYPL0«20 

XVPL0P30 

XYPL0440 
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